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ABSTRACT 


An echosounder is used to probe various atmospheric 
parameters. An acoustic wave is transmitted into the 
atmosphere and information deduced from the backscattered 
energy. 
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excess attenuation, are modeled to predict the performance 
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I. INTRODUCTION 


Many atmospheric turbulence-dependent processes take 
place in the troposphere, the lowest ten to twenty 
kilometers above the earth's surface. Acoustic energy 
interacts with irregularities in the atmosphere more 
strongly than electromagnetic waves and is potentially a 
better probe for these irregularities. Acoustic echoes have 
been used to detect: 

1) wind speed and direction profiles, 

2) HUMiGweEy apTorl les, 

3) strength and location of temperature inversions, 
4) temperature inhomogeneities, 

5) mechanical turbulence. 

The acoustic sounder, also known by the names 
echosounder, echosonde, sodar, and acoustic radar, transmits 
a pulse of acoustic energy into the atmosphere. The various 
atmospheric parameters can be determined based on the 
intensity and frequency of the scattered energy. The range 
to the scattering center is determined from the elapsed time 
between transmission and return of the scattered energy. 

In this thesis we are attempting to improve our 
understanding of the fundamental range limitations of the 


acoustic sounder. In general, the acoustic sounder suffers 


4 


from multiple scattering. Turbulence at shorter ranges 
alters the phase front of the propagating pulse and reduces 
the magnitude of the return from longer ranges. 

We are also exploring the use of return from shorter 
ranges, along with theory for the degradation of the energy 
along the path, to compensate for short range degradation in 
a boot-strap fashion and quantify the return from longer 
ranges. Acoustic sounders at present indicate the presence 
of inhomogeneities at a specific altitude but the magnitude 
of the inhomogeneities may have errors of a factor of four 
Or more. 

The following factors affect the range of the 
echosounder: 

1) Atmospheric: 

a) Pressure, 

b) Temperature, 

c) Temperature Structure Parameter ce 

dad) Velocity Structure Parameter Cv? ’ 

e) Water-vapor Pressure, and 

£) Ambient Noise. 

2) Echosounder: 

a) Antenna Aperture Factor, 

b) Antenna Diameter, 

Go, ELFLcCIeney , 


d) Frequency, 


e) Power Transmitted, and 


f) Pulse Length. 


II. BACKGROUND 


Acoustic energy propagates in the atmosphere as a 
longitudinal wave of pressure variations. The energy is 
scattered, attenuated, and refracted. Energy which is 
scattered constitutes the returned signal. The attenuation 
decreases the energy ensonifying a given volume and 
decreases the returned signal. Energy that is scattered 
through small angles or that is refracted results in further 


degradation or excess attenuation of the signal. 


A. SCATTERING 

The echosounder transmits a pulse of acoustic energy 
that is scattered by temperature and velocity 
inhomogeneities. Information about the inhomogeneities is 
then based on elapsed time between transmission and receipt 
of the return signal, the strength of the returned signal, 
the equations for scattering, and the Doppler shift. 
Variations in the propagation velocity of the wave, which 
are a function of the temperature and velocity variations in 
the air mass, produce the scattering. In addition to 
scattering, the temperature variations cause refraction 
while the velocity turbulence causes both a shift in phase 


and direction of propagation [Ref. l:p. 60]. 


Random temperature or wind structure caused by 
turbulence and uniform gradients of temperature or velocity 
contribute to the scattering. The gradient must change in 
a scale size that is comparable to or smaller than the 
acoustic wavelength in order to contribute to the 
scattering. It appears that the effects of uniform 
gradients are limited to beam bending and possible specular 
reflection for the acoustic frequencies of interest here. 
This refractive structure of the troposphere and 
stratosphere would be better probed with what has come to be 
called infrasound, sound of frequencies below twenty Hertz. 
[Ref «ols pewoe) 


The turbulent atmospheric temperature and velocity 


fluctuations, at high wave number, follow the Kolmogorov 


79/3 


power spectral density law (in one dimension). The 
Fluctuations are correlated spatially on the order of one 
centimeter to tens of meters. The nonzero correlation 
lengths and the declining power spectral density imply a 
different scattering process than for random point 
scatterers, even though the turbulence is a stochastic 
process [Ref. l:p 61]. 


An expression for the power scattered from a unit volume 


per unit incident flux into a unit solid angle is [Ref. 2:p 
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CQ =.055 cos*e (cv* cos?(®) + .13 ct?) (sin(2))711/3 
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where 
is the wavelength of the transmitted sound, 
T is the average atmospheric temperature, 
O is the angle of scattering from the direction of 
propagation, 


Cv? is the velocity structure parameter, 


and Ct? is the temperature structure parameter, 


Ct?= <Tia) > Pi xtr) >? 
1g 


This scattering equation indicates that backscatter 
(6=77) is only a function of the temperature structure. The 
velocity at any point is the vector sum of the phase 
velocity of the sound and the particle velocity of the 
turbulence. This vector sum is always in the forward 
direction because the particle velocity is always less than 
the speed of sound c. Therefore the backscatter is only a 
function of the temperature structure [Ref. l:p 60]. The 
vertical turbulent velocity does cause a Doppler snift of 
the frequency of the backscattered radiation. 

Scattering over a region with correlated scattering 
centers produces constructive and destructive interference. 
The backscattered waves are partially coherent. This 
results in much greater intensities than would be received 


from incoherent scattering such as Rayleigh scattering [Ref. 


l:p 61]. Also the scattering with the interference, over 
the nonzero Conerenee length, acts like an array of 
scattering centers in a regular crystal lattice of spacing 
Leo The Bragg condition 

L= : /sin(®/2), 
is satisfied. The dominant scattering is for scale sizes of 


\/2. [Re£f. l:p 61] 


B. ATTENUATION 
The atmosphere absorbs some acoustic energy that 
propagates through the atmosphere and reradiates this energy 
at different frequencies. Quantitatively the power lost 
eq Over a path length 1 is given by 
P= ae 
where is the attenuation coefficient in Nepers per meter. 
Historically the attenuation coefficient has been considered 
to be the sum of several terms; 
A ei classical viscous losses, 
xe molecular rotation losses, 
eer molecular vibrational losses, No and O° 
Classical and rotational loses are negligible below 
about three kilohertz, the region of interest for 


echosounder operation [Ref. l:p 54 and Ref. 3:p. 18-2], 


consequently vibration of Nj and O, produce most of the 


attenuation. 


Rotational and vibrational losses are referred to as 
relaxation processes. The acoustic energy excites internal 
energy modes of the N4 and O05 molecules. The rate of 
collisions with water vapor determines the rate and 
efficiency of conversion of the energy into translational 
energy (heat). 

The phase is shifted due to the relaxation processes. 
This is one of the reasons some of the microwave radar pulse 
compression techniques cannot be used in echosounders. A 
number of the pulse compression techniques rely on the phase 
not changing during propagation. 

The dependence of the attenuation on the water vapor 
pressure is believed to be due to a resonance process 
Deeween the lewest vibrational stateswof the O» and N9 
molecules with the water molecules. Henderson and Hertfeld 


in Reference 4 [p.986] state that the lowest vibrational 


states of O05, and water vapor are only thirty-nine 
wavenumbers (569K) apart at 1600 cm7!, For this reason 
oxygen was thought to be primarily responsible for the 
humidity dependent absorption of sound. Henderson and 
Herzfeld [Ref. 4], and many others, assumed nitrogen was an 
inactive dilutent Stee no effect on the relaxation 
processes [Ref. 4:p. 987]. Unfortunately the theory did not 


agree with the data for low frequencies and high 


humidities. 


Theory incorporating the relaxation processes involving 
nitrogen and water at low frequencies and relative 
humidities greater than twenty-five percent bring the theory 
into agreement with the data. At high frequencies and low 
relative humidities, theory and data match well with oxygen 
making the main contribution and nitrogen acting as an inert 
dilutent. At low frequencies and high humidities, nitrogen 
seems to make the main contribution and oxygen acts as an 
inert dilutent. For the range of relative humidities in 
which many people live and over a good part of the audible 
Frequency range it appears nitrogen is the major contributer 
to the relaxation processes [Ref. 5:p. 165]. This is in 
contrast to the previous assumption that oxygen was 
responsible [Ref. 6:p. 604]. 

An expression for attenuation due to the molecular 


absorption is given by [Ref. 7:p. 34], 


Oo = Snax (ors+e «(ze ecemy) ” wre 
304.8 fm, 1 +(£/7 im 


in a@B/m 
where 
f is the frequency of the transmitted pulse, 


is the Napier frequency, 


f, =(10+6600 h+4400 h) ( P ) _ _ Sie ) ere 
1014 }) \1.8T+492 


x 
? 
i 


=.0078 Eq( 1-8 r+ag2y 2-5 Seo Tho 1 (a a9 ) 
1,.8T+492 


ines 


where 
T is temperature in degrees Celcius, 


h is the percent mole ratio of water vapor, 


h= 100 (e/p), 
e= water vapor pressure in mb, 
P= atmospheric pressure in mb, 


To convert from dB/m to Nepers/m note that 
10 log(I(x)/I(0))= 10 log(exp(-«x), 
==> 4.34 <(Nepers/m)= X (dB/m). 

The Napier frequency is the frequency of maximum 
absorption per wavelength and “,,, iS the attenuation at 
the Napier frequency. The Napier frequency is shifted to 
higher frequencies by even small amounts of water vapor. 
Plots of attenuation versus water-vapor pressure in 
millibars (Figure 1) and versus relative humidity (Figure 2) 
are shown for a range of frequencies. 

For a given temperature, the attenuation has a maximum 
and decreases for higher or lower humidities. With higher 
temperatures, the maximum attenuation increases and the 
relative humidity at which that maximum occurs decreases. 
Plots of attenuation versus water-vapor pressure in 
millibars (Figure 3) and versus relative humidity (Figure 4) 


are shown for a range of temperatures. 
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Fig. 1. Attenuation versus Water-vapor Pressure (mb.) 
for Frequencies around 1.6 KHz. 
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Fig. 2. Attenuation versus 
Relative Humidity ($%) 
for Frequencies around 1.6 KHz. 
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Fig. 3. Attenuation versus Water-vapor Pressure (mb.) 
for Temperatures Around 20 Celsius. 
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Fig. 4. Attenuation versus Relative Humidity (3%) 
for Temperatures around 20 Celsius. 
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C. EXCESS ATTENUATION 

The atmosphere scatters and attenuates the acoustic 
energy. In the previous section on scattering, we were only 
concerned with the energy scattered from a given volume. 
Excess attenuation accounts for scattering and energy loss 
for the round trip in the atmosphere to and from a given 
volume of interest. Excess attenuation arises because the 
atmosphere degrades the mutual coherence of the acoustic 
wave. The divergent solid angle of the acoustic wave is 
larger than would occur for a coherent, diffraction limited 
wave. 

The excess attenuation, Ze, accounts for this energy 
lost due to small-angle scattering. This turbulent beam 
broadening reduces the on axis intensity. Clifford and 
Brown in Reference 8 [p. 1972] develop the equation 

Ze = 1/(1+N) fOr ING 

ze =1.5/(1+N) for N>1l, 
where 

N= (i 

Do= antenna diameter, 


/ 0e is the atmospheric acoustic coherence length 


= 1.46*k2 ‘o\ ds*cng*(s)*((1 s)°/>+( om aks ac 
0 


aR. aR 


The term varies from a value of one for no excess 
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attenuation to an asymptote of zero, implying the energy 
would be spread over a 2 solid angle. 

There is a step in the functional dependence of the 
excess attenuation when N=l1 that can be seen in the 
equations above and in the plots. There clearly is no 
physical discontinuity but rather a transition between the 
theoretical dependence between two asymptotic regions. 

Figures 5,6, and 7 are plots of excess attenuation 
versus range for antenna sizes of .5, 1, and 1.5 meters, 
respectively, for several frequencies. As can be seen from 
these plots or the equations above, the relationship of the 
antenna size to the coherence length has a significant 
effect on the excess attenuation. A larger antenna will 
not, by itself, produce greater range. The antenna size 
must be matched to the transmission frequency in terms of 
antenna design guidelines and effects upon the excess 


attenuation. 
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Fig. 5. Range versus Excess Attenuation 
for an Antenna Diameter of .5 m. 
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Fig. 6. Range versus Excess Attenuation 
for an Antenna Diameter of 1.0 m. 
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Fig. 7. Range versus Excess Attenuation 
for an Antenna Diameter of 1.5 m. 
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Itt EXPERIMENT 


The acoustic sounder is an excellent tool for probing 
the lower troposnhere. It suffers from several 
shortcomings: 

1) The range has typically been limited to a few hundred 
meters. 

2) It has proven difficult to quantify the measurements 
accurately for a given range. 

In this thesis, we approached the problem from several 
directions. We looked at the hardware to understand the 
Signal to noise limitations and to understand how the 
magnitude of the returned signal could be calibrated. We 
developed a software model which would allow one to estimate 
the power backscattered from the atmosphere at a given 
range, based on profiles of atmospheric characteristics and 


input parameters for the acoustic radar. 


A. HARDWARE 

The acoustic sounder we were working with was the 
Aerovironment Model 300. It consists of an electronic 
module which generates a 1600 Hertz electrical pulse. The 
pulse is converted to acoustic energy by a transducer which 


feeds a 1.25 meter parabolic reflector. Energy 


Ce 


backscattered by the atmosphere is then received by the 
reflector and transducer. The electrical signal is then 
Filtered and amplified. In addition, a ramp amplifier 
compensates for the 1/r decreasing signal amplitude with 
range to decrease the dynamic range requirements of the 
present data display, a strip chart recorder. 

We replaced the various integrated circuit amplifiers 
and filters with more current designs with lower noise. We 
replaced the pre-amplifier with an OPA 111 and the rest with 
LF 356 BN devices. 

We traced the amplifier and filters of the receiver 
board to understand the undocumented choices the 
manufacturer had made. Figures 8 and 9 show a pre- 
amplifier, high and low pass filters and two stage 
amplification. The filters are bi-quad configured with 
notched outputs. Figure 10 shows the notching. Figure 10 
represents the frequency spectrum output of the receiver 
board as measured with the HP 3561A Frequency Spectrum 
Analyser with random noise from the HP 3561A providing the 
input signal before the filters. Figure 1] cepresents 


frequency spectrum of the receiver board with the random 


noise across the transducer. A 10> ohm resistor was used to 


match impedences as shown in Figure 12. 





Fig. 8. Schematic of the Receiver Board 
of the Aerovironment Model 300 
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Fig. 9. Schematic of the High and Low Pass Filters of 
the Receiver Board of Aerovironment Model 300. 
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Fig. 10. Frequency Spectrum of the Receiver Board of 
the Aerovironment Model 300. (Input before filters) 
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Fig. ll. Frequency Spectrum of the Receiver Board of 
the Aerovironment Model 300. (Input across transducer) 
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Fig. 12. Connection of Noise Source to Transducer 
The manufacturer had a filter-oscillator board after the 
receiver board with a switch for three pass bands; narrow, 
medium, and wide. With the same input as in Figure 10 we 
measured the frequency output for the three settings of the 
Filter~oscillator board (Figure 13). As can be seen, there 
is little difference between the three settings. 

From previous measurements, the resistor labeled Ros in 
the schematic (Figure 8) of the receiver board was found to 
load down the input. The resistor was initially 25 kilo- 
ohms. Figures 14, 15, and 16 represent the signal from the 


input transformer with R =25K, 100K, and && ohms 


25 
respectively. Random noise from the HP 3561A was input 
across the transducer as shown in Figure 12. The Q of the 
input transformer was improved by increasing the resistance. 
The lower curve in each Figure is with no noise signal input 
maom the HP 3561A. 

Tt appeared the Q of the filters could also be increased 
by adding the resistors labeled R, in Figure 17. The upper 


curves in Figures 18, 19, and 20 are the output of the 


receiver board with each Sys ae 750, and 560 ohms 
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Fig. 13. Frequency Spectrum of the Filter Oscillator Board 


of the Aerovironment Model 300. 
(Input before the filters of the receiver board) 
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Fig. 14. Frequency Spectrum of the Input Transformer 


(Input across transducer with Ro5_25, 
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Fig. 15. Frequency Spectrum of the Input Transformer 
(input across transducer with Rj5-100k ) 
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Fig. 16. Frequency Spectrum of the Input Transformer 
(input across transducer with R25= o¢) 
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Tig. 17. Schematic of the 
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Fig. 18. Frequency Spectrum of the Receiver Board with Ry= oc 
(input before filters) 
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Fig. 19. Frequency Spectrum of the Receiver Board with 
Ro=/50 - (input before filters) 
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Fig. 20. Frequency Spectrum of the Receiver Board with 
Ry=260 < (input before filters) 
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respectively. The input was random noise from the HP 3561A 
after the pre-amplifier but before the filters of the 
receiver board. The effect of decreasing the resistance 
waS an increase in gain but the pass-band was no longer flat 
on top. However, for the range of vertical wind velocities, 
about two or three meters per second, a pass-band of fifty 
Or sixty Hertz is adequate to allow for the Doppler shift. 
As can be seen, the filters could be tuned for a narrow 
pass-~-band which would have better shape. We would gain 
amplification and still have a pass-band that is wide enough 
to meet our needs. 

The lower curves in Figures 18, 19, and 20 are with 
notching, the upper without. The effect of the notching was 
a Sharper cut off near the bandpass frequency (1600 Hz) but 


a loss at higher and lower frequencies. 


B. SOFTWARE 

The signal received by an acoustic sounder from a given 
range tells us that scattering centers exist but we can say 
little about the size and magnitude of the scattering 
centers. We don't know how much energy was incident on the 
Scattering volume nor how much the signal was degraded on 
the return path. I£, on the other hand, previous signals 


are used to determine the atmospheric characteristics then 


ae 
fl 
4 


an estimate of the degradation of the incident power and the 
degradation of the return signal can be made. 

The software programs I have written here go in the 
opposite direction. Given certain atmospheric 
characteristics, the power returned is estimated. With the 
modules tested, it would then be a matter of turning it 
around to take actual data and estimate atmospheric 
conditions. In the present form they serve our purposes in 
allowing us to explore the effects of parameter changes. 

All the programs are fundamentally built around the 
echosonde equation, also sometimes referred to as the radar 
equation in meteorology [Ref. 9 and Ref. 7:p. 3]. 


Pr = Br*(PE*HhE)* (expl-ceak ee OD) Ok eee )*(A*G) *Ze*, 
2 R 


where 


Pr is th 


iL 


power returned from a range R, 

Pt is the power transmitted at frequency f, 

Er is the efficiency of conversion of acoustic power 
to electrical power by the transducer , 

Et is the efficiency of conversion of electrical power 
to acoustic by the transducer, 


c is the speed of sound in m/sec, 


f is the pulse length, 

A 1s the area of the antenna, 

R is the range, and 

G 1s the effective-aperture factor of the antenna. 
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Transducer efficiencies must be measured for each driver 
and typically range from a few percent [Ref. 7:p. 3] to 
twenty five percent [Ref. 10:p. TI-10]}. 

The scattering cross section per unit volume (0 9) is 
the fraction of incident power backscattered per unit 
distance into a unit solid angle at a given frequency. From 
Reference 7 [p. 4} and Reference 2, 

ame 4.0080 * iet73e% cts 
0 
where 
k is the wavenumber = 2 T/wavelength, 
To is the local ets temperature in Kelvin, 
Ct’ is the temperature structure parameter, 
The power scattered from a scattering volume is 
Pb(1)=(Pt*Et-Pb(I-1))*exp(-2R) *(cV/2) *Ze7* J. , 
The power returned to the antenna is 
Pr=Pb*A*G*Er/R4, 
where the return path attenuation was already included in 
Ps. The excess attenuation Ze* was discussed in the 
background section. 

The dependence with height of the velocity structure 
parameter Cv? was needed for the calculation of the excess 
attenuation. Reference 2 and Reference 11 [p.149] give 


Cy2 = 2 rE 2/3, 


where € is the dissipation rate of turbulent kinetic energy. 


On 
cy 


Reference 11 [p. 154] and Reference (p 194) give 


G&G = (.33m/s)2 41 oe Oe R3/5 3/2, 
k * R 


for a stable surface layer. Figure 21 is a plot of range 
versus this Cv%. 

The programs have four temperature structure parameter 
e) profiles. The operator must choose one. The first is 
based on data as presented in Reference 13 [p. 398] for 
midday clear weather above the Tularosa Basin desert in New 
Mexico. The second choice is for the same data multiplied 
by a factor two to approximate looking up a convective 
plume. The third case for a nocturnal atmosphere assumes a 
dependence proportional to the negative exponential of the 
range as presented in Reference 13 [p. 399] with tower data 
from the same reference used for the first sixty-five 
meters. The fourth case assumes a dependence proportional 
to height to the -4/3 and a surface vertical heat flux of 
.095 [Ref. 7:p. 7]. Also case four allows for the operator 
to input the height of an inversion layer with Ge being 
proportional to height to the -4/3 above the inversion 
layer. Figures 22, 23, 24, and 25 are plots of these four 
profiles. 
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CV aad cr2 were then used to calculate the acoustic 


refractive index factor [Ref. 14:p. 119] 
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Fig. 21. Velocity Structure Parameter (cm) Profile. 
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22. Range versus Temperature Structure Parameter CS) 
for Profile l 
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Fig. 23. Range versus Temperature Structure Parameter Ge) 


for Profile 2 
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Fig. 24. Range versus Temperature Structure Parameter (o>) 
for Profile 3 
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Fig. 25. Range versus Temperature Structure Parameter (Cn) 
for Profile 4 
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Figure 26 is a plot ver ne 

Figure 27 represents an input flow chart. It summarizes 
the effect of each variable or atmospheric parameter. 
Figures 28 and 29 are plots of the range versus the power 
returned to the echosounder for profiles 3 and 4. 

The segments of the programs were tested against 
existing data to verify proper operation. Based on input 
temperature, atmospheric pressure, water-vapor pressure, and 
the frequency of operation of the acoustic sounder the 
programs calculate the attenuation. If the operator 
desires, the first program will plot the attenuation as a 
Function of absolute water-vapor pressure and/or relative 
nhunidity for frequencies at one-third octaves around the 
input frequency and then again for temperatures at ten 
degree Celcius intervals around the input temperature. 

These plots were used to check the attenuation against data 
(Ref. 15] and [Ref. 5]. 

Ambient noise levels of acoustic sounders are found to 
be about ten to forty dB above the theoretical Johnson noise 
limit [Ref. 16:p. 19-4}. This noise level determines the 
maximum range. This maximum range 1S comparable to range 
capabilities of the Aerovironment System 300 when the 
Operating parameters of the Aerovironment are used in the 


program. 
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Fig. 26. Acoustic Index Structure Parameter Cr Profile 
using Cp* Profile 3 
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Fig. 28. Range versus Power Returned to the Antenna 
for Cp* Profile 3, Parameters of the Aerovironment 
Model 300, and S.T.P. 
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Fig. 29. Range versus Power Returned to the Antenna 
for Cr Profile 4, Parameters of the Aerovironment 
Model 300, and S.T.P. 
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Trends, as input parameters where changed, were also 
used to check the program output. As the frequency 
increased, the backscatter of the turbulence in the 
atmosphere decreased slowly and the absorption of the 
atmosphere increased rapidly in agreement with Reference 17 
[p. I-10]. 

The calculation for excess attenuation was difficult to 
check. At ranges of about 40090 to 500 meters the excess 
attenuation of a typical acoustic radar has been found to be 
about .25, in agreement with the program calculation. 
PROGRAM 1 

The first program takes AE WOM MOIST jaeieeNaee and the 
echosounder parameters as inout and outputs nine plots. The 
inputs, by the operator at the keyboard, are: 

Atmospheric parameters; 

Atmospheric pressure in millibars, 

The profile of ct? From four options, 

Temperature in degrees Celsius, 

Water-vapor pressure in millthars, 

Ecnosounder data; 

Antenna diameter in meters, 

Frequency of the echosounder lmederez, 

Power transmitted by the echosounder in Watts, and 


Pulse length of the transmitted acoustic energy. 
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The program outputs the following plots: 


1. Attenuation (1/m) versus water-vapor pressure (mb) for 


five frequencies at one-third octaves around the input 


frequency. 


Z. Attenuation (1/m) versus relative humidity (%) for five 


Frequencies at one-third octaves around the input 


Frequency. 


3. Attenuation (l/m) versus water-vapor pressure (mb) for 


five temperatures at ten degree intervals around the 


input temperature. 


4. Attenuation (l1/m) versus relative humidity (%) for five 


temperatures at ten degree intervals around the inpdut 


temperature. 


5. Range (m) 
6. Range (m) 
7. Range (mn) 
8. Range (m) 

Senuc ture 


9. Range (m) 


versus 


versus 
versus 
versus 
factor 


versus 


echosounder. 


excess attenuation. 


the temperature structure factor Cee. 


thenvelociiymestructune. Lactor cve, 


the acoustic refractive index 
Cne-. 


the power backscattered to the 


All of the programs prompt tiie operator for inputs and 


with a series of yes/no questions allows the operator to re- 


run with the same inputs or change the inputs. Figure 30 1s 


a flow chart of program one and appendix 1 is a copy. 
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Fig. 30. Flowchart of Computer Model One 


46 


Peegures 1, 2, 3,08, 21, 22, 23, 24, 25, 26, 28,0mnd 29 are 
Gutput plots of program one. 

The second program takes the same atmospheric parameters 
used in the first program. With the exception of the 
frequency, it takes the same echosounder parameters also. 
This program outputs a plot of range versus frequency and a 
plot of range versus excess attenuation for various 
frequencies. Figure 31 is a flow chart of the program and 
appendix 2 is a copy. Some output plots of the program are 
included in the conclusions section. 

The third program takes the same atmospheric and 
echosounder Sara ceks as the first. The program outnduts a 
plot of range as a function of efficiency of the transducer, 
assuming the same efficiency for transmit and receives. A} 
the other programs and plots in this tnesis assume 


ictencies of 25%, which is on the high side of typical 


(es 
he 


ma 


NY 


performance. Appendix 3 is a copy and an output plot is 
presented in the conclusions section. 

The fourth program has the same inputs as the first wit’ 
the exception of antenna area, which is the dependent 
variable for the output plot. It plots the range as a 
function of antenna area for several frequencies and is 
presented and discussed in the conclusions. The model does 
not include different efficiencies based on optimum antenna 


design for different frequencies. The output reflects the 
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Fig. 31. Flowchart of Computer Model Two. 
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effect of different antenna diameters on the excess 
attenuation and the echosonde equation. I am uncertain how 
much these drive the considerations for optimum antenna 


design. Appendix 4 is a copy of the program. 


The fifth program has the power to the transducer as the 
dependent variable plotted against range for several 
frequencies. Appendix 5 is a copy and the output plot is 
presented and discussed in the conclusions section. 

The sixth program has the ambient, atmospheric 
background noise as the dependent variable. Appendix 6 
contains a copy of the program and the output plot is 
presented and discussed in the conclusions section. 

These programs, taken together, allow a parametric 


analysis of the effect of different Darameter changes. 


a 


IV CONCLUSIONS 


Reference 10 [{p. II-l1] indicates that as the frequency 


increases: 


i) 


3) 


4) 


Background and wind noises decrease except for marked 
peaks due to fans, etc. This relationship between 
noise and frequency is not included in the this model. 
The reflectivity of the turbulence in the atmosphere 
decreases Slowly and the absorption of the atmosphere 
increases rapidly. The model agrees with this for low 
values of the water-vapor pressure. The effects of 
water-vapor pressure on the attenuation are shown in 
Figures (ito 4. 

Wildlife sounds tend to increase with increasing 
Frequency and dominate the background noise at about 
3000 Hertz. This effect is not included in the model. 
The Doppler shift requires the receiver bandwidth to be 
increased. The model doesn't take this into 

account. The bandwidth is used to calculate the 
Johnson noise but Johnson noise is not significant. 
The noise really should be scaled with frequency to 


depict the gain in range accurately for lower 


frequencies. 


cy 
O 


Transducer efficiencies vary with frequency. This is 
not included in the model. An increase in efficiency will 
have a much greater impact on the potential range than 
increases in, say, power. An increase in efficiency will 
help both with transmission and return, increasing the 
transmitted power and the returned electrical signal 
strength. Figures 32 is a plot of the maximum range for 
efficiencies of .05 to .5. This figure is the output of 
program 3. The range increases quickly with improvements in 
efficiency. The efficiency could be improved by using 
better designed horns, such as catenoidal or exponential. 

Optimum antenna diameters vary with frequency. The 
Output of program 4, Figure 33, demonstrates this effect. 
Not included are the effects on the antenna effective 
aperture factor G, except to the extent it may be effected 
by the excess attenuation. The discontinuities in the 
Curves are due to the discontinuities in the equations for 
the excess attenuation. For a given frequency, the range 
increases and then, as antenna diameters increase beyond an 
optimum, the range decreases. The decrease for larger 
antenna diameters is the result of the excess attenuation. 

One other consideration as to choice of frequency is 
that the resolution is increased with increasing frequency. 
A specific need for the echosounder may drive this 


Gude auot. ands mot considered in this model. 
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Fig. 33. Range versus Antenna Diameter 
for cre Profile 3, Parameters of the Aerovironment 
Model 300, and S.T.P. 
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Figure 34 is a plot of range as a function of frequency 
for an antenna size of 1.5 meters. The water-vapor pressure 
input was ten millibars. Figure 35 is for the same area trioe 
size but with a water-vapor pressure of 2 millibars. The 
step in each curve is not a physical effect but rather due 
to the step of 1.5 in the theoretical equation for the 
excess attenuation. Figure 1 showed the greater value of 
the attenuation at around 2 mb. water-vapor pressure for the 
range of frequencies we are dealing with. These vrevious 
two plots demonstrate the effect. Figures 34 and 35 are 
some of the output plots of program 2. Both figures 
demonstrate the greater attenuation of the acoustic wave as 
the frequency increases. 

Significant increases in range can be achieved by 
Judicious choice of frequency. Considerations are the 
water-vapor pressure (Figures 1 to 4), the excess 
attenuation (Figures 5 to 7), which depends on the antenna 
diameter (Figure 33), the frequency (Figure 33 to 35), and 
the anticipated frequency spectrum of the backyround noise. 

Figure 35 models the effect of increasing the 
transimitted power to increase the range. This is the output 
plot of program 5. The slope of the curve is not very steep 
and becomes less so for higher frequencies. A considerable 
increase in power is required to double the range. Also the 


model does not include nonlinear transducer effects which 
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Fig. 34. Range versus Echosounder Frequency 
for ae Profile 3, Parameters of the Aerovironment 
Model 300, and S.T.P. 
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Fig. 35. Range versus Echosounder Frequency 
for Cr Profile 3, Parameters of the Aerovironment 
Model 300, and S.T.P. except a Water-vapor Pressure of 2 mb. 
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arise when the power is increased. Reference 1 [p. 57] 
points out that at some point the acoustic wave becomes 
distorted, which implies, from Fourier analysis, there is a 
flow of energy out of the fundamental frequency into the 
higher harmonics. Since the higher frequencies are 
attenuated more quickly we will soon reach a point of 
saturation, which Brown [Ref. 1 :p 57] called nonlinear 
saturation to distinguish it from another saturation effect 
he discusses. 

Figure 37 is the output of program 6. This shows the 
effect of decreasing the ambient background noise level. 
This decrease in noise might be achieved by using digital 
processing and fast Fourier transforms to achieve a narrower 
bandwidth. The bandwidth, and therefore the noise, can be 
much smaller. 

These piots demonstrate which parameter changes might 
best improve the range of an echosounder, which was one of 
Our goals in this thesis. The modeled performance changes 
should allow for intelligent decisions of the necessary 
parameters for the expected uses of the echosoundec. 

We also sought to explore the ability to use the 
return#1 signal to quantify atmospheric parameters 
accurately at a given range. As can be seen from Figure 27 
Gunemsenout flow chart), for a single acoustic radar return 


you could calculate en based on an assumed profile of Cy2. 
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Fig. 36. Range versus Power Transmitted 
for Cp Profile 3, Parameters of the Aerovironment 
Model 300, and S.T.P. 
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Fig. 37. Range versus Noise Level 


for Cin? Profile 3, Parameters of the Aerovironment 
Model 300, and S.T.P. 
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However this need to assume a Cv? profile could be 
eliminated with the use of two or more echosounders. 
Reference 11 and many others describe techniques. 

Another approach would be to measure the Doppler width. 
The Doppler width is the spread of frequency around the 
Doppler shifted frequency. Epsilon and therefore Cv? can de 
related to the Doppler width. With the use of fast Fourier 


GIS eae aoe, Vie 


transforms and digital processing Cv 
measured simultaneously with one echosounderc. 

Using the returned signal from lesser ranges the energy 
incident on a given volume could be estimated and the 
degradation of the return signal could be estimated. in 
this way the computer program would allow one to essentially 


boot-strap up to a given range and more accurately depict 


the atmospheric Darameters based on the returned signal. 
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APPENDIX A 
COMPUTER PROGRAM 1 


!FULLER, ROBERT PROG 1 
| 

11@ SEP 85 

| 

| 

[aeeee ee PURPOSE tH ee ean 


THIS PROGRAM WILL TAKE PARAMETERS OF AN ECHOSOUNDER AND 


!ESTIMATE THE RANGE. THE FOLLOWING INPUTS ARE REQUIRED: 


| 
| 
| 
| 
| 
| 
| 
! 
| 
! 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
! 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
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ATMOSPHERIC DATA 
1 AATOMOSPHERIC PRESSURE in millibars 
Z)SELECT ONE OF FOUR ATMOSPHERIC PROFILES FOR THE TEMPERATURE 
STRUCTURE PROFILE. 
a)FOR ONE PROFILE THE HEIGHT OF THE INVERSION LAYER IS ALSO INPUT 
S) TEMPERATURE IN GEGREES CEB SrUs 
4)WATER VAPOR PRESSURE IN millibars 


ECHOSOUNDER DATA 

S)ANTENNA DIAMETER IN METERS 

6 FREQUENCY OF ECHOSOUNDER IN Hz 

7)POWER TRANSMITTED BY ECHOSOUNDER IN WATTS 
8)PULSE LENGTH OF THE TRANSMITTED ACOUSTIC ENERGY 


THE PROGRAM OUTPUTS THE FOLLOWING PLOTS: 
1 YATTENUATION( 1/m.) VERSUS WATER-VAPOR PRESSURE(mb) FOR FIVE 
FREQUENCIES AT ONE-THIRD OCTIVES AROUND THE INPUT FREQUENCY. 
Z)ATTENUATION(C 1/m.) VERSUS RELATIVE HUMIDITY(%) FOR FIVE FREQUENCIES 
AT ONE-THIRD OCTIVES AROUND THE INPUT FREQUENCY. 
SY ATTENUATION( 1/m.) VERSUS WATER-VAPOR PROSSURE( mb) FOR FIVE 
TEMPERATURES AT TEN DEGREE INTERVALS AROUND THE INPUT TEMPERATURE. 
4)ATTENUATION(C 1/m. ) VERSUS RELATIVE HUMIDITY(%) FOR FIVE TEMPERATURES 
AT TEN DEGREE INTERVALS AROUND THE INPUT TEMPERATURE. 
S)RANGE(m.) VERSUS EXCESS ATTENUATION 
6B )RANGE(m.) VERSUS TEMPERATURE STRUCTURE FACTOR. 
7)RANGE(m.) VERSUS VELOCITY STRUCTURE FACTOR. 
B)RANGE(m.) VERSUS ACOUSTIC STRUCTURE FACTOR. 
SJ)RANGE(m.) VERSUS BACKSCATTERED TO ECHOSOUNDER. 
H###242# VARIABLES 44 44% 
Again TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 


WISHES TO MAKE ANOUTHER RUN. 


Ano_change TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 


WISHES TO CHANGE THE VALUE OF INPUT DATA BEFORE MAKING 
ANOUTHER RUN. 


Ant_area ECHOSOUNDER ANTENNA AREA CALCULATED FROM INPUT OF 
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Ant_diam 
Atom_pres 
Atten 
Att_freq 
Att_max 
Att_temp 
Bn 

C 

Cl .Ge€Z 
es 
CneZ(*) 
CteZ(*) 
CveZ(*) 
Dx 

Dy 
Epsilon 
Et 

Er 

Es 
Exc_att(*) 
F 

Fmax 
Freq 
Freq_con 
G 


ANTENNA DIAMETER. 

INPUT OF ECHOSOUNDER ANTENNA DIAMETER IN METERS. 

INPUT OF ATMOSPHERIC PRESSURE IN mb. 

ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 
SUBPROGRAM ATTENUATION. 

VARIABLE USED TO DETERMINE IF OPERATOR 
WANTS TO PLOT ATTENUATION VERSUS WATER_VAPOR 
PRESSURE FOR VARIOUS FREQUENCIES AROUND THE INPUT 
FREQUENCY. IF SO THEN PLOT IS DONE IN SUBPROGRAM 
Att_freq. 

VARIABLE IN SUBPROGRAM ATTENUATION. IT IS THE 
ATTENUATION AT THE FREQUENCY OF THE MAXIMUM 
ATTENUATION FOR THE INPUT CONDITIONS. 

VARIABLE USED TO DETERMINE IF OPERATOR WANTS 
TO PLOT ATTENUATION VERSUS WATER_VAPOR PRESSURE 
FOR VARIOUS TEMPERATURE AROUND THE INPUT TEMPERATURE. 
THE PLOT IS DONE IN THE SUBPROGRAM Att_temp. 

BANDWIDTH OF RECEIVER NEEDED TO RECEIVE SIGNALS 
DOPPLER SHIFFED BY THREE METER PER SECOND 
VERTICAL VELOCITIES. 

VELOCITY OF SOUND CALCULATION AT ZERO CELSIUS. 

VARIABLES USED TO SCALE INVERSION HEIGHT IN 
TEMPERATURE STRUCTURE PROFILE THREE: 

CONSTANT USED IN CALCULATION FOR Cte2 PROFILE 4. 

ARRAY OF VALUES OF THE ACOUSTIC REFRACTIVE INDEX 
PARAMETER. CALCULATED BASED ON SELECTION OF 
PROFILE FOR TEMPERATURE STRUCTURE PARAMETER 
AND ASSUMED PROFILE FOR THE VELOCITY STRUCTURE 
PARAMETER. 

ARRAY OF VALUES OF THE TEMPERATURE STRUCTURE PARAMETER. 
THE OPERATOR SELECTS A PROFILE FROM SEVERAL PROVIDED. 

ARRAY OF VALUES OF THE VELOCITY STRUCTURE PARAMETER. 
VALUES BASED ON CALCULATION USING ASSUMED DISSAPATION 
RATE. 

STEP SIZE FOR X AXIS FOR’ VARIOUS PLOTS. 

STEP SIZE FOR Y WxAIS FOR VARIGWS PLOTS. 

DISSAPATION RATE USED IN CALCULATION OF CveZ. 

TRANSMISSION EFFICIENCY OF ECHOSOUNDER. 

EFFICIENCY OF ECHOSOUNDER WHEN RECEIVING ACOUSTIC 
BACKSCATTER. 

SATURATION VAPOR PRESSURE AT GIVEN TEMPERATURE. 

EXCESS “ATTENUATION” AT GIVEN RANGE. 

VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 

FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
ATTENUATION. 

INPUT FREQUENCY OF ECHOSOUNDER. 

USED TO INSURE VALUE OF FREQUENCY PASSED TO SUBPROGRAM 
Att_freq WAS NOT CHANGED. 

ANTENNA EFFECTIVE APERATURE FACTOR. 
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H 
I 


Interval 
Inver 

J 

K 

L 
Mess_up 


N 
New_va 


Noise 





VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 
MAIN INDEX FOR LOOP TO CHECK RETURNED POWER UNTIL 
POWER WAS LESS THAN MINIMUM DETECTABLE SIGNAL. 
MAXIMUM EFFECTIVE-SCATTERING VOLUME THICKNESS. 
HEIGHT OF INVERSION LAYER. 
FIRST ORDER INDEX FOR ASSORTED LOOPS. 
WAVENUMBER 
SECOND ORDER INDEX FOR ASSORTED LOOPS. 
TEST VALUE TO CHECK RESPONSES WHEN OPERATOR IS INPUTTING 
RESPONSES. 
TEST VALUE IN CALCULATION OF COHERENCE LENTH IN CALCULATION 
OF EXCESS ATTENUATION. 
TEST VALUE IN CHECKING IF OPERATOR WISHES TO CHANGE 
A VARIABLE BEFORE A NEW RUN. 
ASSUMED MINIMUM DETECTABLE SIGNAL. 


Pow_back(#) POWER BACKSCATTERED FROM GIVEN RANGE. 
Pow _ret(#) POWER BACKSCATTERED TO ECHOSOUNDER FROM GIVEN RANGE. 


Pow_trans 


Profile 
Pstar 


Pulse 
R 


Ri 
Range( # ) 


Remainder 


Roe 
Rho 


Sigma( * ) 


INPUT OF POWER SUPPLIED TO TRANSDUCER OF ECHOSOUNDER. 

OPERATOR INPUT OF Cte2 PROFILE FROM AVAILABLE PROFILES. 

VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. | 

TRANSMITTED PULSE LENGTH IN MILLISECONDS. 
OPERATOR. 

RANGE FOR INPUT TO INTEGRAL OF EXCESS ATTENUATION. 

VARIABLE USED IN SCALING INVERSION HEIGHT FOR PROFILE THREE. 

ARRAY OF RANGE VALUES. 

REMAINDER OF MODULO FUNCTION USED TO DECREES THE 
NUMBER OF PASSES THROUGH THE INTEGRATION FOR EXCESS 
ATTENUATION. 

VALUE OF RANGE IN INTEGRAL FOR EXCESS ATTENUATION. 

CORRELATION LENGTH USED IN CALCULATION OF EXCESS 
ATTENUATION. 

FRACTION OF ENERGY BACKSCATTERED FROM GIVEN RANGE. 


INPUT BY 


Speed_sound SPEED OF SOUND IN AIR AT INPUT TEMPERATURE. 
Sumpow_back SUM OF THE BACKSCATTERED ENERGY 


1 
Temp 
Temp_con 


Temps 


Titles 
Tstar 


Var 
Wat_pres 


X 
X$ 


INPUT TEMPERATURE IN DEGREES KELVIN. 

INPUT TEMPERATURE IN DEGREES CELSIUS. 

USED TO INSURE VALUE OF TEMPERATURE PASSED TO SUBPROGRAM 
Att_temp WAS NOT CHANGED. 

VARIABLE STRING USED IN FUNCTION YES. 


STRING PASSED TO SUBPROGRAM PT FOR TITLE OF PLOT. 
INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 
IN SUBPROGRAM ATTENUATION. 

USED TO SELECT THE VARIABLE THE OPERATOR WISHED TO 
CHANGE BEFORE MAKING ANOTHER RUN. 

ATMOSPHERIC WATER PRESSURE IN MILLIBARS. 
OPERATOR. 

THIRD ORDER INDEX USED IN VARIOUS LOOPS. 
STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 


INPUT BY 


oe 


RESPONSE TO YES OR NO QUESTION. 
XlabelS LABEL ON X AXIS PASSED TO SUBPROGRAM Pt FOR PLOTTING. 


Xmax VALUE OF LARGEST VALUE OF X FOR VARIOUS PLOTS. 
Xmin VALUE OF SMALLEST VALUE OF X FOR VARIOUS PLOTS. 
Xrange VALUE OF RANGE OF X VALUES FOR VARIOUS PLOTS. 


Ylabel$ LABEL ON Y AXIS PASSED TO SUBPROGRAM Pt FOR PLOTTING. 


Y¥max LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

Ymin SMALLEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

Yrange RANGE OF Y VARIABLES FOR VARIOUS PLOTS. 

Ze INTERMEDIATE VALUE IN CALCULATION OF EXCESS ATTENUATION. 


| 
! 
{ 
! 
! Xvar(*#) VALUE OF ARRAY OF X VALUES TO BE PLOTTED IN SUBPROGRAM Pt. 
| 
{ 
! 
| 


| 

!DECLARE VARIABLES 

INTEGER I 

DIM Pow_back( 1500) ,Pow_ret( 1500) , Range( 15@@) ,Cte2( 1500) , Sigma’ 1500) 
DIM Cve2( 1500) ,Cne2( 1500) ,Exc_att( 1500) , Xvar( 1500) 

DIM TitleS{50) ,XlabelS{50] ,Ylabel${16] 

| 

PeeiTER eS 70S," HPGL” 

EINES TYPE 1 

I 

!INPUT ATMOSPHERIC DATA 

INPUT “ENTER TEMPERATURE IN DEGREES CELCIUS” , Temp 

INPUT “ENTER ATMOSPHERIC PRESSURE IN MILLIBARS" ,Atom_pres 
INPUT “ENTER WATER VAPOR PRESSURE IN MILLIBARS" ,Wat_pres 
| 

| 

{INPUT ECHOSOUNDER DATA 

INPUT “ENTER FREQUENCY OF ECHOSOUNDER IN HERTZ" ,Freq 
INPUT “ENTER ANTENNA DIAM IN METERS" ,Ant_diam 

INPUT “ENTER ECHOSOUNDER PULSE LENGTH IN MILLISECONDS=" ,Pulse 
INPUT “ENTER POWER TO TRANSMITTER IN WATTS" ,Pow_trans 


Et#.25 ! TRANSMIT EFFICIENCY 
Er=.25 TREGEIVER ERFICIENEY 
G=.40 'ANTENNA EFFECTIVE APERATURE FACTOR 


| 

[TSebEGrettE PROFILE FOR THE TEMPERATURE STRUGTURE 

FRONT “YOU HAVE A CHOICE OF TEMPERATURE STRUCTURE PROFILES” 
Paent “PROFILE TO BE USED.” 

PRINT “YOUR SBEECTIONS ARE” 


PRINT " 1 A TEMPERARURE STRUCTURE PROFILE BASED ON DATA AS ” 
PRINT “ PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 WHICH “ 
PRINT “ GIVES A HEIGHT TO: THE -1.16 PROFILE” 

FRINT? ©.“ 

PRI Nie” 2 THE TEMPERATURE STRUCTURE PROFILE ABOVE BUT" 
ERINT WITH A FACTOR OF TWO TO APPROXIMATE LOOKING" 
PRINT " UP A CONVECTIVE PLUME" 

PRINT =~ 


ow 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
| 

PRINT 
PRINT 
Mess_ 
WHILE 


END W 
OUTPU 
[ 

Again 


iy] 3 

ae 4 
up= | 

Mess_up=! 


A TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS" 
PRESENTED IN WALTERS/KUNDEL 198! PAGE 398 FOR” 
NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT " 
IS EXP(-.@@1*HEIGHT ABOVE 65 METERS AND HEIGHT" 
TO THE -1!.46 BELOW 65 METERS” 


A TEMPERATURE STRUETURE PROEICE TEASED ON AX 
HEIGHT TO THE -4/3 “ 


INPUT "ENTER THE DESIRED PROFILE (1 OR 2 OR 3 OR 4)" Profile 
IF Profile=4 THEN 
INPUT "ENTER THE HEIGHT OF THE INVERSION IN METERS” Inver 


END IF 


IF Profile=! THEN 


Mess_up=® 


ELSE 


IF Profile=2 THEN 
Mess_up=9 


EE oe 


IF Profile=3 THEN 


ESE 


Mess_up=@ 


IF Profile=4 THEN 
Mess_up=@ 

ELSE 
PRINT Profile,” WAS NOT ONE OF THE OPTIONS!!! !" 
Mess_up=! 

END IF 


END IF 


END IF 
END IF 
HILE 
T KBD3"K"; 


=| 


!INITIALIZE ARRAYS FOR SUCCESSIVE RUNS 


WHILE 
FOR 


Again=! 
J=! TO I 


Pow _back(J)=@ 
Pow_ret(J)=@ 


Range( J )=@ 
Cte2(J)=@ 
Sigma(J)=@ 
Cve2(J)=@ 
Cne2( J)=@ 
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Exc_att(J)=@ 

NEXT J 
FOR J=1 TO 15 

PRINT “ = 
NEXT J 
PRINT “DO YOU WISH TO SEE A PLOT OF ATTENUATION VERSUS HUMIDITY" 
PRINT "FOR VARIOUS FREQUENCIES AROUND THE FREQUENCY YOU INPUT?" 
LINPUT “IF YES ENTER Y , IF NO ENTER N"_ XS 
Att _freq=FNYes( X$) 
IF Att_freq=1 THEN 

OUTPUT KBD; "K"; 

CALL Att_freq( Atom_pres ,Freq, Temp) 
END IF 
OUTPUT KBD: "K"; 
I 
FOR J=1 TO 15 

PRINT " " 
NEXT J 
PRINT "DO YOU WISH TO SEE A PLOT OF ATTENUATION VERSUS HUMIDITY" 
PRINT "FOR VARIOUS TEMPERATURES AROUND THE TEMPERATURE YOU INPUT?" 
LINPUT "IF YES ENTER Y , IF NO ENTER N",XS 
Att_temp=FNYes( X$) 
IF Att_temp=1 THEN 

OUTPUT KBD;"K"; 

CALL Att_temp( Atom_pres ,Freq, Temp) 
END IF 
OUTPUT KBD;"K"; 
I 
{CONVERT TEMPERATURE TO KELVIN 
T=Tempt+273 
I 


ICALCULATE SPEED OF SOUND BASED ON INPUT TEMPERATURE‘ CELCIUS) 
Speed_sound=2@.05#(T)*.5 

{CALCULATE THE SPEED OF SOUND AT @ DEGREES CELCIUS 
C=20.05#273% .5 

i 

ICALCULATE THE ATTENUATION COEFFICIENT BASED ON ATMOSPHERIC DATA 
{EQUATION IS FROM NEFF 1975 

CALL Attenuation(Atom_pres ,Atten,Freq,  Temp,Wat_pres) 


I=] | INDEX FOR LOOP 
Pow_back( @)=@ fINITIALIZE VARIABLE FOR POWER BACKSCATERRED 
K=Z*PI*Freq/Speed_sound |WAVENUMBER 


Ant_area=PI#*( Ant_diam/2)*2Z |'ANTENNA AREA 

Interval =( Speed_sound*Pulse#!.E-3)/2 

Bn=2*Freq*( 1-1/(6/Speed_sound+1!))!BANDWIDTH FOR 3M/S VERTICAL VELOCITIES 
Noise=1.38E-23*Bn*(T)+Z2.E-14'IMINIMUM DETECTABLE SIGNAL AS 

| JOHNSON NOISE(NEGLIGABLE)+ ESTIMATED BACKGROUND 

I 


ICALCULATE RETURNED SIGNAL POWER UNTIL IT IS LESS THAN BACKGROUND. 


SS 


ZS a S_-  —_—_- ll 


Range( @)=@ 
Sumpow_ back =@ 


REPEAT 


Range( I)=Range(I-1)+2 
SELECT) Faeotire 


CASE 1 
I'THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
1AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398. 
IWAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
| TULAROSA BASIN, NEW MEXICO. 
1AN ADDITONAL FACTOR OF 1.8 WAS INCLUDED AS PER DIFFERENT 
{AVERAGING TIME. 
Cte2(1)=2.12*Range(I)*(-1.16) 
Chse Zz 
ITHIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
1AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398. 
IWAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
|! TULAROSA BASIN, NEW MEXICO. 
!AN ADDITONAL FACTOR OF {.8 WAS INCLUDED AS PER DIFFERENT 
|\AVERAGING TIME AND AN ADDITIONAL FACTOR OF TWO WAS 
|! INCLUDED TO APPROXIMATE LOOKING UP A CONVECTIVE PLUME 
Cte2(I)=2#*2.12*Range(I)*(-1.16) 
EASEao 
ITHIS TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS” 
IPRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 FOR" 
INIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT " 
11S EXP(-.@01*HEIGHT ABOVE 65 METERS AND HEIGHT" 
1G THe = 1246 BELOW Sa eneteRs: 
IF Range(I)<65 THEN 
Cte2(I)=75.5*Range(I)>*(-1.46) 
ELSE 
CteZ(I)=3.66E-Z2*EXP(-.@01*Range(I)) 
END IF 
CASE 4 
ICALCULATE THE TEMPERATURE STRUCTURE FACTOR BASED ON ASSUMPTION 
| THAT IT'S PROPORTIONAL TO HEIGHT TO THE -4/3 POWER 
| THAT THE SURFACE VERTICAL HEAT FLUX IS .@95 C.m./sec 
ITHE DEPENDENCE WITH HEIGHT ABOVE THE INVERSION I PULLED FROM 
Whe SS Siale 
LTEQUATION (-ROMENEFF Fala7S 
SELECT Range(I)/Inver 
EASE «ag 
C3=((.024)*(T)*¢(.667) >) 
Cte2(I )=C3*( Range(I))*(-1.33) 
Gfice 328 10m! 
Cte2Z(1)=Ctez( f-ip 
R1=Range( I) 
Ci=CteZ( 1) 
CASESIa TG 2123 
CteZ(1)=10° (({LGT( €3#Ranget 1) 7-UGT( C3*RI) 4G Cl) 
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RZ=Range( I) 
C2=Cte2(1) 
CASE GEESE 
Cte2Z( I >)»C34#( Ragged I)*(-1.33)-R2Z2°(-1.33))+4+E2 
ENDP SELECT 
END SELECT 
{CALCULATE VELOCITY STRUCTURE FACTOR. FORMULA FROM GAYNOR 77 
1!AND TETARSKI 
Epsilon=(.2866/Range(I))#( 1+.07#( Range(I))*(.6))*(1.5) 
Cve2(1I)=24#Epsilon*’(.667) 
'CALCULATE ACOUSTIC REFRACTIVE INDEX STRUCTURE FACTOR. 
!FORMULA FROM TETARSKI 1961 
Cne2(1)=(Cte2(1)/(2.98E+5) )4+( Cve2(1)/(C#C)) 
!CALCULATE THE FRACTION OF POWER BACKSCATTERED FROM INTERVAL. 
!'EQUATION FROM NEFF, 1975 
Sigma( I )=(.O0@39"(K*%( 1/3) )#Cte2(1))/(T)%2 
{CALCULATE THE EXCESS ATTENUATION. EQUATION FROM CLIFFORD, 19890 
1'THE EXCESS ATTENUATION IS Exc_att(I) 
!'THE MODULO STRUCTURE IS TO SKIP SOME OF THE INTEGRALS ONCE THE 
{RESOLUTION IS LESS IMPORTANT. 
IF I<S@ THEN 
Remainder=@ 
ELSE 
IF I<15@ THEN 
Remainder=Range(I) MODULO 19 
BSe 
Remainder=Range(I) MODULO 29 
END IF 
END IF 
IF Remainder=@ THEN 
Rho=@ 
L=@ 
H=@ 
Rgoe=Range( I) {CONSTANT IN INTEGRAL 
R=Q 
FOR J=@ TO 2¢I 
F=Cne2(INT( J/2+1)) 
F=Fe(1-R/(Roe)*(1.67)4+0(R/( Roe) )*(1.67)) 
IF J>@ THEN 
IF J<2*I THEN 
IF INT(J/2)=J3/2 THEN 
L=L+F 
F=0 
ja Lis |S 
H=H+F 
F=Q 
END Te 
END IF 
END IF 


©) 


Rho=RhotF 

R=R+1 
NEXT J 
Rho=Rhot4#L+2+*H 
Rho=(((Rho#.33)*#K*#K*1.46)°(-.6)) 
N=(Ant_diam/Rho)*2 
IF N<=1 THEN 


Ze=1/(14N) 
ERS 
Ze=1.5/(1+N) 
{STEP OF 1.S5==>SEE CLIFFORD 1989 
END IF 
Exc_att(I)=Ze+Ze 
ELSE 
Exc_att(I)=Exc_att(I-1) 
END IF 


ICALCULATE THE POWER BACKSCATTERED 
Pow _back( I )=(Pow_trans*Et-Sumpow_back )#*EXP(-Atten*Range( I )) 
Pow back(I)=Pow_back(I)*#Interval*Exc_att(1)*Sigma(I) 
Sumpow_back=Sumpow_back+P ow_back( I) 
!CALCULATE THE POWER BACKSCATTERED TO THE ANTENNA 
!'TWO WAY PATH ATTENUATION ACCOUNTED FOR ABOVE. 
Pow _ret(I)=Pow_back(I)*EXP(-~Atten*Range( 1) )*Ant_area*G*Er/Range(I1)*2 
PRINT “RANGE=",Range( I ) 
PRINT “POWER RETURNED=" ,Pow_ret(I) 
T=I+1 
UNTIL Pow _ret(1I-1)>@ AND Pow _ret(I-1)<=Noise 
OUTPUT KEDi ak 
! 
PRINT 


“INPUT CONDITIONS” 
i 
PRINT " “ 
PRINT USING (Ky EhRERA TURE " Lempe CEECIUS: 
PRINT USING "K";"ATMOSPHERIC PRESSURE= ",Atom_pres, mb" 
PRINT USING "K";"WATER VAPOR PRESSURE= ", Wat_pres,“mb" 
PRINT USING “K"; “PULSE LENGTH= "Pulse, "ms" 
PRINT USING "“K"; "TRANSMITTED FREQUENCY= ° Freq," Hz" 
PRINT USING "K"; "ANTENNA DIAMETER= ",Ant_diam," m." 
PRINT USING "K";"POWER TRANSMITTED= ",Pow_trans," WATTS" 
PRINT “TEMPERATURE STRUCTURE PROFILE USED " ,Profile 


IF Profile=4 THEN 
PRINT “INVERSION HEIGHT °, Inver 


END IF 
SELECT Peoiite 
CAseml 
PRINT " 1 A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
PRINT “ AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 “" 
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PRINT “ WHICH GLYES A HEIGHT TOsTRHE =i... 16 .RRGEIWE: 
Ese Z 
PRINT " A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
PRINT " AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 " 
PRINT " WHICH GIVES A HEIGHT TO THE -1.16 PROFILE” 
PRINT “ BUT WITH A FACTOR OF TWO TO APPROXIMATE" 
PRINT " LOOKING UP A CONVECTIVE PLUME" 
CASE 3 
PRINT " A TEMPERATURE STRUCTURE PROFILE BASED ON DATA" 
PRINT " AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398" 
PRINT " FOR NIGHT CONDITIONS. THE DEPENDENCE WITH" 
PRINT " HEIGHT 15 EXP(-.@01#HEIGHT ABOVE SS WETERS” 
PRINT " AND HEIGHT TO THE -1.46 BELOW 65 METERS” 
CASE 4 
PRINT " 4 A TEMPERATURE STRUCTURE PROFILE BASED ON A” 
PRINT " HEIGHT TO THE -4/3 " 
END SELECT 
{ 
PRINT USING "K";"THE MINIMUM DETECTABLE SIGNAL WAS SET AT” ,Noise,"WATTS” 
{ 
eeiINT ~ ° 
PRINT “OUTPUT CONDITIONS” 
{ 
PRINT USING “K";"RANGES=” ,Range(I-2)," m." 
! 
PRINT “HIT CONTINUE TO CONTINUE" 
PAUSE 
OUTPUT KBD; "K"; 
PRINT “THE FOLLOWING GRAPHS WILL NOW BE PLOTTED WITH RANGE VERSUS" 
PRINT “EXCESS ATTENUATION" 
PRINT "CTEZ (TEMPERATURE STRUCTURE PARAMETER)” 
men CVEZ AVECOCTTY STRUCTURE PARAMETER) * 
PRINT "“CNEZ (ACOUSTIC REFRACTIVE INDEX STRUCTURE PARAMETER)" 
PRINT “POWER RETURNED" 
PRINT “WHEN YOU ARE READY FOR THE FIRST GRAPH HIT CONTINUE” 
PAUSE 


OUTPUT KBD;"K"; 


GRAPHICS ON 
VIEWPORT 15,120,10,70 


| 
lEGor 


Xmin=@ 
Xmax= 1 
Ymin=@ 


RANGE VERSUS EXCESS ATTENUATION 


IF Range(I-Z2)<=1000 THEN 
Ymax=3 


EESE 


Ymax =4 


oe 


END IF 
Dx=.1 
Dy=1 
Xrange=Xmax- Xmin 
Yrange=Ymax-Ymin 
WINDOW Xmin, Xmax , Ymin, Ymax 
AXES Dx ,Dy,Xmin, Ymin,1,1 
CEI Pr aGRE 
| 
[LABEL PLOT 
CSI ZEV4 
LDIR @ 
LORG 5S 
MOVE .5,Ymax#l. 1 
LABEL “RANGE VERSUS EXCESS ATTENUATION" 
ILABEL HORIZONTAL AXES 
LDIR @ 
LORG S 
FOR J=@ TO Xmax STEP Dx 
MOVE J,-.1*#Dy 
LABELS 
NEXT J 
MOVE .S*Xrange,-.3#Dy 
CSIZE- 43-5 
LABEL “EXCESS ATTENUATION" 
ILABEL VERTICAL AXES 
LORG 8 
FOR J=@ TO Ymax STEP Dy 
CSIP ZE=45 
MOVE =23*Dx- 
UABEE. 10" 
CSIiZE 2 
MOVE -.1#Dx,J+.@5*Dy 
LABEL USING "K";J 
NEXT J 
LDIR PI/Z 
LORG 6 
MOVE -Dx,.5#Yrange 
CSIZEG4 765 
LABEL "RANGE ( METERS)" 
CLIP ON 
FOR Jz t=] 
PLOT Exc_att(J),LGT( Range( J) ) 
NEXT J 
PRINT “HIT CONTINUE FOR NEXT PLOT" 
i 
IPLOT RANGE VERSUS TEMPERATURE STRUCTURE PARAMETER 
PAUSE 
GCLEAR 


Ze 


| 
IF Range(I-2)<=1@@@ THEN 
Xmin=-4 
Ymax?3 
ELSE 
Xmin=-5 
Ymax=4 
END IF 
Xmax=@ 
Ymin=@ 
TitleS="RANGE vs TEMPERATURE STRUCTURE PARAMETER" 
XlabelS=""TEMPERATURE STRUCTURE PARAMETER" 
YlabelS="RANGE ( METERS)" 
CALL Log_log( Xmin, Xmax , Ymin, Ymax ,Cte2(*),I-1,Range(*),I-1,Title$,Xlabel$ 
11$) 
i 
| 
IPLOT RANGE VERSUS VELOCITY STRUCTURE PARAMETER 
PAUSE 
GCLEAR 
! 
Xmin=-2 
Xmax=@ 
Ymin=@ 
IF Range(I-2)<=1000 THEN 
Ymax=3 
BIESE 
Ymax=4 
ENG IF 
TitleS="“RANGE VERSUS VELOCITY STRUCTURE PARAMETER" 
XlabelS="VELOCITY STRUCTURE PARAMETER" 
CALL Log_log( Xmin, Xmax , Ymin, Ymax ,Cve2(#) ,I-1,Range(*#),I-1,TitleS,Xlabel$S 
219) 
| 
I 
IPLOT RANGE VERSUS ACOUSTIC REFRACTIVE INDEX STRUCTURE PARAMETER 
PAUSE 
GCLEAR 
Xmin=-7 
Xmax=-5 
Ymin=Q 
IF Range(I-2)<=1@0@ THEN 
Ymax=3 
PSE 
Ymax =4 
END IF 
TitleS="RANGE vs ACOUSTIC INDEX STRUCTURE PARAMETER" 
Xlabel$="ACOUSTIC REFRACTIVE INDEX STRUCTURE PARAMETER" 
CALL Log_log( Xmin, Xmax , Ymin, Ymax ,Cne2(*) ,I-1,Range(*) ,I-1,TitleS,XlabelS 
elS) 


(| 


! 
| 
{PLOT RANGE VERSUS POWER RETURNED TO ANTENNA 
PAUSE 
GCLEAR 
IF Range(I-2)<=1000 THEN 
Xmin=-15 
Xmax=-5 
Ymax=3 
EESE 
Xmin=-17 
Xmax=-4 
Ymax=4 
ENO IF 
Ymin=@ 
TitleS="RANGE VERSUS POWER RETURNED TO ANTENNA” 
Xlabel$="POWER RETURNEO TO ANTENNA” 
YlabelS="RANGE ( METERS)" 
CALL Log_log( Xmin, Xmax , Ymin, Ymax ,Pow_ret(*#),I-1,Range(*),I-1,TitleS,Xlab 
YlabelS) 
PAUSE 
GCLEAR 
| 
| 
LINPUT “WOULD YOU LIKE TO MAKE ANOTHER RUN (CENTER Y OR N)?" XS 
Again=FNYes( XS) 
Mess_up=1 
WHILE Mess _up=1 
SELECT Again 
GASEW 
LINPUT "WOULD YOU LIKE TO CHANGE A VARIABLE (ENTER Y OR N)?" XS 
New_va=FNYes( XS) 
Ano_change=1 
WHILE Ano_change=1 
SELECT New_va 


CASE } 

PRINT “ VARIABLE CURRENT VALUE" 

PRINT USING "K";"1 TEMPERATURE “, lemp, CEUSTMss 
PRINT USING "K"3;"2 ATOMOSPHERIC PRESSURE ",Atom_pres,"mb" 
PRINT USING "K";"3 WATER VAPOR PRESSURE ",Wat_pres,"mb" 
PRINT “ " 

PRINT USING “K";"4 FREQUENCY OF ECHOSOUNDER "Freq," Hz" 
PRINT USING "K";"5 ANTENNA DIAMETER ' Ant eianiee mo: 
PRINT USING “K"3 36. “PULSE EENGii Slate Bie” 
PRINT USING "“K";"7 POWER TRANSMITTED ",Pow_trans," WATTS" 
PRINT USING "K";"8 ATMOSPHERIC PROFILE ",Profile 

PRINT " “ 

PRINT "ENTER THE NUMBER OF THE VARIABLE YOU WISH TO" 

PRINT “CHANGE” 


CZ 


INPUT Var 


SELECT Var 
CASE } 
INPUT “TEMPERATURE=" , Temp 
CASE 2 
INPUT “ATMOSPHERIC PRESSURE IN mb="Atom_pres 
CASE 3 
INPUT “WATER VAPOR PRESSURE IN mb=" Wat_pres 
CASE 4 
INPUT “FREQUENCY IN Hz=" Freq 
CASE 5 
INPUT “ANTENNA DIAMETER IN m.=" ,Ant_diam 
CASE. 6 
INPUT “PULSE LENGTH IN ms=" Pulse 
CASE 7 
INPUT "POWER TRANSMITTED IN WATTS=" ,Pouw_trans 
CASE 8 
PRINT “ 1 =e Ct*Z PROFILE OF Z“(-1.15)" 
PRINT “ FROM WALTERS/KUNDEL 13981" 
PRINT " 2 ==> SAME AS ONE BUT WITH FACTOR" 
PRINT " OF TWO TO APPROXIMATE LOOKING" 
PRINT " UP A CONVECTIVE PLUME" 
PRINT " 3 ==> Ct*Z2 PROFILE OF EXP{(-.@@01#Z)" 
PRINT “ FROM WALTERS/KUNDEL 1981." 


PRINT ~ 4 == Gt? IROFIWe OF Z2°(-4/3)" 
INPUT “ENTER NUMBER OF DESIRED PROFILE" ,Profile 
IF Profile=4 THEN 
INPUT “HEIGHT OF INVERSION IN METERS=" Inver 
END IF 
CASE@EL SE 
PRINT Var,"IS NOT ONE OF THE OPTIONS" 
ENOP SELECT 
LINPUT “MADE ANOTHER CHANGECENTER Y OR N)?" XS 
New_va=FNYes( X$) 
Mess_up=2 
CASESZ 
Mess_up=2 
Ano_change=Z 
CASE, Etoe 
PRINT "YOUR CHOICES WHERE Y OR N !lllt" 
Ano_change= | 
LINPUT “CHANGE A VARIABLE (ENTER Y OR N)?", XS 
New_va=FNYes( X$) 
END SELECT 
END WHILE 
CASE 2 
Mess_up=2 
Gace ELSE 
PRINT “YOUR CHOICES WHERE Y OR N fletee" 
Mess_up=1 


{® 


LINPUT “MAKE ANOTHER RUN (ENTER Y OR N)?" | X& 
Again=FNYes< XS) 
END SELECT 
END WHILE 
END WHILE 
PRINT “THAT’S ALL, FOLKS" 
end 
l 
ICALCULTE THE ATTENUATION 
sub Attenuation(Atom_pres,Atten,Freq, Temp,Wat_pres) 
ITHIS SUBPROGRAM CALCULATES THE ATTENUATION OF SOUND 
1IN AIR BASED UPON EQUATIONS IN NEFF 1975 
i 
| INPUT ATMOSPHERIC PRESSURE IN MILLIBARS 
FREQUENCY OF SOUND WAVE IN HERTZ 
| TEMPERATURE IN DEGREES CELCIUS 
WATER-VAPOR PRESSURE IN MILLIBARS 
| 


OUTPUT ATTENUATION IN 1/METERS 


VARIABLES 
Atom_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 
Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 


| 

| 

! 

| 

| 

| 

SUBPROGRAM ATTENUATION. 

| = Att_max VARIABLE IN SUBPROGRAM ATTENUATION. IT IS THE 
‘ ATTENUATION AT THE FREQUENCY OF THE MAXIMUM 

| ATTENUATION FOR THE INPUT CONDITIONS. 

| E VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
| 

| 

| 

| 

| 

! 

| 

| 

| 


RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 

Fmax FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
ATTENUATION. 

Freq INPUT FREQUENCY OF ECHOSOUNDER. 

H VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 

Pstar VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 

Temp INPUT TEMPERATURE. ENGDEGREES CEESIMS. 

Tstar INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 


IN SUBPROGRAM ATTENUATION. 
Wat_pres ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 

! OPERATOR. 
H=1@0*Wat_pres/Atom_pres 
Tstar=(1.8*Tempt492)/519 
Pstar=Atom_pres/1014 
Fmax=( 1046600*H+444002H#H)*Pstar/Tstar*.8 
Att_max=.Q078*#Fmax*Tstar*(-2.5)*EXP(7.77*( 1-1/Tstar)) 
F=Freq/Fmax 
Atten=(Att_max/304.8)#((.18#F)*24+(2*FeF/( 14F*#F))°2)° 25 
Atten=(Attent+1.74E-10*Freq*Freq)/4. 35 
SUBEND 
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| 

def FNYes( XS) 

(THIS FUNCTION INTERPRETS THE OPERATORS RESPONSES TO YES NO QUESTIONS 
! INPUT XS 

l 

| OUTPUT FNYes 


!VARIABLES 
! Temps VARIABLE STRING USED IN FUNCTION YES. 
! X$ STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 


| RESPONSE TO YES OR NO QUESTION. 
DIM TempS[ 1] 
Temp$[1,1])=TRIMS( XS) 
SELECT Temps 
CASE ay . y~ 
RETURN 1 
CASE N° an: 
RETURN 2 
CASE * * 
RETURN 1 
CASE SElse 
RETURN -2 
BNOMoELECE 
FNEND 
sub Att_freq(Atom_pres,Freq, Temp) 
| 
| THIS SUBPROGRAM PLOTS THE ATTENUATION OF SOUND VERSUS WATER- 
[VAeUeerReESSURE FOR FIVE OIFFERENT FREQUENCIES AT 1/73 OCTIVE 
!'INTERVALS AROUND THE INPUT FREQUENCY. 
| 
INPUT ATOMOSPHERIC PRESSURE IN MILLIBARS 
FREQUENCY IN HERTZ 
MEMPERAMURE IN VCEECIUS 


OUTPUT PLOT OF ATTENUATION VERSUS WATER-VAPOR PRESSURE 


VARIABLES 
Atom_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 
Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 


| 

| 

| 

SUBPROGRAM ATTENUATION. 

'  Att_freq VARIABLE USED TO DETERMINE IF OPERATOR 

l WANTS TO PLOT ATTENUATION VERSUS WATER_VAPOR 

l PRESSURE FOR VARIOUS FREQUENCIES AROUND THE INPUT 
| FREQUENCY. IF SO THEN PLOT IS DONE IN SUBPROGRAM 
l 

l 

| 

l 


Att _freq. 
Es SATURATION VAPOR PRESSURE AT GIVEN TEMPERATURE. 
Freq INPUT FREQUENCY OF ECHOSOUNDER. 


Freq_con USED TO INSURE VALUE OF FREQUENCY PASSED TO SUBPROGRAM 
Att_freq WAS NOT CHANGED. 


ig 


1 J FIRST ORDER INDEX FOR ASSORTED LOOPS. 

le SECOND ORDER INDEX FOR ASSORTED LOOPS. 

1 Temp INPUT. TEMPERAIURE IN DEGREES CEESEUS: 

! Temp_con USED TO INSURE VALUE OF TEMPERATURE PASSED TO SUBPROGRAM 
| Wat_pres ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
i 
1 
\ 
i 
1 


OPERATOR. 

X THIRD ORDER INDEX USED IN VARIOUS LOOPS. 

X$ STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 
RESPONSE TO YES OR NO QUESTION. 

Ymax LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 


Freq_con=Freq 
Temp_con=Temp 
Ymax=@ 
Freq=DROUND( Freq#2%(2/3) ,3) 
FOR J=@ TO 4 STEP .1 
Wat_pres=J 
CALL Attenuation( Atom_pres ,Atten,Freq, Temp, Wat_pres) 
IF Atten>Ymax THEN 
Ymax=Atten 
END IF 
NEXT J 
Ymax=PROUND( Ymax+ .@@5, -Z ) 
GRAPHICS ON 
VIEWPORT 15,120,190, 7 
WINDOW @,15,@, Ymax 
AXES .5,2005 C2 c 
EEIP OFF 
(LASERS 
ESI 7Ee4 6 
LDIR @ 
MOVE 7.5, Ymax*#1.1 
LORG 5 
LABEL “ATTENUATION VERSUS WATER-VAPOR PRESSURE" 
SS50Ze 496 
'LABEL HORIZONTAL AXES 
LDIR @ 
LORG 6 
FOR J=@ TO 15 
MOVE J,@ 
LABEL J 
NEXT J 
MOVE 7.5,-Ymax*.1 
CSIZE 4,.6 
LABEL "WATER-VAPOR PRESSURE mb" 
[LABEL VERTICAL AXES 
LORG 8 
CSIZE 4,.6 
FOR J=@ TO Ymax STEP .@} 
MOVE = 252) 


(6 


BABEL 4) 
NEXT J 
ERBER PI/2Z 
LORG 6 
MOVE -2,Ymax#.S 
CSIZE 4,.6 
LABEL "ATTENUATION 1/m" 
IPRINT INPUT CONDITIONS 


LDIR @ 

LORG 1 

MOVE S.5, Ymax 
CSIZE 4..6 


LABEL USING "K";"ATMOSPHERIC PRESSURE=",Atom_pres," mb" 
MOVE 5.5, Ymax+#.9S 
GSr7e 4-26 
LABEL USING "K";"TEMPERATURE=",Temp,"” CELSIUS" 
J=@ 
FOR L=i TO S 
FOR X=@. TO 6.5+J STEP .@S 
Wat_pres=X 
CALL Attenuation(Atom_pres ,Atten,Freq, Temp, Wat_pres) 
IF INT(L/2)=L/2 THEN 
EY NESE 2 
al sys 
DEVE TYPE J 
END IF 
PLOT X,Atten 
NEXT X 
LDIR @ 
LORG Z 
GSILZE 3,.6 
MOVE X,Atten 
LINE TYPE | 
LABEL USING "K";Freq,"Hz” 
Freq=DROUND( Freq#2°(-1/3),3) 
J=J+1.5 
NEXT L 
| 
PRINT “HIT CONTINUE TO CONTINUE" 
PAUSE 
GCLEAR 
FOR J=@ TO 15 
PRINT " “ 
NEXT J 
PRINT "WOULD YOU LIKE TO SEE THE SAME PLOT QNLY" 
PRINT "WITH ATTENUATION VERSUS RELATIVE HUMIDITY?" 
LINPUT “IF YES ENTER Y, ELSE ENTER N” XS 
Att _freq=FNYes( X$) 
IF Att_freq=! THEN 
OUTPUT KBD3;"K"; 


= 


WINDOW @, 100,@, Ymax 
CLIP ON 
AXES 5.,.005,0.072e2 
Cele “ORE 
IEABEL FLOT 
CSIZE 4,.6 
LDIR @ 
MOVE 5@,Ymaxl. | 
LORG 5S 
LABEL “ATTENUATION VERSUS RELATIVE HUMIDITY" 
CSIZE 4.265 
!LABEL HORIZONTAL AXES 
LDIR @ 
LORG 6 
FOR J=@ TO 100 STEP 19 
MOVE J,@ 
LABEL J 
NEXT J 
MOVE 5@0,-Ymax#. | 
CSI ZE V4 
LABEL “RELATIVE HUMIDITY so 
'LABEL VERTICAL AXES 


LORG 8 

CSIZE 4,.6 

FOR J=@ TO Ymax STEP .O} 
MOVE 1,J 
LABELS 

NEXT J 

Ree by Z 

LORG 6 

MOVE = tS. %mMax* 25 

CSlZE 4G 


LABEL "ATTENUATION }/m" 
IPRINT INPUT CONDITIONS 
LDIR @ 
LORG 1 
MOVE 4@, Ymax 
CSIZE 4.6 
LABEL USING "K";"ATMOSPHERIC PRESSURE=" ,Atom_pres," mb" 
MOVE 40,Ymax«*.95 
CStZE a= 6 
LABEL USING "K";"TEMPERATURE=",Temp," CELSIUS" 
CLIP ON 
T=Temp+273 
Es=10°(9.4-2353/T) 
J=@ 
Freq=DROUND( Fregq_con*#2*(-2/3),3) 
FOR L=1 TO 5 
FOR X=@. TO 4@4+J STEP .5 
Wat_pres=Es*Atom_pres*X/((Atom_pres-Es+#(1+X/10@) )*1@@) 


VS 


CALL Attenuation(Atom_pres ,Atten,Freq, Temp, Wat_pres) 
IF INT(L/2)=L/2 THEN 
ine, | YPE -Z 
FIESE 
LINE TYPE 1 
END IF 
PLOT X,Atten 
NEXT X 
LDIR @ 
LORG 2 
Csive- 3s, . 6 
MOVE X,Atten 
LINE TYPE 1 
LABEL USING “K";Freq,”Hz” 
Freq=DROUND( Freq*#2*( 1/3) ,3) 
J=J+11 
NEXT L 
PRINT "HIT CONTINUE TO CONTINUE” 
PAUSE 
GCLEAR 
OUTPUT® KBD: "K’’: 
EESE 
GCLEAK 
OUTPUT KBD;"K"; 
END IF 
Temp=Temp_con 
Freq=Freq_con 
SUBEND 
! 
l 
sub Att_temp(Atom_pres ,Freq, Temp) 
| 
| 
| THIS SUBPROGRAM PLOTS THE ATTENUATION OF SOUND VERSUS WATER- 
!'VAPOR PRESSURE FOR FIVE DIFFERENT TEMPERATURES AT 1@ DEGREE 


! INTERVALS AROUND THE INPUT TEMPERATURE. 
| 


{ 
| INPUT ATOMOSPHERIC PRESSURE IN MILLIBARS 

! FREQUENCY OF ACOUSTIC ENERGY IN HERTZ 
| TEMPERATURE TNSGGESIUS 

| 


OUTPUT PLOTS OF ATTENUATION VERSUS WATER-VAPOR PRESSURE 


VARIABLES 


Atom_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 
Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 
SUBPROGRAM ATTENUATION. 


t 
I 
| 
| 
\ 
| 
eewtt temp VARIABLE USED TO DETERMINE IF OPERATOR WANTS 
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l 


Es 

Freq 

I 

J 

L 

Temp 
Temp_con 


Wat_pres 


X 
XS 


Ymax 


Ymax =Q 


TO PLOT ATTENUATION VERSUS WATER_VAPOR PRESSURE 

FOR VARIOUS TEMPERATURE AROUND THE INPUT TEMPERATURE. 
THE PLOT IS DONE IN THE SUBPROGRAM Att_temp. 

SATURATION VAPOR PRESSURE AT GIVEN TEMPERATURE. 

INPUT FREQUENCY OF ECHOSOUNDER. 

MAIN INDEX FOR LOOP TO CHECK RETURNED POWER UNTIL 
FIRST ORDER INDEX GOR ASSORT EEeLOGOrs: 

SECOND ORDER INDEX FOR ASSORTED LOOPS. 

INPUT TERRBERATURE IN DEGREES CEESTUS: 

USED TO INSURE VALUE OF TEMPERATURE PASSED TO SUBPROGRAM 
Att_temp WAS NOT CHANGED. 

ATMOSPHERIC WATER PRESSURE IN@MTERI BARS. 
OPERATOR. 

THIRD ORDER INDEX USED IN VARTOUS LOOPS. 
STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 
RESPONSE TO YES OR NO QUESTION. 
LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 


INFUT BY 


Temp_con=Temp 
Temp=Tempt+2@ 
FOR J=@ TO 4 STEP .1 
Wat_pres=J 
CALL Attenuation( Atom_pres ,Atten,Freq, Temp,Wat_pres) 
IF Atten>Ymax THEN 
Ymax=Atten 
END IF 
NEXT J 
Ymax=PROUND( Ymax+.@05,-2Z) 
GRAPHICS ON 
VIEWPORT 15,120,10,70 
WINDOW @,15,0, Ymax 
AXES .S 200s 0404.2 
CLIP. OFF 
(LABEL PLOT 
CSi Ze 42 oO 
LDIR @ 
MOVE 7.5, Ymax*1. 1 
LORG 5S 
LABEL “ATTENUATION VERSUS WATER-VAPOR PRESSURE" 
CSI 7B 4=G6 
!LABEL HORIZONTAL AXES 
LDIR @ 
LORG 6 
FOR J=@ TO 15 
MOVE J,@ 
EAB ard 
NEXT J 
MOVE 7.5,-Ymax*.1 


BO 


CSiVeE 4 ..6 

LABEL “WATER-VAPOR PRESSURE mb" 

!'LABEL VERTICAL AXES 

LORG 8 

GS17E 4-76 

FOR J=@ TO Ymax STEP .@!1 
MOVE .25.J 
LABEL J 

NEXT J 

LDIR PI/2 

LORG 6 

MOVE -2,Ymax+#.5 

CSIZE 4,.6 

LABEL “ATTENUATION We 

'PRINT INPUT CONDITIONS 


LDIR @ 
LORG 1 
MOVE 6, Ymax 
eS1ZE 4.6 


LABEL USING “K";"ATMOSPHERIC PRESSURE=" ,Atom_pres," mb" 
MOVE 6,Ymax*.95 
GS ZEe4 6 
LABEL USING "K";"“FREQUENCY=" Freq," HZ" 
CLIP ON 
J=@ 
FOR L=1 T0 S 
FOR X=@. TO 6+J STEP .@5 
Wat_pres=X 
IF INT(L/2)=L/2 THEN 
LINE TYPE 3 
EUSE 
LENE SRYEE. 1 
END IF 
CALL Attenuation(Atom_pres,Atten,Freq, Temp,Wat_pres) 
PLOT X,Atten 
NEXT X 
LDIR @ 
EGRG Z 
CSrZe 3. .6 
MOVE X,Atten 
LINE, TYPE 1 
LABEL USING "K";Temp,"C" 
Temp=Temp- 10 
J=J+1.6 
NEXE L 
PRINT "HIT CONTINUE TO CONTINUE" 
PAUSE 
GCLEAR 
OUTPUT KBD; "°K": 
FOR J=0 70 15 


8 | 


PRINT “ “ 
NEXT J 
| 
PRINT “WOULD YOU ARE LIKE TO SEE THE SAME PLOT ONLY" 
PRINT “WITH ATTENUATION VERSUS RELATIVE HUMIDITY" 
LINPUT "IF SO ENTER Y, IF NOT ENTER N",XS 
Att_temp=FNYes( XS) 
IF Att_temp=! THEN 
OUTPUT KBO;"K"; 
Temp=Temp_con 
WINDOW @,1900,9, Ymax 
AXES 10,.905,90,9,2,2 
CLIP OFF 
{LABEL PLOT 
CSIZE 4226 
LDIR @ 
MOVE 5@, Ymax+#l.1 
LORG S 
LABEL “ATTENUATION VERSUS RELATIVE HUMIDITY” 
1LABEL HORIZONTAL AXIS 
LDIR @ 
LORG 6 
FOR J=@ TO 10@ STEP 1 
MOVE J,@ 
LABEL J 
NEXT J 
MOVE 5@,-Ymax*.] 
CSIZE 456 
PABER ene enit ve HUM mT TY” 
[EABED VERTICAL AXES 


LORG 8 

CS 2204s 

FOR J=@ TO Ymax STEP .@!1 
MOVE Z,J 
LABEL J 

NEXT J 

EDIR  Pi/Z 

LORG 6 

MOVE -10,Ymax+#.5 

CS Ze eae 


LABEL “ATTENUATION 1/m" 

'IPRINT INPUT CONDITIONS 

LDIR @ 

LORG 6 

MOVE 5@, Ymax 

CSIZE 4.776 

LABEL USING "K";"ATMOSPHERIC PRESSURE="Atom_pres," mb" 
MOVE 5@,Ymax+#.95 

CSIZE 4,.6 

[LABEL USING “kK”; PREQUENGY= {Freq 2. 


S2 


CLIP ON 
J=@ 
Temp=Temp-Z@ 
FOR L=! TO S 
T=Temp+273 
Es=10*(9.4-2353/T) 
FOR X=@ TO 9@-J STEP .S 
Wat_pres=Es*Atom_pres#X/((Atom_pres-Es#( 1+X/1@@) )#100) 
IF INT(L/2)=L/2 THEN 
LINE TYPE 3 
ELSE 
ENE TYPE 1 
END IF 
CALL Attenuation( Atom_pres,Atten,Freq, Temp ,Wat_pres) 
PLOT X,Atten 
NEXT X 
LDIR @ 
LORG 2 
Corze 3.5 
MOVE X,Atten 
ING TYPE +t 
LABEL USING "“K"“;Temp," C" 
Temp=Temp+1@ 
J=J+15 
NEXT L 
PRINT "HIT CONTINUE TO CONTINUE” 
PAUSE 
GCLEAR 
OUTPUT KBD3"K"; 
Eieoe 
OUTPUT KBD; "K"; 
END IF 
Temp=Temp_con 
SUBEND 
! 
| 
sub Log_log( Xmin, Xmax , Ymin, Ymax , Xvar(+#),J,Range(*#),L,Title$,Xlabel$,Ylabel 


THIS SUBROUTINE MAKES A LOG-LOG PLOT OF DATA PASSED FROM THE 
MAIN PROGRAM. 


| 
| 
| 
| 
! INPUT MINIMUM VALUES OF X AND Y FOR PLOT 
| MAXIMUM VALUES OF X AND Y FOR PLOT 
| X AND Y VALUES TO BE PLOTTED 
| TITRE SOE eon 

| LABELS FOR X AND Y AXIS 

OUTPUT LOG IESG PLOT 

| 


OS 


!VARIABLES 


Ox 

Dy 

nj 

ia 
Titles 
Xlabel$ 
Xmax 
Xmin 
Xrange 
Xvar(*) 
Ylabel$ 
Ymax 
Ymin 
Yrange 


Dy=1 
Xrange=ABS( Xmax- Xmin) 
Yrange=ABS( Ymax-Ymin) 
WINDOW Xmin, Xmax , Ymin, Ymax 
AXES Dx ,Dy,Xmin,Ymin,1!,1 
CLIP Gak 


'LABEL PLOT 
CSIZEe4 ac 


LDIR @ 
LORG 5 
MOVE Xmint+.5*Xrange,1.1#*Yrange 


STEP SIZE FOR X AXIS FOR VaRTOUSSEeoic. 

STEP SIZE FOR Y AXIS FOR VARI CUSS EEGs. 

FIRST ORDER INDEX FOR ASSORTED LOOPS. 

SECOND ORDER INDEX FOR ASSORTED LOOPS. 

STRING PASSED TO SUBPROGRAM PT FOR TITLE OF PLOT. 
LABEL ON X AXIS PASSED TO SUBPROGRAM Pt FOR PLOTTING. 
VALUE OF LARGEST VALUE OF X FOR VARIOUS PLOTS. 

VALUE OF SMALLEST VALUE OF X FOR VARIOUS PLOTS. 

VALUE OF RANGE OF X VALUES FOR VARIOUS PLOTS. 

VALUE OF ARRAY OF X VALUES TO BE PLOTTED IN SUBPROGRAM Pt. 
LABEL ON Y AXIS PASSED TO SUBPROGRAM Pt FOR PLOTTING. 
LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 
SMALLEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

RANGE OF Y VARIABLES FOR VARIOUS PLOTS. 


TitleS=TRIMS( Titles) 

LABEL Titles 

!'LABEL HORIZONTAL AXES 

LDIR @ 

LORG 5S 

FOR J=Xmin TO Xmax STEP Ox 
CSiZEs 4.0 
MOVE J-.@13*xXrange,-.9067#Yrange 
LABEL “10° 
MOVE J+.013*Xrange,-.033*Yrange 
CSiZE sz 
ABE 

NEXT J 

MOVE Xmint+.5*Xrange,-.12*Yrange 

LORG 5S 

CSONZE V4 as 

Xlabel $=TRIMS( Xlabel®&) 

LABEL Xlabel$ 

!LABEL VERTICAL AXES 

LORG 8 

FOR J=@ TO Ymax STEP Dy 


8 4 


CSrzZe 4,.6 
MOVE Xmin-.025*Xrange, J 
PABE “10” 
CSIZE 2 
MOVE Xmin- .Q@0@25* Xrange ,J+.@3*Yrange 
EABEL J 
NEXT J 
LDIR PI/2 
LORG 6 
MOVE Xmin-.1#xXrange,.S*#Yrange 
CSlLZE 4,.6 
Ylabel $=TRIMS( YlabelS) 
LABEL Ylabel$ 
CLIP ON 
FOR J=1 TOL 
PLOT LGT( Xvar(J)),LGT( Range( J) ) 
NEXT J 
SUBEND 
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leaeeeeaePURPOQSEseeuzae 


| THIS PROGRAM WILL TAKE PARAMETERS OF AN ECHOSOUNDER AND 


fESTIMATE THE 


!ARE REQUIRED: 


RANGE AS A FUNCTION OF FREQUENCY. THE FOLLOWING INPUTS 


ATMOSPHERIC DATA 


7) SELEGY 


Again 


Ano_change 


Ant_area 


Ant_diam 
Atom_pres 
Atten 


Bn 
C 
eI 


cs 


! 
! 
| 
! 
! 
! 
| 
! 
| 
{ 
| 
l 
! 
! 
| 
! 
| 
| 
| 
| 
| 
! 
| 
| 
| 
l 
| 
| 
| 
{ 
! 
l 
! 
| CneZ2¢(#) 


1 )ATOMOSPHERIC PRESSURE in millibars 


ONE OF FOUR ATMOSPHERIC PROFILES FOR THE TEMPERATURE 


STRUCTURESPROETEE- 
a)FOR ONE PROFILE THE HEIGHT OF THE INVERSION LAYER IS ALSO INPUT 
3) TEMPERATURE IN DEGREES CEESTUS 
4)WATER VAPOR PRESSURE IN millibars 


ECHOSOUNDER DATA 

S)ANTENNA DIAMETER IN METERS 

6)POWER TRANSMITTED BY ECHOSOUNDER IN WATTS 
7)PULSE LENGTH OF THE TRANSMITTED ACOUSTIC ENERGY 


THE PROGRAM OUTPUTS THE FOLLOWING GRAPHS TO AN EXTERNAL PLOTTER: 


1 )RANGE(m.) VERSUS FREQUENCY OF ECHOSOUNDER 
2 )RANGE(m.) VERSUS EXCESS ATTENUATION FOR VARIOUS FREQUENCUIES 


*#e#44#VARITABLES# 424% 


TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 
WISHES TO MAKE ANOUTHER RUN. 

TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 
WISHES TO CHANGE THE VALUE OF INPUT DATA BEFORE MAKING 
ANOUTHER RUN. 

ECHOSOUNDER ANTENNA AREA CALCULATED FROM INPUT OF 
ANTENNA DIAMETER. 

INPUT OF ECHOSOUNDER ANTENNA DIAMETER IN METERS. 

INPUT OF ATMOSPHERIC PRESSURE IN mb. 

ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 
SUBPROGRAM ATTENUATION. 

BANDWIDTH OF RECEIVER NEEDED TO RECEIVE SIGNALS 
DOPPLER SHIFFED BY THREE METER PER SECOND 
VERTICAL VEEGETIIEs- 

VELOCITY OF SOUND CALCULATION AT ZERO CELSIUS. 

VARIABLE USED TO SCALE INVERSION HEIGHT IN 
TEMPERATURE STRUCTURE PROFILE THREE. 

VARIABLE USED IN CALCULATING THE FOURTH CteZ PROFILE. 

ARRAY OF VALUES OF THE ACOUSTIC REFRACTIVE INDEX 
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PARAMETER. CALCULATED BASED ON SELECTION OF 
PROFILE FOR TEMPERATURE STRUCTURE PARAMETER 
AND ASSUMED PROFILE FOR THE VELOCITY STRUCTURE 


PARAMETER. 

Ctez VALUE OF THE TEMPERATURE STRUCTURE PARAMETER. 
THE OPERATOR SELECTS A PROFILE FROM SEVERAL PROVIDED. 

CveZz VALUE OF THE VELOCITY STRUCTURE PARAMETER. 
VALUES BASED ON CALCULATION USING ASSUMED DISSAPATION 
RATE. 

Dx STEP SIZE FOR X AXIS FOR VARIOUS PLOTS. 

Dy STEP SIZE FOR Yonnts FORP VARIOUS PEOTS. 

Epsilon DISSAPATION RATE USED IN CALCULATION OF Cvez2. 

3 TRANSMISSION EFFICIENCY OF ECHOSOUNDER. 

EI EFFICIENCY OF ECHOSOUNDER WHEN RECEIVING ACOUSTIC 
BACKSCATTER. 

Exc_att(*,#) EXCESS “ATTENUATION” AT GIVEN RANGE. 

F VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 

Fmax FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
ATTENUATION. 

eed INPUT FREQUENCY OF ECHOSOUNDER. 

. G6 ANTENNA EFFECTIVE APERATURE FACTOR. 

H VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 

I MAIN INDEX FOR LOOP TO CHECK RETURNED POWER UNTIL 
POWER WAS LESS THAN MINIMUM DETECTABLE SIGNAL. 

IND INDEX USED FOR LOOP FOR DIFFERENT FREQUENCIES 

Interval MAXIMUM EFFECTIVE-SCATTERING VOLUME THICKNESS. 

Inver HEIGHT OF INVERSION LAYER. 

J FIRST ORDER INDEX FOR ASSORTED LOOPS. 

K WAVENUMBER 

L SECOND ORDER INDEX FOR ASSORTED LOOPS. 

Mess_up TEST VALUE TO CHECK RESPONSES WHEN OPERATOR IS INPUTTING 
RESPONSES. 

N TEST VALUE IN CALCULATION OF COHERENCE LENTH IN CALCULATION 
OF EXCESS ATTENUATION. 

Neu_va TEST VALUE IN CHECKING IF OPERATOR WISHES TO CHANGE 
A VARIABLE BEFORE A NEW RUN. 

Noise ASSUMED MINIMUM DETECTABLE SIGNAL. 


Pow_back(*) POWER BACKSCATTERED FROM GIVEN RANGE. 

Pow_ret(*) POWER BACKSCATTERED TO ECHOSQUNDER FROM GIVEN RANGE. 
Pow_trans INPUT OF POWER SUPPLIED TO TRANSDUCER OF ECHOSOUNDER. 
Profile OPERATOR INPUT OF Cte2 PROFILE FROM AVAILABLE PROFILES. 


Pstar VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 

Pulse TRANSMITTED PULSE LENGTH IN MILLISECONDS. INPUT BY 
OPERATOR. 

R RANGE FOR INPUT TO INTEGRAL OF EXCESS ATTENUATION. 

R1 VARIABLE USED IN SCALING INVERSION HEIGHT FOR PROFILE THREE. 


Range( * } ARRAY OF RANGE VALUES. 
Remainder REMAINDER OF MODULO FUNCTION USED TO DECREES THE 


s) 


NUMBER OF PASSES THROUGH THE INTEGRATION FOR EXCESS 
ATTENUATION. 


Rge VALUE OF RANGE IN INTEGRAL FOR EXCESS ATTENUATION. 

Rho CORRELATION LENGTH USED IN CALCULATION OF EXCESS 
ATTENUATION. 

Sigma FRACTION OF ENERGY BACKSCATTERED FROM GIVEN RANGE. 


Speed_sound SPEED OF SOUND IN AIR AT INPUT TEMPERATURE. 
Sumpow_back SUM OF BACKSCATTERED ENERGY 


| 

! 

| 

| 

| 

! 

Le INPUT TEMPERATURE IN DEGREES KELVIN. 

! Temp INPUT TEMPERATURE IN DEGREES CELSIUS. 

! Temps VARIABLE STRING USED IN FUNCTION YES. 

{ 

| Titles STRING PASSED TO SUBPROGRAM Semi_log FOR TITLE OF PLOT. 
! star INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 
IN SUBPROGRAM ATTENUATION. 

! Var USED TO SELECT THE VARIABLE THE OPERATOR WISHED TO 

! CHANGE BEFORE MAKING ANOTHER RUN. 

| Wateeges ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPU!? BY 

! OPERATOR. 

1 Xx THIRD ORDER INDEX USED IN VARIOUS LOOPS. 

1 X% STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 

! RESPONSE TO YES OR NO QUESTION. 

! Xlabel$ LABEL ON X AXIS PASSED TO SUBPROGRAM Semi_log FOR PLOTTING. 
{ max VALUE OF LARGEST VALUE OF X FOR VARIOUS PLOTS. 

1 Xmin VALUE OF SMALLEST VALUE OF X FOR VARIOUS PLOTS. 

| range VALUE OF RANGE OF X VALUES FOR VARIOUS PLOTS. 

| Xvar¢(*) VALUE OF ARRAY OF X VALUES TO BE PLOTTED IN SUBPROGRAM Semi_ 
! Ylabel& LABEL ON Y AXIS PASSED TO SUBPROGRAM Semi_log FOR PLOTTING. 
! Ymax LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

! Ymin SMALLEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

| Yrange RANGE OF Y VARIABLES FOR VARIOUS PLOTS. 

! Ze INTERMEDIATE VALUE IN CALCULATION OF EXCESS ATTENUATION. 
! 


!DECLARE VARIABLES 
INTEGER I 
DIM Pow_back(15@@) ,Pow_ret( 15@@) ,Range( 1500) 

DIM Cne2( 1500) ,Exc_att( 1500, 30) , Xvar( 1500) ,Ran( 3@) 
DIM TitleS$(5Q0] ,Xlabel$(5Q] , Ylabel3%[16] 


PLOTTER 1S 7Os=s RCE. 
LINE iy Gees 


!'INPUT ATMOSPHERIC DATA 

INPUT “ENTER TEMPERATURE IN DEGREES CELCIUS" , Temp 

INPUT “ENTER ATMOSPHERIC PRESSURE IN MILLIBARS" ,Atom_pres 
INPUT "ENTER WATER VAPOR PRESSURE IN MILLIBARS" , Wat_pres 
| 

' 


sie 


! INPUT ECHOSOUNDER DATA 


INPUT 
INPUT 


“ENTER ANTENNA DIAM IN METERS" ,Ant_diam 
"ENTER ECHOSOUNDER PULSE LENGTH IN MILLISECONDS=" Pulse 


INPUT "ENTER POWER TO TRANSMITTER IN WATTS" ,Pow_trans 
By=.75 !ITRANSMIT EFFICIENCY 

Ers.25 !RECEIVER EFFICIENCY 

G=.40 'ANTENNA EFFECTIVE APERATURE FACTOR 


ISELECT THE PROFILE FOR THE TEMPERATURE STRUCTURE 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
| 

PRINT 
PRINT 
Mess_ 
WHILE 


"YOU HAVE A CHOICE OF TEMPERATURE STRUCTURE PROFILES" 

orPrRoriLew Wo "Se USED.” 

“YOUR SELECTIONS ARE" 

1 A TEMPERARURE STRUCTURE PROFILE BASED ON DATA AS " 
e PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 WHICH " 
GIVES A HEIGHT TO THE -1.16 PROFILE” 

‘ L THE TEMPERATURE STRUCTURE PROFILE ABOVE BUT" 

re WITH A FACTOR OF TWO TO APPROXIMATE LOOKING" 

7 UP A THERMAL PLOOM" 

pe 3 A TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS" 
i PRESENTEO IN WALTERS/KUNDEL 1981 PAGE 398 FOR" 

: NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT " 

2 IS EXP(-.@@1*HEIGHT ABOVE 65 METERS AND HEIGHT" 

g TO) jHGe—1.46 s8EHL0U 65 METERS” 

oe 4 A TEMPERATURE STRUCTURE PROFILE BASED ON A" 
HEIGHT 10 THE -4/3 ° 
up=1 

Mess_up=] 
INPUT "ENTER THE DESIRED PROFILE (1 OR 2 OR 3 OR 4)" Profile 


IF Profile=4 THEN 
INPUT “ENTER THE HEIGHT OF THE INVERSION IN METERS" , Inver 
ENO IF 
IF Profile=1 THEN 
Mess_up=@ 
ELSE 
IF Profile=Z THEN 
Mess_up=@ 
EES 
IF Profile=3 THEN 
Mess _up=@ 
ELSE 
IF Profile=4 THEN 
Mess_up=@ 
ELSE 
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PRINT Profile,” WAS NOT ONE OF THE GPiaoNSI ITT” 
Mess_up=1 
END IF 
END IF 
ENO IF 
END IF 
END WHILE 
OUTPUT KBD;"K"; 
j 
| 
Again=1 
WHILE Again=1 
FOR J=1 TO I 
Pow_back(J)=@ 
Pow _ret(J)=@ 
Range( J)=@ 
Cne2(J)=@ 
FOR L=f TO Ind 
Exe att( J), ,Ld=0 
NEXT L 
NEXT J 
| 
ICONVERT TEMPERATURE TO KELVIN 
T=Tempt+273 
[CALCULATE SPEED OF SOUND BASED ONS oi TEMPERATURE, CELCIUS) 
Speed_sound=2@. @5#(T)*.5 
|CALCULATE THE SPEED OF SOUND Al @sSEGREES CELeruUsS 
C=20-05* 7273-45 
| 


I=1 !INDEX FOR LOOP 
Ind=1i 
Pow_back(@)=@ lINITIALIZE VARIABLE FOR POWER BACKSCATERRED 


Ant_area=PI*(Ant_diam/2Z)*2 !'ANTENNA AREA 
Interval=( Speed_sound*Pulse*#l.E-3)/2 
I 
lLINITIALIZE PLOTTING PARAMETERS SO CAN COMPARE DURING RUN TO SCALE 
IPLOTS 
Ymin=100@ 
Ymax=@ 
ICALCULATE RETURNED SIGNAL POWER UNTIL IT IS LESS THAN BACKGROUND 
1FOR EACH FREQUENCY. 
Range( @)=@ 
FOR Freq=100 TO 3000 STEP 140 
Bn=2*Freq#( 1-1/(6/Speed_sound+!))! BANDWIDTH FOR 3M/S VERTICAL VELOCITY 
Noise=1.38E-23*Bn*(T)+2.E-14!MINIMUM DETECTABLE SIGNAL AS 
| JOHNSON NOISE( NEGLIGABLE)+ ESTIMATED BACKGROUND 
ICALCULATE THE ATTENUATION COEFFICIENT BASED ON ATMOSPHERIC DATA 
CALL Attenuation( Atom_pres ,Atten,Freaq, Temp,Wat_pres) 
K=2*PI*Freg/Speed_sound |! WAVENUMBER 


S© 


I=] 
Sumpouw_back =@ 
REPEAT 
Range( I )=Range(I-1)+2 
SEEBCT Profile 
CASE 1 
I'THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
1AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398. 
IWAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
1TULAROSA BASIN, NEW MEXICO. 
1'AN ADDITONAL FACTOR OF 1.8 WAS INCLUDED AS PER DOIFFERENT 
lAVERAGING TIME. 
Cte2=2.12*#Range(I)*(-1.16) 
CASE Z 
(THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
1AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398. 
IWAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
1 TULAROSA BASIN, NEW MEXICO. 
1AN ADDITONAL FACTOR OF 1.8 WAS INCLUDED AS PER DIFFERENT 
'AVERAGING TIME AND AN ADDITIONAL FACTOR OF TWO WAS 
1 INCLUDED TO APPROXIMATE LOOKING UP A THERMAL PLOOM 
Cte2=2%*2.12*Range(I)*(-1.16) 
CASE 3S 
ITHIS TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS" 
IPRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 FOR" 
INIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT “ 
11S EXP(-.@@1*HEIGHT ABOVE 65 METERS ANDO HEIGHT" 
1'TO THE -1.46 BELOW 6S METERS” 
IF Range(I)<65 THEN 
CteZ=75.5*Range(I)*(-1.46) 
EIESE 
CteZ=3.66E-2*EXP(-.901*#*Range(I)) 
ENO IF 
CASE 4 
ICALCULATE THE TEMPERATURE STRUCTURE FACTOR BASED ON ASSUMPTION 
| THAT IT’S PROPORTIONAL TO HEIGHT TO THE -4/3 POWER 
| THAT THE SURFACE VERTICAL HEAT FLUX IS .095 C.m./sec 
I THE DEPENDENCE WITH HEIGHT ABOVE THE INVERSION I PULLEO FROM 
IMY LEFT EAR 
!EQUATION FROM NEFF, 1975 
SELECT Range(I)/Inver 
GASES. 3 
C3=(( .0@24)#(T)*(.667)) 
Cte2=C3*#( Range(I))*(-1.33) 
CASEe.3, 10 | 
CteZ2=Ctez 
C1=CteZ 
R1l=Range( I) 
CASE 1 TO 1.3 
Cte2=10*((LGT(C3*#Range(I))-LGT(C3*#R1I))4+LGT(CI)) 


S| 


R2=Range( I ) 
C2=Cte2 
GASE ESE 
CteZ=C3*( Range(1I)*(-1.33)-RZ7(-1.33))+4€Z 
END SELECT 

ENG SELECT 
{CALCULATE VELOCITY STRUCTURE FACTOR. FORMULA FROM GAYNOR 77 
Epsilon=(.2866/Range(I))*(1+.07#( Range(I))*(.6))*(1.5) 
Cve2=2*Epsilon*( .667) 
!CALCULATE ACOUSTIC REFRACTIVE INDEX STRUCTURE FACTOR. 
'FORMULA FROM TETARSKI 
Cne2(1I)=(Cte2/(2.98E15) )+( Cve2/(C#C) ) 
'CALCULATE THE FRACTION OF POWER BACKSCATTERED FROM INTERVAL. 
'EQUATION FROM NEFF, 1975 
Sigma=( .0039#(K°%( 1/3) )#Cte2)/(T)*2 
i 


ICALCULATE THE EXCESS ATTENUATION. EQUATION FROM CLIFFORD, 1980 
'THE EXCESS ATTENUATION IS Exc_att(I,L,X) 
IF I<5@ THEN 
Remainder =@ 
SAS S 
IF 19150 THEN 
Remainder=Range(I) MODULO 190 
ELSE 
Remainder=Range(I) MODULO 2@ 
END IF 
END IF 
IF Remainder=@ THEN 
Rho=@ 
L=@ 
H=@ 
Roe=Range( I) 1'CONSTANT IN INTEGRAL 
R=@ 
FOR J=@ TO 2#I 
F=Cne2( INT(J/2+1)) 
F=F#e( 1-R/(Roe)*(1.67)4¢(R/(Rge) )*(1.67)) 
IF J>@ THEN 
IF J<2#I THEN 
IF INT( J/2)=J3/2 THEN 
L=L+F 
F=Q 
ELSE 
H=H+F 
F=Q 
END IF 
END IF 
END IF 
Rho=RhotF 
R=R+ 1 
NEXT J 


Se 


Rho=Rho+4#L+Z4#H 

Rho=((( Rhot.33)*#K#K#1.46)°(-.6)) 
N=(Ant_diam/Rho)*2 

IF N<=1 THEN 


Ze=1/( 14N) 
ESE 
Ze=1.5/(14N) 
'STEP OF 1.5==>SEE CLIFFORD 198@ 
END IF 
Exc_att(I ,Ind)=Ze#Ze 
ELSE 
Exc_att(I ,Ind)=Exc_att(I-1,Ind) 
END UIF 


ICALCULATE THE POWER BACKSCATTERED 
Pow_back(I)=(Pow_trans*#Et-Sumpow_back )#*EXP(-Atten*Range(!I )) 
Pow_back(I)=Pow_back(I)#Interval#*Exc_att(I,Ind)*Sigma 
Sumpow_back=Sumpow_back+Pow_back( TI) 

[CALCULATE THE POWER BACKSCATTERED TO THE ANTENNA 
Pow_ret(I)=Pow_back(I)#EXP( -Atten*Range(I))#Ant_areatG#tEr/Range(I)*2 
PRINT "“RANGE=" ,Range( I ) 

PRINT "POWER RETURNED=" ,Pow_ret(I) 

P=l+ 


UNTIL Pow _ret(I-1)>@ AND Pow_ret(I-1)<=Noise 
Ran( Ind)=Range(I-2Z) 
IF Range(I-2)>Ymax THEN 


Ymax=Range(I-2Z) 


ENS IF 
IF Range(I-2)<Ymin THEN 


Ymin=Range(I-2) 


END IF 
Ind=Ind+1 
Again=1 

FOR J=!1 TO I 


Pow_back( J)=@ 
Pow_ret( J)=@ 
Range( J)=@ 
Cne2( J)=@ 


NEXT J 

PRINT "FREQUENCY=" ,Freq 
NEXT Freq 

! 


OUTPUT KBD;"K"; 
PRINT “INPUT CONDITIONS" 


PRINT " “ 


PRINT USING "“K";" TEMPERATURE= " Temp, “CELCIUS" 
PRINT USING "K";"ATMOSPHERIC PRESSURE= " Atom_pres, "mb" 
PRINT USING "K";"WATER VAPOR PRESSURE= " Wat_pres, "mb" 


i), 


PRINT USING "K";"PULSE LENGTH= 
PRINT USING "K"; "ANTENNA DIAMETER= " Ant_diam,"“ m. 
PRINT USING "K";"POWER TRANSMI TTED= " Pow_trans," WATTS" 
PRINT "TEMPERATURE STRUCTURE PROFILE USED ",Profile 
IF Profile=4 THEN 

PRINT “INVERSION HEIGHT “, Inver 


" Pulse,"“ms 


END IF 
SELECT Profile 
CASE! 
PRINIs= A TEMPERARURE STRUCTURE PROFILE BASED ON DATA” 
PRINT " AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 " 
PRINT " WHICH GIVES A HEIGHT IG Tee +17 leeeROFILe: 
CASE ae 
PRINT “ A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
PRINT “ AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 “ 
PRINT: WHICH GIVES A HEIGHT TOmiHEe- 1. |GePROBILE” 
PRINT “ BUT WITH A FACTOR OF TWO TO APPROXIMATE" 
PRINT " LOOKING UP A THERMAL PLOOM" 
CASE = 
PRINT " A TEMPERATURE STRUCTURE PROFILE BASED ON DATA” 
PRINT “ AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 338" 
PRINT “ FOR NIGHT CONDITIONS. THE DEPENDENCE WITH" 
PRINT HEIGHT IS EXP(-.@@1*HEIGHT ABOVE 6&5 METERS" 
PRINT ° AND HEIGHT TO THE -1!.46 BELOW 6&5 METERS | 
CASE 4 
PRINT A TEMPERATURE STRUCTURE PROFILE BASED ON A 
PRINT " HEIGHT TO THE -4/3 " 
ENG SEEEET 


PRINT USING "K"3;"THE MINIMUM DETECTABLE SIGNAL WAS SET AT" ,Noise,"WATTS" 
PRINT em 

PRINT “HIT CONTINUE TO CONTINUE” 
PAUSE 

OUTPUT KBD;"K"; 

! 

GRAPHICS ON 

VIEWPORT 15,120,190, 70 

! 

! 

IPLOT RANGE VERSUS FREQUENCY 
Xmin=O 

Xmax = 300@ 

Ymin=PROUND( Ymin-5@, 2) 
Ymax=PROUND( Ymax+5@, 2) 

Dx=50@ 

Dy=100 

Xrange=Xmax-~-Xmin 
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Yrange=Ymax-Ymin 

WINDOW Xmin, Xmax , Ymin, Ymax 

AXES Dx,Dy,Xmin,Ymin,Z,] 

ELIP OFF 

! 

!LABEL PLOT 

CSI7E.4,.& 

LDIR @ 

LORG 5 

MOVE .S5*#Xrange, Ymax*#1.@05 

LABEL “RANGE VERSUS FREQUENCY" 

ILABEL HORIZONTAL AXES 

LDIR @ 

LORG 5 

FOR J=@ TO Xmax STEP Dx 
SSI 7E 425 
MOVE J, Ymin-.@5* Yrange 
CABEL J 

NEXT J 

MOVE .S+*Xrange, Ymin-.1+#Yrange 

Come 4.215 

LABEL “FREQUENCY” 

1'LABEL VERTICAL AXES 

LORG 8 

FOR J=Ymin TO Ymax STEP Dy 
Goetze 4,.5 
MOVE Xmin- .0@025* Xrange, J 
EABEE USING mak 3a 

NEXT J 

LDIR PI/Z 

LORG 6 

MOVE Xmin-.15*Xrange, .5* YrangetYmin 

ESIZE 4..6 

LABEL “RANGE ( METERS)” 

CLIP ON 

Ind=1 

FOR Freq=100 TO 3000 STEP 120@ 
PLOT Freq,Ran( Ind) 
Ind=Ind+1 

NEXT Freq 

I 


PRINT “HIT CONTINUE TO CONTINUE" 
PAUSE 

GCLEAR 

OUTPUT KBD; "K"; 

l 

l 

IPLOT RANGE VERSUS EXCESS ATTENUATION 
Xmin=@ 
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Ymax=4 
Xmax= | 
Ymin=@ 
Range( @)=@ 
TitleS="RANGE vs EXCESS ATTENUATION" 
XlabelS="EXCESS ATTENUATION" 
YlabelS="RANGE ( METERS)" 
Dx=.1 
Dy=1 
Xrange=ABS( Xmax- Xmin) 
Yrange=ABS( Ymax-Ymin) 
WINDOW Xmin, Xmax , Ymin, Ymax 
AXES Dx ,Dy,Xmin,Ymin,1!, 1 
CLIP Cie 
| 
{LABEL PLOT 
CSIZE 4385 
LDIR @ 
LORG S 
MOVE Xmint+.5#Xrange,!.1#Yrange 
TitleS=TRIMS( Titles) 
LABEL Titles 
'LABEL HORIZONTAL AXES 
LDIR @ 
LORG’S 
FOR J=Xmin TO Xmax STEP Dx 
MOVE J,-.033#Yrange 
CSIZe 4.6 
LABEL J 
NEXT J 
MOVE Xmint.5#Xrange,-.12*Yrange 
LORG S 
G5) 76.4 36 
XlabelS=TRIMS( Xlabel $ ) 
LABEL XlabelS 
ILABEL VERTICAL AXES 
LORG 8 
FOR J=@ TO Ymax STEP Dy 
CSTZE «4 2.6 
MOVE Xmin-.@25*Xrange, J 
rales 7 ee 
CSI ZBa2 
MOVE Xmin- .0025* Xrange, J+.@3* Yrange 
LABEL J 
NEXT J 
LDIR PI7Z 
LORG §& 
MOVE Xmin-.1*Xrange,.5*Yrange 
CSiWZEs4 a 
YlabelS=TRIMS( YlabelS) 


ote 


LABEL Ylabel$ 
FOR L=5 TO Ind-1! STEP 5 
X= | 
WHILE Exc_att(X,L)>@ 
Xvar(X)=Exc_att(X,L) 
Range( X)=Range( X-1)+2 
X=X+ | 
END WHILE 
IF INT(L/2)=L/2 THEN 
PINE wee 
EESE 
LINE TYPE 3 
ENO IF 
Freq=L+ 100 
CALL Semi_log( Freq, Xvar(#),X-1,Range(#),X-1) 
NEXT L 
PRINT “HIT CONTINUE TO CONTINUE" 
PAUSE 
! 
GCLEAR 
OUTPUT KBD;"K"; 
LINPUT “WOULD YOU LIKE TO MAKE ANOTHER RUN (CENTER Y OR N)?" XS 
Again=FNYes( X$) 
Mess_up=] 
WHILE Mess_up=1 
SELECT Again 
Choe! 

LINPUT "WOULD YOU LIKE TO CHANGE A VARIABLE (ENTER Y OR N)?" XS 
New_va=FNYes( X$) 
Ano_change=! 

WHILE Ano_change=! 
SELECT New_va 


CASE®! 
PRINT “ VARIABLE CURRENT VALUE" 
PRINT USING "K";"1 TEMPERATURE P lemp. CELSIUS: 


PRINT USING "K";"2 ATOMOSPHERIC PRESSURE ",Atom_pres,"mb”" 
PRINT USING “K";"3 WATER VAPOR PRESSURE “,Wat_pres, mb” 


PRINT " " 
PRINT USING "K";"4 ANTENNA DIAMETER " Ant_diam," m." 
PRINT USING "K"3;"5 PULSE LENGTH " Pulse,” ms” 


PRINT USING "K";"6 POWER TRANSMITTED “,Pow_trans," WATTS" 
PRINT USING: (seni MCOPHERICSPROM@ILE ° Pret vie 
PRIN 

PRINT "ENTER THE NUMBER OF THE VARIABLE YOU WISH TO" 

PRINT “CHANGE” 


INPUT Var 
SELECT Var 
GASe | 
INPUT "TEMPERATURE=" , Temp 
CASE 2 


a 


INPUT “ATMOSPHERIC PRESSURE IN mb=",Atom_pres 


CASE 3 
INPUT “WATER VAPOR PRESSURE IN mb=",Wat_pres 
CASE 4 
INPUT “ANTENNA DIAMETER IN m.=",Ant_diam 
CASE 5 
INPUT “PULSE LENGTH IN ms=",Pulse 
CASE 6 
INPUT “POWER TRANSMITTED IN WATTS=" ,Pow_trans 
CASE 7% 
PRINT “ 1 ==> Ct*Z PROFILE OF Z*(-1.16)" 
PRINT " FROM WALTERS/KUNDEL 1981" 
PRINT " 2 ==> SAME AS ONE BUT WITH FACTOR" 
PRINT " OF TWO TO APPROXIMATE LOOKING" 
PRINT " UP A THERMAL PLOOM" 
PRINT “ 3 ==> Ct°Z PROFILE OF EXPC-200l2 7 
PRINT " FROM WALTERS/KUNDEL 1981." 


PRINT “ 4. ==> Ci Z PROFEEE OF 2 C=47-30— 
INPUT "ENTER NUMBER OF DESIRED PROFILE” ,Profile 
IF Profile=4 THEN 
INPUT “HEIGHT OF INVERSION IN METERS=" Inver 
END IF 
CASE "EESE 
PRINT Var,"IS NOT ONE OF THE OPTIONS" 
END SELECT 
LINPUT “MADE ANOTHER CHANGECENTER Y OR N)?" XS 
New_vacFNYes( XS) 
Mess_up=2Z 
GASe 2 
Mess_up=2 
Ano_change=Z 
CASE EESE 
PRINT "YOUR CHOICES WHERE Y OR N !!t!!" 
Ano_change=! 
LINPUT “CHANGE A VARIABLE (ENTER Y OR N)?" | X& 
New_va=FNYes( X$) 
END SELECT 
END WHILE 
CASE  Z 
Mess_up=2 
CASE EESE 
PRINT “YOUR CHOICES WHERE Y OR N J!!/Et1" 
Mess _up=! 
LINPUT “MAKE ANOTHER RUN (ENTER Y OR N)?" XS 
Again=FNYes( X$) 
END SELECT 
END WHILE 
END WHILE 
PRINT “THAT'S ALLY FOLKS 
end 


aS 


! 

ICALCULTE THE ATTENUATION 

sub Attenuation(Atom_pres ,Atten, Freq, Temp, Wat_pres) 

!1'THIS SUBPROGRAM CALCULATES THE ATTENUATION OF SOUND 

fIN AIR BASED UPON EQUATIONS IN NEFF 1975 

INPUT ATMOSPHERIC PRESSURE IN MILLIBARS 

FREQUENCY OF SOUND WAVE IN HERTZ 

TEMPERATURE IN DEGREES CELCIUS 
WATER-VAPOR PRESSURE IN MILLIBARS 

OUTPUT ATTENUATION IN 1/METERS 


VARIABLES 
Atom_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 
Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 


l 

| 

l 

| 

| 

| 

| SUBPROGRAM ATTENUATION. 

i Att_max VARIABLE IN SUBPROGRAM ATTENUATION. IT IS THE 
| ATTENUATION AT THE FREQUENCY OF THE MAXIMUM 

l ATTENUATION FOR THE INPUT CONDITIONS. 

| 48 VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
| 

i 

| 

| 

i 

i 

{ 


RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 

Fmax FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUSPROGRAM 
ATTENUATION. 

Freq INPUT FREQUENCY OF ECHOSOUNDER. 

H VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 

Pstar VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 

Temp INPUT TEMPERATURE IN DEGREES CELSIUS. 

Tstar INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 


IN SUBPROGRAM ATTENUATION. 

Wat pres ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
| OPERATOR. 
H=100*Wat_pres/Atom_pres 
Tstar=(1.8*Temp+492 )/519 
Pstar=Atom_pres/1014 
Fmax=( 104+6600*H+44400*#H#*H)*Pstar/Tstar*.8 
Att_max=.@078*#Fmax*Tstar*(-2.5)#*#EXP(7.77#(1-1/Tstar) ) 
F=Freq/F max 
Atten=(Att_max/304.8)#((.18#F)*24(Z24F#F/( 14F*#F))0°2)*.5 
Atten=(Atten+1.74E-10*#Freq#Freq)/4. 35 
SUBEND 
j 
def FNYes( XS) 


['THIS FUNCTION INTERPRETS THE OPERATORS RESPONSES TO YES NO QUESTIONS 
’ 


- 


| INPUT XS 
{ 
| OUTPUT FNYes 


oN 


{VARIABLES 
! Temps VARIABLE STRING USED IN FUNCTION YES. 
{ =6XS STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 
RESPONSE TO YES OR NO QUESTION. 
DIM TempS[ 1] 
Temp$S(1,tJ=TRIMS( XS) 
SELECT Temps 
CASE “You yy: 
RETURN | 
CASE “NN” “nn” 
RETURN 2 
Caco 
RETURN 1 
CASE Reise 
RETURN -2 
ENO sere 
FNEND 
| 
sub Semi_log(Freq,Xvar(+#),L,Range(*#) ,J) 
| 
! THIS SUBROUTINE MAKES A SEMI-LOG PLOT OF DATA PASSED FROM THE 
!MAIN PROGRAM. 
l 
| INPUT 
! X AND Y VALUES TO BE PLOTTED 
OUTPUT SEMI-E0GePEOT 


VARIABLES 
IL SECOND ORDER INDEX FOR ASSORTED LOOPS. 


| 
Poe FIRST ORDER INDEX FOR ASSORTED LOOPS. 
| 


Xvar(#) VALUE OF ARRAY OF X VALUES TO BE PLOTTED IN SUBPROGRAM Seni_ 


CLIP ON 
FOR J=)" (Os esiTees 
PLOT Xvar(J),LGT( Range(J)) 
NEXT J 
(DIR eels 
LORG 2 
CSl7E325 
LINE TYPE | 
MOVE Xvar(J-1),LGT( Range( J-1)) 
LABEL USING@ Kes treea, az 
SUBEND 


LOO 


APPENDIA C 
COMPUTER PPOGRAM 3 


!FULLER, ROBERT PROG_3 EFFICIENCY 
| 
{10 SEP 8S 
| 
| 
leeeeeee PURPOSE se enanae 
| 
| THIS PROGRAM WILL TAKE PARAMETERS OF AN ECHOSOUNDER AND 
[ESTIMATE THE RANGE AS A FUNCTION OF EFFICIENCY OF THE THE TRANSDUCER. 
!1THE FOLLOWING INPUTS ARE REQUIRED: 
| ATMOSPHERIC DATA 
1 )ATOMOSPHERIC PRESSURE in millibars 
Z)SELECT ONE OF FOUR ATMOSPHERIC PROFILES FOR THE TEMPERATURE 
STRUCTURE PROFILE. 
a)FOR ONE PROFILE THE HEIGHT OF THE INVERSION LAYER IS ALSO INPUT 
3) TEMPERATURE IN DEGREES CELSIUS 
4)WATER VAPOR PRESSURE IN millibars 


ECHOSOUNDER DATA 

S)ANTENNA DIAMETER IN METERS 

6) FREQUENCY 

7)POWER TRANSMITTED BY ECHOSOUNDER IN WATTS 
8)PULSE LENGTH OF THE TRANSMITTED ACOUSTIC ENERGY 


THE PROGRAM OUTPUTS THE FOLLOWING GRAPH TO AN EXTERNAL PLOTTER: 
1 )RANGE(m.) VERSUS EFFICIENCY OF ECHOSOUNDER 


tee2eaVARTABLES#a44# 


Again TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 
WISHES TO MAKE ANOUTHER RUN. 

Ano_change TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 
WISHES TO CHANGE THE VALUE OF INPUT DATA BEFORE MAKING 
ANOUTHER RUN. 

Ant_area ECHOSOUNDER ANTENNA AREA CALCULATED FROM INPUT OF 
ANTENNA DIAMETER. 

Ant_diam INPUT OF ECHOSOUNDER ANTENNA DIAMETER IN METERS. 

Atom_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 


Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 
SUBPROGRAM ATTENUATION. 
Bn BANDWIDTH OF RECEIVER NEEDED TO RECEIVE SIGNALS 


DOPPLER SHIFFED BY THREE METER PER SECOND 
VERTICAL VELOCITIES. 


C VELOCITY OF SOUND CALCULATION AT ZERO CELSIUS. 

C1 VARIABLE USED TO SCALE INVERSION HEIGHT IN 
TEMPERATURE STRUCTURE PROFILE THREE. 

C3 VARIABLE USED IN CALCULATING THE FOURTH CteZ PROFILE. 
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Cnez(*) ARRAY OF VALUES OF THE ACOUSTIC REFRACTIVE INDEX 
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o—_— oS co — com —_ — — — —— em sui a — ¢—— —— ——_ — — ca — — — — — sa — — — so o— — o— — sao — — — —— o— saa (or ee ee eee ——— cit Sot 


Ctez 


Cve2 


Ox 

Dy 
Epsilon 
Eff 

Eff 


ExClawt 
E 


Fmax 


Freq 
G 
H 
i 


IND 
Interval 
Inver 

J 

K 

2 
Mess_up 


N 
New_va 


Noise 


PARAMETER. CALCULATED BASED ON SELECTION OF 

PROFILE FOR TEMPERATURE STRUCTURE PARAMETER 

AND ASSUMED PROFILE FOR THE VELOCITY STRUCTURE 
PARAMETER. 

VALUE OF THE TEMPERATURE STRUCTURE PARAMETER. 

THE OPERATOR SELECTS A PROFIEE FROM SEVERGE PRO TpED- 
VALUE OF THE VELOCITY STRUCTURE PARAMETER. 

VALUES BASED ON CALCULATION USING ASSUMED DISSAPATION 
RATS 

STEP SIZE FOR X AXIS FOR VARIOUS PLOTS. 

STEP SIZE FOR Y AXIS FOR VARIOUS PEOTse 

DISSAPATION RATE USED IN CALCULATION OF Cvez. 
TRANSMISSION EFFICIENCY OF ECHOSOUNDER. 

EFFICIENCY OF ECHOSOUNDER WHEN RECEIVING ACOUSTIC 
BACKSCATTER. 
EXCESS “ATTENUATION” AT GIVEN RANGE. 
VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 

RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 
FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
ATTENUATION. 

INPUT FREQUENCY OF ECHOSOUNDER. 

ANTENNA EFFECTIVE APERATURE FACTOR. 

VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 
MAIN INDEX FOR LOOP TO CHECK RETURNED POWER UNTIL 
POWER WAS LESS THAN MINIMUM DETECTABLE SIGNAL. 

INDEX USED FOR LOOP FOR DIFFERENT EFFICIENCIES 

MAXIMUM EFFECTIVE-SCATTERING VOLUME THICKNESS. 

HEIGHT OF INVERSION LAYER. 

FIRST ORDER INDEX FOR ASSORTED LOOPS. 

WAVENUMBER 

SECOND ORDER INDEX FOR ASSORTED LOOPS. 

TEST VALUE TO CHECK RESPONSES WHEN OPERATOR IS INPUTTING 
RESPONSES. 

TEST VALUE IN CALCULATION OF COHERENCE LENTH IN CALCULATION 
OF EXCESS ATTENUATION. 
TEST VALUE IN CHECKING IF OPERATOR WISHES TO CHANGE 

A VARIABLE BEFORE A NEW RUN. 
ASSUMED MINIMUM DETECTABLE SIGNAL. 


Pow_back(*) POWER BACKSCATTERED FROM GIVEN RANGE. 


Pow _ret(#) 
Pow_trans 
Profile 
Pstar 


Pulse 


R 

RI 
Range( * ) 
Remainder 


POWER BACKSCATTERED TO ECHOSOUNDER FROM GIVEN RANGE. 
INPUT OF POWER SUPPLIED TO TRANSDUCER OF ECHOSOUNDER. 
OPERATOR INPUT OF Ctez PROFILE FROM AVAILABLE PROFILES. 
VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 

TRANSMITTED PULSE LENGTH IN MILLISECONDS. INPUT BY 
OPERATOR. 

RANGE FOR INPUT TO INTEGRAL OF EXCESS ATTENUATION. 
VARIABLE USED IN SCALING INVERSION HEIGHT FOR PROFILE THREE. 
ARRAY OF RANGE VALUES. 

REMAINDER OF MODULO FUNCTION USED TO DECREES THE 
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NUMBER OF PASSES THROUGH THE INTEGRATION FOR EXCESS 


ATTENUATION. 
Rage VALUE OF RANGE IN INTEGRAL FOR EXCESS ATTENUATION. 
Rho CORRELATION LENGTH USED IN CALCULATION OF EXCESS 
ATTENUATION. 
Sigma FRACTION OF ENERGY BACKSCATTERED FROM GIVEN RANGE. 


Speed_sound SPEED OF SOUND IN AIR AT INPUT TEMPERATURE. 
Sumpow_ back 


if INPUT TEMPERATURE IN DEGREES KELVIN. 

Temp INPUT TEMPERATURE INe DEGREES CELSIUS. 

Temp$ VARIABLE STRING USED IN FUNCTION YES. 

Tstar INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 
IN SUBPROGRAM ATTENUATION. 

Var USED TO SEPECT THESVeRI ABLE TRESOPERATOR WISHED TO 


Wat_pres ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 


OPERATOR. 

x THIRD ORDER INDEX USED IN VARIOUS LOOPS. 

XS STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 
RESPONSE TO YES OR NO QUESTION. 

Xmax VALUE OF LARGEST VALUE OF X FOR VARIOUS PLOTS. 

Xmin VALUE OF “SMALEBST VALUE OF X FOR VARIOUS PLOTS. 

Xrange VALUE OF RANGE OF X VALUES FOR VARIOUS PLOTS. 

Ylabel$ LABEL ON Y AXIS PASSED TO SUBPROGRAM Semi_log FOR PLOTTING. 

Ymax LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

Ymin SMALLEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

Yrange RANGE OF Y VARIABLES FOR VARIOUS PLOTS. 

we INTERMEDIATE VALUE IN CALCULATION OF EXCESS ATTENUATION. 


l 
l 
l 
l 
l 
l 
l 
| 
l 
l 
! CHANGE BEFORE MAKING ANOTHER RUN. 
l 
l 
’ 
l 
| 
l 
| 


| 

(DECLARE VARIABLES 

INTEGER I 

DIM Pow_back( 1500) ,Pow_ret( 150@) , Range( 1500) 

DIM Cne2( 1500) ,Ran({ 30) 

| 

PLOTTER IS 705,"HPGL" 

LINE TYPE 1 

! INPUT ATMOSPHERIC DATA 

INPUT "ENTER TEMPERATURE IN DEGREES CELCIUS" , Temp 

INPUT “ENTER ATMOSPHERIC PRESSURE IN MILLIBARS" ,Atom_pres 
INPUT “ENTER WATER VAPOR PRESSURE IN MILLIBARS" ,Wat_pres 

| 

| 

1! INPUT ECHOSOUNDER DATA 

INPUT “ENTER ANTENNA DIAM IN METERS" ,Ant_diam 

INPUT “ENTER ECHOSOUNDER FREQUENCY IN HERTZ" Freq 

INPUT “ENTER ECHOSOUNDER PULSE LENGTH IN MILLISECONDS=" ,Pulse 
INPUT “ENTER POWER TO TRANSMITTER IN WATTS" ,Pow trans 
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G=.40@ 
| 
| 


'ANTENNA EFFECTIVE APERATURE FACTOR 


Hj ee) THE PROFILE FOR THE TEMPERATURE STRUCTURE 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
1 

PRINT 
PRINT 
Mess_ 
WHILE 


“YOU HAVE A CHOICE OF TEMPERATURE STRUCTURE PROFILES” 

“PROFILE TO BE USED.” 

“YOUR SELECTIONS ARE” 

‘i 1 A TEMPERARURE STRUCTURE PROFILE BASED ON DATA AS “ 
i PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 WHICH " 
i. GIVES A HEIGHT TO THE -1.16 PROFILE" 

a Z THE TEMPERATURE STRUCTURE PROFILE ABOVE BUT" 

" WITH A FACTOR OF TWO TO APPROXIMATE LOOKING" 

: UP A THERMAL PLOOM" 

" 3 A TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS“ 
. PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 FOR" 

"3 NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT " 

7 IS EXP(-.@@1*HEIGHT ABOVE 65 METERS AND HEIGHT" 

" TO THE -1.46 BELOW 65 METERS" 

x 4 A TEMPERATURE STRUCTURE PROFILE BASED ON A” 

Gj HEIGHT TO THE -4/3 " 
up= 1 

Mess _up=1 


INPUT “ENTER THE DESIRED PROFILE (1 OR 2 °OR 3 Ohe4>  -Protile 
IF Profile=4 THEN 

INPUT “ENTER THE HEIGHT OF THE INVERSION IN METERS", Inver 
END IF 
IF Profile=1 THEN 


Mess_up=0 
Suse 
IF Profile=2 THEN 
Mess_up=0 
ELSE 
IF Profile=3 THEN 
Mess_up=@ 
ELSE 
IF Profile=4 THEN 
Mess_up=@ 
ELSE ; 
PRINT Profile," WAS NOT ONE OF THE OPTIONS! !!1!" 
Mess_up=] 
END IF 
END IF 
END IF 


END IF 
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ENO WHILE 
OUTPUT KBD; "KK; 
| 
Again=1 
WHILE Again=1 
FOR J=1 TO I 
Pow_back(J)=@ 
Pow_ret(J)=@ 
Range( J)=@ 
Cne2(J)=@ 
NEXT J 
| 
!CONVERT TEMPERATURE TO KELVIN 
T=Temp+273 
| 


!'CALCULATE SPEED OF SOUND BASED ON INPUT TEMPERATURE( CELCIUS ) 
Speed_sound=2@. 05*(7T)*.5 

PCALECUEANE TheeSreee GF SOUND AT @ GEGREES CEECIUS 
C=20.05+*273* .5 

| 


Le) !INDEX FOR INNER LOOP (POWER RETURNED TO NOISE) 
Ind=! !INDEX FOR OUTER LOOP (EFFICIENCY) 
Pow_back(@)=@ !INITIALIZE VARIABLE FOR POWER BACKSCATERRED 


Ant_area=PI+*(Ant_diam/2)*°*2 !'ANTENNA AREA 
Interval=( Speed_sound*Pulse+*#!.E-3)/2 
lINITIALIZE PLOTTING PARAMETERS SO CAN COMPARE DURING RUN TO SCALE 
HT eeLO a pS 
Ymin=1000@ 
Ymax=@ 
1CALCULATE RETURNED SIGNAL POWER UNTIL IT IS LESS THAN BACKGROUND 
'FOR EACH FREQUENCY. 
Range( @)=@ 
POReet7=.05 10 .S STEP .95 
Bn=2*#Freq*( 1-1/(6/Speed_sound+1))!BANDWIOTH FOR 3M/S VERTICAL VELOCITY 
Noise=1.38E-23#Bn#(T)+2.E-14!MINIMUM DETECTABLE SIGNAL AS 
JOHNSON NOISE(NEGLIGABLE)+ ESTIMATED BACKGROUND 
ICALCULATE THE ATTENUATION COEFFICIENT BASED ON ATMOSPHERIC DATA 
CALL Attenuation( Atom_pres,Atten,Freq, Temp, Wat_pres) 
K=Z*PI*Freq/Speed_sound 1! WAVENUMBER 
I=] 
Sumpouw_back =@ 
Kereta 
Range( I )=Range(I-1)+2 
SELECT Prof ve 
Criss 
'THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
1'AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398. 
'WAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 


}O5 


! TULAROSA BASIN, NEW MEXICO. 
!'AN ADDITONAL FACTOR OF 1.8 WAS INCLUDED AS PER DIFFERENT 
!AVERAGING TIME. 
Cte2=2.12*#Range(I)*(-1.16) 
EAoe 
'THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
!AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398. 
!'WAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
!TULAROSA BASIN, NEW MEXICO. 
!AN ADDITONAL FACTOR OF 1.8 WAS INCLUDED AS PER DIFFERENT 
!AVERAGING TIME AND AN ADDITIONAL FACTOR OF TWO WAS 
!INCLUDED TO APPROXIMATE LOOKING UP A THERMAL PLOOM 
Cte2=2+#*2.12*Range(I)*(-1.16) 
CASESS 
!'THIS TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS" 
IPRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 FOR" 
!'NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT “ 
11S EXP(-.@@1*«HEIGHT ABOVE 65S METERS AND HEIGHT" 
(TO THE -1746 SELOU ESeyETERS 
IF Range(1I)<65 THEN 
Cte2=75.5*Range(I)*(-1.46) 
ELSE 
Cte2=3.66E-Z2*EXP(-.001*Range(I)) 
END IF 
CASE 4 
!'CALCULATE THE TEMPERATURE STRUCTURE FACTOR BASED ON ASSUMPTION 
THAT IT’S PROPORTIONAL TO HEIGHT TO THE -4/3 POWER 
THAT THE SURFACE VERTICAL HEART FLUX 15290395 €.m.7sec 
!'THE DEPENDENCE WITH HEIGHT ABOVE THE INVERSION I PULLED FROM 
thy LEFT EAR 
{EQUATION FROM NEFF, 1975 
SELECT Range(I)/Inver 
GASE-<.3 
C3=((.024)*(T)*( .667)) 
CteZ=C3*#( Range(I))*(-1.33) 
Css eel ey) 
Ctez=Ctez 
C1=CteZ2 
R1=Range( I) 
CASE 170 1.3 
CteZ=10*°( (LGT( C3*Range( I) )-LGT( C3#R1))+LGT( C1) ) 
R2Z=Range( I) 
C2=Cte2 
CASE ELSE 
CteZ=C3#( Range(I)*(-1.33)-RZ°(-1.33))+C2 
END SELECT 


END SELECT 

!CALCULATE VELOCITY STRUCTURE FACTOR. FORMULA FROM GAYNOR 77 
Enpsilon=( .2866/Range( 1) )#*(1+.07#( Range(I))*¢(.6))*¢1.5) 
Cve2Z=2*Epsilon*( .667) 
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{CALCULATE ACOUSTIC REFRACTIVE INDEX STRUCTURE FACTOR. 
!FORMULA FROM TETARSKI 
Cne2(1I )=(Cte2/(2.98E+S) )+( Cve2/(C#C) ) 
{CALCULATE THE FRACTION OF POWER BACKSCATTERED FROM INTERVAL. 
{EQUATION FROM NEFF, 1975 
Sigma=( .Q@039#(K*(1/3))#Cte2)/(T)%*2 
{ 
{CALCULATE THE EXCESS ATTENUATION. EQUATION FROM CLIFFORD, 1980 
|THE EXCESS ATTENUATION IS Exc_att 
IF I1<S@ THEN 
Remainder =@ 
EVoe 
IF I<tS@ THEN 
Remainder=Range(I) MODULO 1@ 
ESE 
Remainder=Range(I) MODULO 20 
END IF 
END IF 
IF Remainder=@ THEN 
Rho=@ 
L=@ 
H=@ 
Rge=Range( I) !CONSTANT IN INTEGRAL 
R=@ 
FOR J=@ TO 2#I 
F=Cne2( INTC J/2+1)) 
F=Fe( 1-R/(Rge)*(1.67)4¢(R/( Rage) )*(1.67)) 
IF J>@ THEN 
IF J<2*I THEN 
IF INT(J/2)=J/2 THEN 
L=L+F 
F=@ 
Bese 
H=H+F 
F=Q@ 
ENO IF 
END IF 
END IF 
Rho=RhotF 
R=R+1 
NEXT 
Rho=Rhot4#L4+2+*H 
Rho=(((Rho#.33)*#K#K#1.46)°(-.6)) 
N=(Ant_diam/Rho)*2 
IF N<=! THEN 
Ze=1/¢(14+N) 
BESE 
Ze=1.S/(14N) 
!'STEP OF 1.S5==>SEE CLIFFORD 1980 
END IF 


Exc_att=Ze+#Ze 
ELSE 
Exc_att=Exc_att 
END IF 
I 
ICALCULATE THE POWER BACKSCATTERED 
Pow_back(1I)=(Pow_trans#Ef  f-Sumpou_back )*EXP(-Atten#Range(I) ) 
Pow_back(I)=Pow_back(I)*Interval#Exc_att*#Sigma 
Sumpow_back=Sumpow_back +Pow_back( 1) 
ICALCULATE THE POWER BACKSCATTERED TO THE ANTENNA 
Pow_ret(1)=Pow_back(I)*#EXP(-Atten*Range( I) )*Ant_area#G#Ef f/Range(I)*Z 
PRINT "RANGE=",Range(I ) 
PRINT “POWER RETURNED=" ,Pow_ret(I) 
I=I+1 
UNTIL Pow_ret(I-1!1)>@ AND Pow_ret(I-1)<=Noise 
Ran( Ind)=Range( I-2Z) 
IF Range(I-2)>Ymax THEN 
Ymax=Range(I-Z ) 
END IF 
IF Range(I-2)<Ymin THEN 
Ymin=Range(1I-Z) 
END IF 
Ind=I nd+1 
Again=1| 
FOR J=! TO Il 
Pow_back(J)=@ 
Pow _ret(J)=@ 
Range( J)=@ 
Cne2Z(J)=@ 
NEXT J 
NE Bat 
OUTPUT KBD;"K"; 


PRINT “INPUT CONDITIONS” 

{ 

PRINT " “ 

PRINT USING "K";"TEMPERATURE= " jiieipe CELCIUS” 
PRINT USING “K";"ATMOSPHERIC PRESSURE= ",Atom_pres, “mb” 
PRINT USING "K"; "WATER VAPOR PRESSURE= " Wat_pres,"mb" 
PRINT USING “K"( {SUb SE Genet h= ", Puligey moe 
PRINT USING "K"; "ANTENNA DIAMETER= Werenigns op w= (Je 1 am 
PRINT USING “K";"ECHOSOUNDER FREQUENCY= ",Freq,"H2"“ 
PRINT USING “K"; "POWER TRANSMITTED= " Pou_trans,”" WATTS” 
PRINT “TEMPERATURE STRUCTURE PROFILE USED ",Profile 


IF Profile=4 THEN 
PRINT “INVERSION HEIGHT ", Inver 


END IF 


SELECT Profile 
CASE I 


1OS 


PRINT 
PRINT 
PRINT 
CASE. Z 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
EASE 3 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
CASE 4 
PRINT 
PRINT 


ENO SELECT 
| 


A TEMPERARURE STRUCTURE PROFILE BASED ON DATA” 
AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 298 " 
WHICH GIVES APHETGHT TO THE -1.1G6 PROFILE” 


A TEMPERARURE STRUCTURE PROFILE BASED ON DATA” 
AS PRESENTED IN WALTERS/KUNDEL 1398! PAGE 398 " 
WHICH GIVES A HEPGHT TO THE -1.16 PR@F ELE” 

BUT WITH A FACTOR OF TWO TO APPROXIMATE" 
LOOKING UP A THERMAL PLOOM" 


A TEMPERATURE STRUCTURE PROFILE BASED ON DATA" 
AS PRESENTED IN WALTERS/KUNDEL 1381! PAGE 3938" 
FOR NIGHT CONDITIONS. THE DEPENDENCE WITH" 
HEIGHT 1S EXP(-.@O1*HEIGHT ABOVE "GS METERS” 
AND HEIGHT TO THE -1.46 BELOW 65 METERS" 


A TEMPERATURE STRUCTURE PROFILE BASED ON A" 
HEIGHT TO THE -4/3 " 


PRINT USING "K";"THE MINIMUM DETECTABLE SIGNAL WAS SET AT" ,Noise,"WATTS" 


| 
PRINT “ “ 
' 


PRINT "HIT CONTINUE TO CONTINUE” 


PAUSE 


OUTPUT KBD;"K"; 


GRAPHICS ON 


VIEWPORT 15,120,10,70 


IPLOT RANGE VERSUS FREQUENCY 


Xmin=@ 
Xmax=.5 


Ymin=PROUND( Ymin-5@, 2) 
Ymax =PROUND( Ymax+5@, Z) 


Dx=.@5 
Dy=100 


Xr ange=Xmax- Xmin 
Yrange=Ymax- Ymin 

WINDOW Xmin, Xmax , Ymin, Ymax 
AXES Ox Oy . Xmin- Yuin.) 


EETP ORE 

| 

feneel PLOT 
GSIZE -4...6 
LDIR @ 

LORG 5 


FOS 


MOVE .5*Xrange, Ymax*1.@5 
LABEL “RANGE VERSUS TRANSDUCER EFFICIENCY" 
!LABEL HORIZONTAL AXES 


LDIR @ 
LORG 5 
FOR J=@ TO Xmax STEP Dx 
CSI ZE 42.5 
MOVE J,Ymin-.@5* Yrange 
CABELL J 
NEXT J 
MOVE .5#Xrange,Ymin-.1#Yrange 
CSIZE 4226 


LABEL “EFFIGRENGY 
{LABEL VERTICAL AXES 
LORG 8 
FOR J=Ymin TO Ymax STEP Dy 
CSIZE 4465 
MOVE Xmin-.@0025*Xrange, J 
LABEL USING "K";J 
NEXT J 
LDOIR PI/2Z 
LORG 6 
MOVE Xmin-.15*Xrange, .5*Yrange+Ymin 
CSIZEV4...5 
LABEL “RANGE (METERS)" 
CLIP ON 
Ind=} 
FOR Eff-—205 05-5 cTEP 2.05 
PLOT Eff ,Ran( Ind) 
Ind=Ind+ 1 
NEXT Ett 
I 
PRINT “HIT CONTINUE TO CONTINUE" 
PAUSE 
GCLEAR 
OUTPUT KBD;"K"; 
' 
LINPUT “WOULD YOU LIKE TO MAKE ANOTHER RUN (ENTER Y OR N)?" XS 
Again=FNYes( X$) 
Mess_up=] 
WHILE Mess_up=! 
SELECT Again 
CASE 1 
LINPUT “WOULD YOU LIKE TO CHANGE A VARIABLE (ENTER Y OR N)?" XS 
New_va=FNYes( XS) 
Ano_change= 1 
WHILE Ano_change=!1 
SELECT New_va 
CASE 1 
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PRINT " VARIABLE CURRENT VALUE" 
PRINT USING "“K";"1 TEMPERATURE * Tempo, “CELSIUS” 
PRINT USING "“K";"2 ATOMOSPHERIC PRESSURE “ ,Atom_pres, "mb" 
PRINT USING "“K";"3 WATER VAPOR PRESSURE “, Wat_pres, "mb" 
PRINT * 
PRINT USING "“K";"4 ANTENNA DIAMETER " Ant_diam," m." 
PRINT USING=@°K';"S PULSE LERGTH * Pulse.” cms” 
PRINT USING "K";"6 POWER TRANSMITTED " ,Pow_trans," WATTS" 
PRINT USUNG = "K°:"7  AMMOSPHERIC PROFILE “,Protitle 
PRIN = —= 
PRINT “ENTER THE NUMBER OF THE VARIABLE YOU WISH TO" 
PRINT “CHANGE” , 
INPUT Var 
SELECT Var 
CASE 1 
INPUT “TEMPERATURE=" , Temp 
CASE 2 
INPUT “ATMOSPHERIC PRESSURE IN mb=",Atom_pres 
GAce. 3 
INPUT “WATER VAPOR PRESSURE IN mb=",Wat_pres 
CASE 4 
INPUT "ANTENNA DIAMETER IN m.=",Ant_diam 
GASES 
INPUT "PULSE LENGTH IN ms=" ,Pulse 
cass ts 
INPUT “POWER TRANSMITTED IN WATTS=",Pow_trans 
CASEaa 
PRINT = | =S=)%@022 PROFILE OF Z*(-1.16)" 
PRINT “ FROM WALTERS/KUNDEL 1981" 
PRINT " 2 ==> SAME AS ONE BUT WITH FACTOR" 
PRINT " OF TWO TO APPROXIMATE LOOKING" 
PRINT “ UP A THERMAL PLOOM" 
PRINT “" 3 ==> Ct*2 PROFILE OF EXP(-.@01#Z)" 
BRI Ni FROM WALTERS/KUNDEL 1981." 
PRINT “ 4 «==> Ci*Z PROF IeeeOF Z*(-4/3)" 
INPUT "ENTER NUMBER OF DESIRED PROFILE" Profile 
IF Profile=4 THEN 
INPUT “HEIGHT OF INVERSION IN METERS=", Inver 
END IF 
CAS Ew EIESE 
PRINT Var ,"“IS NOT ONE OF THE OPTIONS" 
END SELECT 
LINPUT "MADE ANOTHER CHANGEC ENTER Y OR N)?" XS 
New_va=FNYes( XS) 
Mess _up=2 
CASE 2 
Mess_up=2Z 
Ano_change=2 
EASES IESE 
PRINT "YOUR CHOICES WHERE Y OR N /till" 


Ano_change=!1 
LINPUT "CHANGE A VARIABLE (ENTER Y OR N)?", X$ 
New _va=FNYes( XS) 
END? SEEeGr 
END WHILE 
CASE Z 
Mess_up=2 
CASeeeesr 
PRINT “YOUR CHOICES WHERE Y OR N Jlldti" 
Mess_up=1 
LINPUT "MAKE ANOTHER RUN CENTER Y OR N)?" | X$ 
Again=FNYes( X$) 
END Seeee 
END WHILE 
END WHILE 
PRINT “THAT'S ALL, FOLKS” 
end 
| 
ICALCULTE THE ATTENUATION 
sub Attenuation( Atom _pres,Atten,Freq, Temp, Wat_pres) 
!THIS SUBPROGRAM CALCULATES THE ATTENUATION OF SOUND 
'IN AIR BASED UPON EQUATIONS IN NEFF 1975 
| 
INPUT ATMOSPHERIC PRESSURE IN MILLIBARS 
FREQUENCY OF SOUND WAVE IN HERTZ 
TEMPERATURE IN DEGREES CELCIUS 
WATER-VAPOR PRESSURE IN MILLIBARS 


GUTPUT ATTENUATION IN°1/METERS 


VARIABLES 
Atom_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 
Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 


| 

l 

| 

| 

| 

| 

| SUBPROGRAM ATTENUATION. 

| Att_max VARIABLE IN SUBPROGRAM ATTENUATION. IT IS THE 
ATTENUATION AT THE FREQUENCY OF THE MAXIMUM 
| ATTENUATION FOR THE INPUT CONDITIONS. 

| 
| 
| 
| 
| 
l 

l 
l 

| 


iP VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 

Fmax FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
ATTENUATION. 

Freq INPUT FREQUENCY OF ECHOSOUNDER. 

H VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 

Pstar VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 

Temp INPUT TEMPERATURE IN DEGREES CEbSotuG. 

Tstar INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 


IN SUBPROGRAM ATTENUATION. 
Wat _pres ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 


bi 2 


| OPERATOR. 

H=1@0#Wat_pres/Atom_pres 

Tstar=(1.8*Temp+492)/519 

Pstar=Atom_pres/1@Q14 

Fmax=( 10+6600#H+44400#H#H)*#Pstar/Tstar*.8 
Att_max=.Q@078*#F max*Tstar*(-2.5)#EXP(7.77#(1-1/Tstar) ) 
F=Freq/F max 
Atten=(Att_max/304.8)#((.184F)°2+(2#F#F/( 14+F#F))°2)°.5 
Atten=(Attent!.74E-10#Freq#Freq)/4.35 

SUBEND 

| 

def FNYes( XS) 

{THIS FUNCTION INTERPRETS THE OPERATORS RESPONSES TO YES NO QUESTIONS 
| 

| INPUT X$ 

| 

| OUTPUT FNYes 


!VARIABLES 
! Temps VARIABLE STRING USED IN FUNCTION YES. 
| XS STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 


| RESPONSE TO YES OR NO QUESTION. 
DIM TempS[ 1] 
TempSl1,1]=TRIMS( X$) 
SELECT Temps 
GBS “Y". 3a 
RETURN 1 
CASE NN] ae 
RETURN 2 
GASE ~~ 
RETURN 1 
CASEVELSE 
RETURN -2 
ENDs SELBET 
FNEND 
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APPENDIX D 


COMPUTER PROGRAM 4 


!FULLER, ROBERT PROG_4 ANTENNA SIZE 
ig sep Tes 
{ 
| 
beeeeeee2 PURPOSE RBH RHEE 
! 
! THIS PROGRAM WILL TAKE PARAMETERS OF AN ECHOSOUNDER AND 
'ESTIMATE THE RANGE AS A FUNCTION OF ANTENNA SIZE OF THE THE ECHOSOUNDER. 
!'THE FOLLOWING INPUTS ARE REQUIRED: 
| ATMOSPHERIC DATA 
1 )ATOMOSPHERIC PRESSURE in millibars 
Z)SELECT ONE OF FOURSATMOSPHERIG PROFILES FOR TRE TENPERRmwae 
STRUCTURE PROPIEE. 
a)FOR ONE PROFILE THE HEIGHT OF THE INVERSION LAYER IS ALSO INPUT 
SI TENPERATURE IN DEGREES CEEonES 
4)WATER VAPOR PRESSURE IN millibars 


ECHOSOUNDER DATA 

5 )FREQUENCY 

6)POWER TRANSMITTED BY ECHOSOUNDER IN WATTS 
7)PULSE LENGTH OF THE TRANSMITTED ACOUSTIC ENERGY 


THE PROGRAM OUTPUTS THE FOLLOWING GRAPH TO AN EXTERNAL PLOTTER: 
})RANGE(m.) VERSUS ANTENNA DIAMETER FOR VARIOUS FREQUENCIES 


*e#e2 4 VARI ABLESH# 444% 


WISHES TO MAKE ANOUTHER RUN. 

Ano_change TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 
WISHES TO CHANGE THE VALUE OF INPUT DATA BEFORE MAKING 
ANOUTHER RUN. 

Ant_area ECHOSOUNDER ANTENNA AREA CALCULATED FROM INPUT OF 
ANTENNA DIAMETER. 

Ant_diam INPUT OF ECHOSOUNDER ANTENNA DIAMETER IN METERS. 

Atom_pres INPUT OF ATMOSPHERIC SP RESSURE SI Nem. 


Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 
SUBPROGRAM ATTENUATION. 
Bn BANDWIDTH OF RECEIVER NEEDED TO RECEIVE SIGNALS 


DOPPLER SHIFFED "Sy eToREEsMEreR PERBSECONG 
VERTICAL SVELOCITIES. 


C VELOCITY OF SOUND CAECULATION AT ZEROVCEESTUS- 

Cl VARIABLE USED TO SCALE INVERSION HEIGHT IN 
TEMPERATURE STRUGMURE PROF IEE STRReE- 

C3 VARIABLE USED IN CALCULATING THE FOURTH GteZ PROFICE- 


Cnez(*) ARRAY OF VALUES OF THE ACOUSTIC REFRACTIVE INDEX 


! 
! 
! 
! 
| 
! 
| 
l 
I 
| 
I 
! 
| 
| 
| 
| 
| Again TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 
| 
| 
l 
’ 
| 
| 
| 
l 
l 
PARAMETER. CALCULATED BASED ON SELECTION OF 
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o— cam ‘amp cam cm ¢umh comm cam tam 4 ocamm came came camp 1 == tom _— re eee eee eile ee cee 


Ctez 


Cve2 


Dx 

Dy 
Epsilon 
Eff 

Eff 


Excwvatt 
F 


F max 
Freq 
G 
H 
I 
IND 


Interval 
Inver 


New_va 


Noise 


PROFILE FOR TEMPERATURE STRUCTURE PARAMETER 
AND ASSUMED PROFILE FOR THE VELOCITY STRUCTURE 
PARAMETER. 
VALUE OF THE TEMPERATURE STRUCTURE PARAMETER. 
THETGEERATGR SEPECTS AlLPROFILE FROM SEVER@E PROVIDED. 
VALUE GF TRE VEEGOCITY STRUCTURE PARAMETER. 
VALUES BASED ON CALCULATION USING ASSUMED DISSAPATION 
nanee 
STEP SIZE FOR Xeaas FOR VARIOUS PLOTS. 
STEP “SIZESROR Y AXES FOR VARIOUS PLOTS. 
DISSAPATION RATE USED IN CALCULATION OF Cve2. 
TRANSMISSION EFFICIENCY OF ECHOSOUNDER. 
EFFICIENCY OF ECHOSOUNDER WHEN RECEIVING ACOUSTIC 
BACKSCATTER. 
EXCESS “ATTENUATION” AT GIVEN RANGE. 
VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 
FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
ATTENUATION. 
INPUT FREQUENCY OF ECHOSOUNDER. 
ANTENNA EFFECTIVE APERATURE FACTOR. 
VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 
MAIN INDEX FOR LOOP TO CHECK RETURNED POWER UNTIL 
POWER WAS LESS THAN MINIMUM DETECTABLE SIONAL. 
INDEX USED FOR LOOP FOR DIFFERENT EFFICIENCIES 
MAXIMUM EFFECTIVE-SCATTERING VOLUME THICKNESS. 
HEIGHT OF INVERSION LAYER. 
FIRS ORDER -ENOEX FOR ASSORTES LOOPS. 
WAVENUMBER 
SECOND ORDER INDEX FOR ASSORTED LOOPS. 
TEST VALUE TO CHECK RESPONSES WHEN OPERATOR IS INPUTTING 
RESPONSE 5. 
TEST VALUE IN CALCULATION OF COHERENCE LENTH IN CALCULATION 
OF EXCESS ATTENUATION. 
TEST VALUE IN CHECKING IF OPERATOR WISHES TO CHANGE 
A VARIABLE BEFORE A NEW RUN. 
ASSUMED MINIMUM DETECTABLE SIGNAL. 


Pow_back(*) POWER BACKSCATTERED FROM GIVEN RANGE. 


Pow_ret(*) 
Pow_trans 
Profile 
Pstar 


Pulse 


R 

RI 
Range( # ) 
Remainder 


POWER BACKSCATTERED TO ECHOSOUNDER FROM GIVEN RANGE. 
INPUT OF POWER SUPPLIED TO TRANSDUCER OF ECHOSOUNDER. 
OPERATOR INPUT OF Cte2 PROFILE FROM AVAILABLE PROFILES. 
VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 
TRANSMITTED PULSE LENGTH IN TWRLLISECONOS. |. PNPUT@ByY 
OPE RAIGR: 
RANGE FOR INPUT TO INTEGRAL OF EXCESS ATTENUATION. 
VARIABLE USED IN SCALING INVERSION HEIGHT FOR PROFILE THREE. 
ARRAY OF RANGE VALUES. 
REMAINDER OF MODULO FUNCTION USED TO DECREES THE 
NUMBER OF PASSES THROUGH THE INTEGRATION FOR EXCESS 


eS 


Rge 
Rho 


Sigma 


ATTENUATION. 

VALUE OF RANGE IN INTEGRAL FOR EXCESS ATTENUATION. 
CORRELATION LENGTH USED IN CALCULATION OF EXCESS 
ATTENUATION. 

FRACTION OF ENERGY BACKSCATTERED FROM GIVEN RANGE. 


Speed_sound SPEED OF SOUND IN AIR AT INPUT TEMPERATURE. 


Sumpow_back 


+ 
Temp 
Temp$S 


Tstar 
Var 
Wat_pres 


x 
XS 


Xmax 
Xmin 
Xrange 
Ylabel$S 
Ymax 
Ymin 
Yrange 
Ze 


INPUT TEMPERATURE IN DEGREES KELVIN. 
INPUT TEMPERATURE IN BEGREES GEESTUS. 
VARIABLE STRING USED IN FUNCTION YES. 


INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 
IN SUBPROGRAM ATTENUATION. 

USED TO SELECT THE VARIABLE THE OPERATOR WISHED TO 
CHANGE BEFORE MAKING ANOTHER RUN. 

ATMOSPHERIC WATER PRESSURE IN MILLIBARS. 
OPERATOR. 

THIRD ORDER INDEX USED IN VARIOUS LOOPS. 
STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 
RESPONSE TO YES OR NO QUESTION. 

VALUE OF LARGEST VALUE OF X FOR VARIOUS PLOTS. 

VALUE OF SMALLEST VALUE OFSX FOR VARIOUS PEOns. 
VALUE OF RANGE OF X VALUES FOR VARIOUS PLOTS. 

LABEL ON Y AXIS PASSED TO SUBPROGRAM Semi_log FOR PLOTTING. 
LARGEST VALUE OF Y VeERTRAEEE FOR VaRIGUS PEOTS. 

SMAELEST VALUE OF Y VYARIRBEE FOR VARLOUS SPEOms . 

RANGE OF Y VARIABLES FOR VARIOUS PLOTS. 

INTERMEDIATE VALUE IN CALCULATION OF EXCESS ATTENUATION. 


INPUT BY 


!DECLARE VARIABLES 
INTEGER I 
DIM Pow back (150) ,Pow_ret( 150@) ,Range( 1500) 
DIM Cne2( 1500) , Ran( 50,6) 


PLOTTER 25" 7053 free. 
LINE ee 


! INPUT ATMOSPHERIC DATA 

INPUT “ENTER TEMPERATURE IN DEGREES CELCIUS"”, Temp 

INPUT “ENTER ATMOSPHERIC PRESSURE IN MILLIBARS",Atom_pres 
INPUT “ENTER WATER VAPOR PRESSURE IN MILLIBARS" ,Wat_pres 


!TNPUT ECHOSOUNDER DATA 

INPUT “ENTER ECHOSOUNDER PULSE LENGTH IN MILLISECONDS=" Pulse 
INPUT “ENTER POWER TO TRANSMITTER IN WATTS" ,Pow_trans 

G=.4@ 
Eff=.25 


'ANTENNA EFFECTIVE APERATURE FACTOR 


PSEEECT THE PROEVEE FOR THE TEMPERATURE STRUCTURE 

“YOU HAVE A CHOICE OF TEMPERATURE STRUCTURE PROFILES" 
VPROPIL ES Oe GE SUSED.” 

“YOUR SEEEC TIONS sare. 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
| 

PRINT 


PRINT 


Mess_up=1 
WHILE Mess_up=1 


INPUT “ENTER THE DESIRED PROFILE (1 


1 


A TEMPERARURE STRUCTURE PROFILE BASED ON DATA AS “ 
PRESENTED IN WALTERS/KUNDEL 1981! PAGE 398 WHICH " 
Gives m HElon) 10 THE =!.16 PROFILE” 


THE TEMPERATURE STRUCTURE PROFILE ABOVE BUT" 
WITH A FACTOR OF TWO TO APPROXIMATE LOOKING" 
UP A THERMAL PLOOM" 


A TEMPERATURE STRUCTURE EROFLLE BASED Of@DATA AS” 
PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 FOR” 
NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT " 
IS EXP(-.@01*#HEIGHT ABOVE 65 METERS AND HEIGHT" 
PORinew | 45072 ROU GS METERS” 


A TEMPERATURE STRUCTURE PROFILE BASED ON A” 
HEIGHT TO THE -4/3 " 


OR@Z OR 2a0R 4)” Prof ike 


IF Profile=4 THEN 
INPUT “ENTER THE HEIGHT OF THE INVERSION IN METERS", Inver 


END IF 
IF Profile=1 THEN 
Mess_up=® 
ELSE 
IF Profile=2 THEN 
Mess_up=@ 
ESE 
IF Profile=3 THEN 
Mess_up=@ 
ELSE 
IF Profile=4 THEN 
Mess_up=@ 
ELSE 
PRINT Profile," WAS NOT ONE OF THE OPTIONS! ! II" 
Mess_up=]1 
END IF 
END IF 
END IF 
END IF . 
END WHILE 


OUTPUT KBD;"K"; 


|| / 


| 

| 

Again=1 

WHILE Again=! 

FOR J=1 TO I 
Pow_back(J)=@ 
Pow_ret(J)=@ 
Range( J)=@ 
Cne2(J)=@ 

NEXT J 

I 

ICONVERT TEMPERATURE TO KELVIN 

T=Tempt+273 


ICALCULATE SPEED OF SOUND BASED ON INPUT TEMPERATURE( CELCIUS) 
Speed_sound=Z2@.@05*(T)*.5 
ICALCULATE THE SPEED OF SOUND Al @ BEGREES EEPELUS 
C=20.05* 273° .5 
| 
I=] lINDEX FOR INNER LOOP (POWER RETURNED TO NOISE) 
Ind=1 !INDEX FOR SECOND LOOP (EFFICIENCY) 
Ifreq=! !INDEX FOR THIRD LOOP (FREQUENCY) 
Pow_back(@)=@ lINITIALIZE VARIABLE FOR POWER BACKSCATERRED 
Interval =( Speed_sound*Pulse#1.E-3)/2Z 
l'INITIALIZE PLOTTING PARAMETERS SO CAN COMPARE DURING RUN TO SCALE 
IPLOTS 
Ymin=1000 
Ymax=@ 
ICALCULATE RETURNED SIGNAL POWER UNTIL IT IS LESS THAN BACKGROUND 
'FOR EACH FREQUENCY. 
Range( @)=@ 
FOR Freq=5@0 TO 1500 STEP 50@ 
Sumpow_back =Q 
FOR Ant_diam=.S5 TO 6&.25-Ifreq STEP .25 
Ant_area=PI*(Ant_diam/Z)°Z 
Bn=2*Freq*#( 1-1/(6/Speed_sound+1))! BANDWIDTH FOR 3M/S VERTICAL VELOCITY 
Noise=1.38E-Z23*Bn#(T)+Z.E-14!1MINIMUM DETECTABLE SIGNAL AS 
| JOHNSON NOISE(NEGLIGABLE)+ ESTIMATED BACKGROUND 
ICALCULATE THE ATTENUATION COEFFICIENT BASED ON ATMOSPHERIC DATA 
CALL Attenuation(Atom_pres ,Atten,Freq, Temp, Wat_pres) 
K=Z2*PI*Freq/Speed_sound |WAVENUMBER 
I=] 
REPEAT 
Range( I )=Range(I-1)+2Z 
SELECT VPeGtage 
CASE a 
'THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
l'AS PRESENTED IN WALTERS/KUNDEL 1981 PAGES S9e- 
'WAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
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I TULAROSA BASIN, NEW MEXICO. 
'AN ADDITONAL FACTOR OF 1.8 WAS INCLUDED AS PER DIFFERENT 
|tAVERAGING TIME. 
Cte2=2.12#Range(I)*(-1.15) 
EASE Z 
'THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
'AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398. 
IWAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
!TULAROSA BASIN, NEW MEXICO. 
'AN ADDITONAL FACTOR OF 1.8 WAS INCLUDED AS PER DIFFERENT 
{AVERAGING TIME AND AN ADDITIONAL FACTOR OF TWO WAS 
!INCLUDED TO APPROXIMATE LOOKING UP A THERMAL PLOOM 
CteZ=2#2.12#Range(I)*(-1.16) 
EASEes . 
I'THIS TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS" 
(PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 FOR” 
!'NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT " 
(IS EXP(-.@@1#HEIGHT ABOVE 65 METERS AND HEIGHT" 
1TO THE -1.46 BELOW 65 METERS” 
IF Range(I)<65 THEN 
CteZ=75.5#Range(I)*(-1.46) 
; ESE 
Cte2Z=3.66E-2*EXP(-.001#Range(I)) 
END IF 
CASE 4 
!CALCULATE THE TEMPERATURE STRUCTURE FACTOR BASED ON ASSUMPTION 
! THAT IT’S PROPORTIONAL TO HEIGHT TO THE -4/3 POWER 
| THAT THE SURFACE VERTICAL HEAT FLUX IS .095 C.om./sec 
{THE DEPENDENCE WITH HEIGHT ABOVE THE INVERSION I PULLED FROM 
My sBEFT ERR 
(EQUATION FROM NEFF, 1975 
SELECT Range(I)/Inver 
CASE Y<29 
C3=((.024)+#(7T)*(.657)) 
CteZ2=C3#( Range(I))*(-1.33) 
CASE .9 TO 1 
Cte2=Ctez 
C1=Cte2 
R1l=Range( I) 
CASE, I 10 1235 
Cte2=10°((LGT(C3*Range(I))-LGT(C3*R1))4+LGT(CI)) 
RZ2=Range( I) 
C2=Cte2 
CASE ELSE 
GiLeZ=05*( Ranoe()) (-1.330- R20 -1. 33) )+CZ 
END SELECT 
END, SELECT 
|CALCULATE VELOCITY STRUCTURE FACTOR. FORMULA FROM GAYNOR 77 
Ensilon=(.2866/Range(I))#(1+.07#( Range(I))*(.6))*(1.5) 
Cve2Z=2*Epsilon*( .667) 


Bic 


ICALCULATE ACOUSTIC REFRACTIVE INDEX STRUCTURE FACTOR. 
!FORMULA FROM TETARSKI 

Cne2(1I)=(Cte2/(2.98E+5) )4+( Cve2/(C#C) ) 

!CALCULATE THE FRACTION OF POWER BACKSCATTERED FROM INTERVAL. 
!EQUATION FROM NEFF, 1975 

Sigma=( .@039#(K° (1/3) )#Cte2)/(T)*2 

I 


[CALCULATE THE EXCESS ATTENUATION. EQUATION FROM CLIFFORD, 1980 
!'THE EXCESS ATTENUATION IS Exc_att 
IF I1<S5@ THEN 
Remainder = 
ELSE 
IF I<15@ THEN 
Remainder=Range(I) MODULO 10 
ELSE 
Remainder=Range(I) MODULO 290 
END IF 
END IF 
IF Remainder=0 THEN 
Rho=@ 
L=@ 
H=@ 
Rge=Range( I ) [CONSTANT IN INTEGRAL 
R=@ 
FOR J=@ TO 2+*I 
F=Cne2¢(INT(J/2+1)) 
F=F*( 1-R/(Rge)*(1.67)4+¢(R/( Roe) )*(1.67)) 
IF J>@ THEN 
IF J<2*I THEN 
IF INT(J/2)=J/2 THEN 
L=L+F 
F=@ 
ELSE 
H=H+F 
F=Q 
END IF 
END IF 
END Sir 
Rho=Rho+F 
R=R+] 
NEXT J 
Rho=Rhot4#lL+2#H 
Rho=(((Rho#.33)#K#K#1].46)°%(-.6)) 
N=(Ant_diam/Rho)*2 
IF N<=1 THEN 
Ze=1/(14N) 
EES 
Ze=1.5/(14+N) 
ISTEP OF 1.5==>SEE CLIFFORD 1980 
END IF 


Ze 


Exc_att=Ze*Ze 
ELSE 
Exc_att=Exc_att 
ENE E 
{ 
ICALCULATE THE POWER BACKSCATTERED 
Pow_back(I)=(Pou_trans*Eff-Sumpow_back )#*EXP(-Atten*Range( I ) ) 
Pow_back(I)=Pow_back(I)#*Interval#Exc_att*Sigma 
Sumpow_back =Sumpow_back+Pow_back( I) 
(CALCULATE THE POWER BACKSCATTERED TO THE ANTENNA 
Pow_ret(I)=Pow_back(I)*EXP(-Atten*Range(I))*Ant_areatGeEf f/Range(I)*Z 
PRINT “RANGE=" , Range( I) 
PRINT "POWER RETURNED=" ,Pow_ret(I) 
PRINT "FREQUENCY=" Freq 
I=I+] 
UNTIL Pow_ret(I-1)>@ AND Pow_ret(I-1)<=Noise 
Ran( Ind, Ilfreq)=Range(I-2) 
IF Ran(Ind,Ifreq)>Ymax THEN 
Ymax=Ran( Ind, lfreq) 
END IF 
IF Ran(Ind,Ifreq)<Ymin THEN 
Ymin=Ran( Ind, lfreq) 
END IF 
Ind=Ind+1 
Again=1 
FOR J=®@ TO I 
Pow_back( J)=@ 
Pow _ret(J)=@ 
Range( J)=@ 
Cne2(J)=@ 
NEXT J 
NEXT Ant_diam 
Ifreq=lfreqtl 
Ind=1 
NEXT Freq 


OUTPUT KBD; "K"; 


PRINT "INPUT CONDITIONS” 

{ 

PRINT " " 

PRINT USING "K";" TEMPERATURE= »,bemp, CELCIUS” 
PRINT USING "“K";"ATMOSPHERIC PRESSURE= ", Atom_pres, "mb" 
PRINT USING "K";"WATER VAPOR PRESSURE= " Wat_pres,"mb" 
PRINT USING "K"; "PULSE LENGTH= Tanlse, MS 
Primus mnoe Kk 3 TRANSDUCER EFFICIENCY= “Ett 

PRINT USING "“K";"POWER TRANSMITTED= " Pow_trans," WATTS" 
PRINT “TEMPERATURE STRUCTURE PROFILE USED " Profile 


IF Profile=4 THEN 


PRINT 


“INVERSION HEIGHT 


" Inver 


aa 


END IF 
SELECT Profile 


CASE 1 
PRINT ° 1 A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
PRINT " AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 " 
PRINT " WHICH GIVES A HEIGHT TO THE = 1-7 6egGR0FILE™ 

CASE 2 
PRINT " Zz A TEMPERARURE STRUCTURE PROFILE BASED ON DATA” 
PRINT “ AS PRESENTED IN WALTERS/KUNOEL 1381) PAGE 398 ” 
PRINT " WHICH GIVES A HEZGHT TO THe -TIGVPROFILE™ 
PRINT “ BUT WITH A FACTOR OF TWO TO APPROXIMATE" 
PRINT " LOOKING UP A THERMAL PLOOM" 

CASESS 
PRINT. = 3 A TEMPERATURE STRUCTURE PROFILE BASED ON DATA" 
PRINT “ AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398" 
PRINT " FOR NIGHT CONDITIONS. THE DEPENDENCE WITH” 
PRINT “ HEIGHT IS EXP(-.0@1*HEIGHI@ABOVE GSeMETERS” 
PRINT " AND HEIGHT TO THE -1.46 BELOW 65 METERS” 

CASE 4 
PRINT " 4 A TEMPERATURE STRUCTURE PROFILE BASED ON A" 
PRINT " HEIGHT 10° THES=473 5) 

END SELEGT 


! 
PRINT USING “K";"THE MINIMUM DETECTABLE SIGNAL WAS SET AT" ,Noise,” WATTS” 
| 
PRINT " ” 
| 
| 
PRINT “HIT CONTINUE TO CONTINUE" 
PAUSE 
OUTPUT KBD; "K"; 
GRAPHICS ON 
VIEWPORT 15,120,10,79 
| 
| 
'PLOT RANGE VERSUS ANTENNA DIAMETER FOR VARIOUS FREQUENCIES 
Xmin=@ 
Xmax=6 
Ymin=PROUND( Ymin-5@, 2) 
Ymax=PROUND( Ymax+5@, 2) 
Dx=1} 
Dy=10@ 
Xr ange=Xmax- Xmin 
Yrange=Ymax-Ymin 
WINDOW Xmin, Xmax , Ymin, Ymax 
AXES Dx ,Dy,Xmin,Ymin,1,1 
CLIP CEFF 
| 


NCABEL) PEGT 


}22 


earze 4.6 

LDIR @ 

LORG 5 

MOVE .5*#xXrange, Ymax#1.@5 

LABEL “RANGE VERSUS ANTENNA DIAMETER" 
1'LABEL HORIZONTAL AXES 


LDIR @ 
LORG 5 
FOR J=@ TO Xmax STEP Dx 
CSIZE 4,.6 
MOVE J,Y¥min-.@5#Yrange 
LABEL J 
NEXT J 
MOVE. .5#Xrange,Ymin-.1*Yrange 
CSI ZE 4.6 


LABEL “ANTENNA DIAMETER (m.)" 
ILABEL VERTICAL AXES 
LORG 8 
FOR J=Ymin TO Ymax STEP Dy 
CSIZE 4,.6 
MOVE Xmin- .@025* Xrange, J 
LABEL USING "K";J 
NEXT J 
EDLR IPE Z 
LORG 6 
MOVE Xmin-.15*Xrange,.5#Yranget+Ymin 
CS5lZE 4-6 
LABEL “RANGE (METERS)" 
CLIP ON 
Ind=] 
Ifreq=1 
FOR Freq=500 TO 1500 STEP 5S0® 
FOR Ant_diam=.5 TO 6.25-Ifreq STEP .25 
IF INT(I1freq/2)=Ifreq/2 THEN 
LINE TYPE 3 
EESE 
LINE TYPE 1 
END IF 
PLOT Ant_diam,Ran(Ind,I freq) 
Ind=Ind+1 
NEXT Ant_diam 
LDIR @ 
LORG 2 
CSIZE 3,...6 
MOVE Ant_diam-.1,Ran(Ind-1,I freq) 
LINE TYPE 1 
LABEL USING "K";Freq,"Hz" 
Ifreq=Ifreqt | 
Ind=1 
NEXT Freq 


ae 


PRINT “HIT CONTINUE TO CONTINUE” 
PAUSE 
GCLEAR 
OUTPUT KBD;"K"; 
| 
LINPUT “WOULD YOU LIKE TO MAKE ANOTHER RUN (CENTER Y OR N)?" XS 
Again=FNYes( XS) 
Mess_up=1 
WHILE Mess_up=1 
SELECT Again 
CASE 
LINPUT “WOULD YOU LIKE TO CHANGE A VARIABLE (ENTER Y OR N)?" | XS 
New_va=FNYes( XS) 
Ano_change=1 
WHILE Ano_change=! 
SELECT New_va 


GASE | 
PRINT “ VARIABLE CURRENT VALUE" 
PRINT USING "“K";"1 TEMPERATURE ~ hema “Cees lus” 
PRINT USING “K";"2 ATOMOSPHERIC PRESSURE " Atom_pres,'mb" 
PRINT USING "“K";"3 WATER VAPOR PRESSURE “,Wat_pres, "mb" 
PRINT " “" 
PRINT USING “K";"4 ANTENNA DIAMETER " Artec ham.” mae 
PRINT USING “K";"5 PULSE LENGTH " Pulse," ms" 
PRINT USING "K";"6 POWER TRANSMITTED " Pow _trans," WATTS" 
PRINT USING “K";"7 ATMOSPHERIC PROFILE “Profile 
PRINT “ “" 
PRINT “ENTER THE NUMBER OF THE VARIABLE YOU WISH TO" 
PRINT “CHANGE" 
INPUT Var 
SELEGieVvar 
CAsaa 
INPUT "TEMPERATURE=" , Temp 
CASeez 
INPUT “ATMOSPHERIC PRESSURE IN mb=",Atom_pres 
CASE 3 
INPUT "WATER VAPOR PRESSURE IN mb=",Wat_pres 
CASE 4 
INPUT “ANTENNA DIAMETER IN m.=",Ant_diam 
Gases = 
INPUT "PULSE LENGTH IN ms=" Pulse 
CASE 6 
INPUT “POWER TRANSMITTED IN WATTS=" ,Pow_trans 
CASE 7 
PRINT “ 1 ==> Ct*2 PROFILE OF Z*(-1.16)" 
PRINT " FROM WALTERS/KUNDEL 1981" 
PRINT “ 2 ==> SAME AS ONE BUT WITH FACTOR" 
PRINT " OF TWO TO APPROXIMATE LOOKING" 


}24 


PRINT “ UP” A THERRAL PLOOM” 
PRINE ~ 3” ==> Ct°ZePROFILE OF EXP(-.ODilmez)" 
PRINT “ FROM WALTERS/KUNDEL 1981." 
PRINT “~ 4 ==>. CT°Z PROP MBE"Or 7°(-4/3)" 
INPUT “ENTER NUMBER OF DESIRED PROFILE” ,Profile 
IF Profile=4 THEN 
INPUT “HEIGHT OF INVERSION IN METERS=", Inver 
END IF 
CASEVBESE 
PRINT Var,"IS NOT ONE OF THE OPTIONS" 
END@oELEE TS 
LINPUT "MADE ANOTHER CHANGE(C ENTER Y OR N)?", XS 
New_va=FNYes( XS) 
Mess_up=2 
CASE 22 
Mess_up=2Z 
Ano_change=Z 
CASESEESE 
PRINT “YOUR CHOICES WHERE Y OR N !!IUI" 
Ano_change= |! 
LINPUT “CHANGE A VARIABLE (ENTER Y OR N)?" XS 
New_va=FNYes( XS) 
END SELECT 
END WHILE 
EASE SZ 
Mess_up=Z 
CASE 2USE 
PRINT “YOUR CHOICES WHERE Y OR N JU!IttHtt" 
Mess_up=| 
LINPUT "MAKE ANOTHER RUN (ENTER Y OR N)?" XS 
AgainsFNYes( X$) 
END SELECT 
END WHILE 
END WHILE 
Been?) “THAT"S ALCL, FOLKS” 
end 
| 
| 
ICALCULTE THE ATTENUATION 
sub Attenuation(Atom_pres ,Atten,Freq, Temp,Wat_pres) 
'THIS SUBPROGRAM CALCULATES THE ATTENUATION OF SOUND 
'IN AIR BASED UPON EQUATIONS IN NEFF 1975 
| 
| INPUT ATMOSPHERIC PRESSURE IN MILLIBARS 
FREQUENCY OF SOUND WAVE IN HERTZ 
| TEMPERATURE IN DEGREES CELCIUS 
WARER-VAFOR PRESSURE IN MIPELIBARS 
| 


OUTPUT ATTENUATION IN 1/METERS 
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!VARIABLES 
Atom_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 
Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 
SUBPROGRAM ATTENUATION. 
Att_max VARIABLE IN SUBPROGRAM ATTENUATION. IT IS THE 
ATTENUATION AT THE FREQUENCY OF THE MAXIMUM 
ATTENUATION FOR THE INPUT CONDITIONS. 


F VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 
ATTENUATION. 

Freq INPUT FREQUENCY OF ECHOSOUNDER. 

H VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 

Pstar VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 

Temp INPUT TEMPERATURE IN DEGREES CELSIUS. 

Tstar INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 


IN SUBPROGRAM ATTENUATION. 
Wat _pres ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
| OPERATOR. 
H=1@0*#Wat_pres/Atom_pres 
Tstar=(1.8* Temp+492)/519 
Pstar=Atom_pres/1@O14 
Fmax=( 10+6600#H+444Q00#H#H)*#Pstar/Tstar*.8 
Att _max=.Q0078*#Fmax*Tstar*(-2.5)#EXP(7.77#( 1-1/Tstar) ) 
F=Freq/Fmax 
Atten=(Att_max/304.8)8((. 188F )°24+0Z24FeF/(C 1 4F*#F))°2)°.5 
Atten=(Attentl. 74E- 1O#Freq#eFreq)/4. 35 
SUBEND 
| 
def FNYes( S$) 
1'THIS FUNCTION INTERPRETS THE OPERATORS RESPONSES TO YES NO QUESTIONS 
| 
| INPUT X$ 
| 
| OUTPUT FNYes 


| 
| 
| 
| 
| 
| F max FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
l 
| 
| 
| 


!'VARIABLES 
| Temps VARIABLE STRING USED IN FUNCTION YES. 
1 X$ STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 


l RESPONSE TO YES OR NO QUESTION. 
DIM TempS( 1] 
TempS(1,1)=TRIMS( XS) 
SELECT Temps 
CASE ees 

RETURN } 
CASE No 2 ae 

RETURN 2 
CASEI 

RETURN } 
CASE VERSE 


JZ2e 


RETURN -2 
END SELECT 
FNEND 


jeg 


APEENPIXs = 
COMPUTER PROGRAM 5 


!FULLER, ROBERT PROG_S POWER 
| 
11@ SEP 8S 
| 
! 
[eeeeeee PURPOSES #4 EES 
| : 
THIS PROGRAM WILL TAKE PARAMETERS OF AN ECHOSOUNDER AND 
{ESTIMATE THE RANGE AS A FUNCTION OF POWER TRANSMITTED BY THE ECHOSOUNDER 
{FOR A RANGE OF FREQUIENCIES 
'THE FOLLOWING INPUTS ARE REQUIRED: 
| ATMOSPHERIC DATA 
1 ATOMOSPHERIC PRESSURE in millibars 
Z)SELECT ONE OF FOUR ATMOSPHERIC PROFILES FOR THE TEMPERATURE 
STRUCTURE PROFILE. 
a)FOR ONE PROFILE THE HEIGHT OF THE INVERSION LAYER IS ALSO INPUT 
S) TEMPERATURE IN DEGREESeGEESIUG 
4)WATER VAPOR PRESSURE IN millibars 


ECHOSOUNDER DATA 

5 ) FREQUENCY 

B)ANTENNA DIAMETER IN METERS 

7)PULSE LENGTH OF THE TRANSMITTED ACOUSTIC ENERGY 


'THE PROGRAM OUTPUTS THE FOLLOWING GRAPH TO AN EXTERNAL PLOTTER: 
1) RANGE(m. ) VERSUS POWER TRANSMITTED IN WATTS FOR VARIOUS FREQUENCIES 


1.902 + VARTRELESSaeer 


WISHES TO MAKE ANOUTHER RUN. 

Ano_change TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 
WISHES TO CHANGE THE VALUE OF INPUT DATA BEFORE MAKING 
ANOUTHER RUN. 

Ant_area ECHOSOUNDER ANTENNA AREA CALCULATED FROM INPUT OF 
ANTENNA DIAMETER. 

Ant_diam INPUT OF ECHOSOUNDER ANTENNA DIAMETER IN METERS. 

Atom_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 


Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 
SUBPROGRAM ATTENUATION. 
Bn BANDWIDTH OF RECEIVER NEEDED TO RECEIVE SIGNALS 


GORPLER SHIFFED BY THREE METER PER SEcone 
VERTICAL VELOCITIES. 


iC VELOCITY OF SOUND CALCULATION AT ZERO CELSIUS. 

ei VARIABLE USED TO SCALE INVERSION HEIGHT IN 
TEMPERATURE STRUCTURE PROFILE THREE. 

Be VARIABLE USED IN CALCULATING THE FOURTH Ctez2 PROFILE. 


l 
| 
| 
; 
l 
l 
l 
l 
l 
| 
| 
l 
| 
! Again TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 
| 
| 
| 
| 
l 
l 
l 
| 
l 
l 
l 
; 
l 
l 


CneZz(*) ARRAY OF VALUES OF THE ACOUSTIC REFRACTIVE INDEX 
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cam 30 comm 40mm imi te i ieeemiieieeeeeeeeeeeeteeleoemee eeeee 


CteZ 


Cve2 


Ox 

Dy 
Epsilon 
Eff 

ert 


Exc_att 
F 


F max 


Freq 
G 
H 
It 


Ifreq 
IND 


Interval 
Inver 


New_va 


Noise 


Primnnme DER. CALCULATED BASED OR? SELECTION’ OF 
PROFILE FOR TEMPERATURE STRUCTURE PARAMETER 
AND ASSUMED PROFILE FOR THE VELOCITY STRUCTURE 
PARAMETER. 
VALUE OF THE TEMPERATURE STRUCTURE PARAMETER. 
THE OPERATOR SELECTS A PROFILE FROM SEVERAL PROVIDED. 
VALUE OF TRE VELOCITY STRUCTURE PARAMETER. 
VALUES BASED ON CALCULATION USING ASSUMED DISSAPATION 
adel 
STEP ST ZEVFOR Xonxts FOR VARIOUS PLODS. 
STEP SIZE EGR Yaak ssg0n VARTOUS =FlEOms. 
DISSAPATION RATE USED IN CALCULATION OF Cve2. 
TRANSMISSION EFFICIENCY OF ECHOSOUNDER. 
EFFICIENCY OF ECHOSOUNDER WHEN RECEIVING ACOUSTIC 
BACKSCATTER. 
EXCESS “ATTENUATION” AT GIVEN RANGE. 
VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 
FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
ATTENUATION. 
INPUT FREQUENCY OF ECHOSOUNDER. 
ANTENNA EFFECTIVE APERATURE FACTOR. 
VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 
MAIN INDEX FOR LOOP TO CHECK RETURNED POWER UNTIL 
POWER WAS LESS THAN MINIMUM DETECTABLE SIGNAL. 
INDEX FOR INCREMENTS OF FREQUENCY 
INDEX USED FOR LOOP FOR DIFFERENT TRANSMITTED POWER 
MAXIMUM EFFECTIVE-SCATTERING VOLUME THICKNESS. 
HEIGHT OF INVERSION LAYER. 
FIRST ORDER INDEX FOR ASSORTED LOOPS. 
WAVENUMBER 
SECOND ORDER INDEX FOR ASSORTED LOOPS. 
TEST VALUE TO CHECK RESPONSES WHEN OPERATOR IS INPUTTING 
RESPONSES. 
TEST VALUE IN CALCULATION OF COHERENCE LENTH IN CALCULATION 
OF EXCESS ATTENUATION. 
TEST VALUE IN CHECKING IF OPERATOR WISHES TQ CHANGE 
A VARIABLE BEFORE A NEW RUN. 
ASSUMED MINIMUM DETECTABLE SIGNAL. 


Pow _back(#) POWER BACKSCATTERED FROM GIVEN RANGE. 
Pow_ret(#) POWER BACKSCATTERED TO ECHOSOUNDER FROM GIVEN RANGE. 


Pow trans 


Profile 
Pstar 


Pulse 
R 


R1 
Range( * ) 


INPUT OF POWER SUPPLIED TO TRANSDUCER OF ECHOSOUNDER. 
OPERATOR INPUT OF CteZ PROFILE FROM AVATLABLE PROFILES. 
VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 

TRANSMITTED PULSE LENGTH IN MILLISECONDS. 
OPERATOR. 

RANGE FOR INPUT TO INTEGRAL OF EXCESS ATTENUATION. 
VARIABLE USED IN SCALING INVERSION HEIGHT FOR PROFILE THREE. 
ARRAY OF RANGE VALUES. 


INPUT BY 


fas 


Remainder REMAINDER OF MODULO FUNCTION USED TO DECREES THE 
NUMBER OF PASSES THROUGH THE INTEGRATION FOR EXCESS 
ATTENUATION. 

Rge VALUE OF RANGE IN INTEGRAL FOR EXCESS ATTENUATION. 

Rho CORRELATION LENGTH USED IN CALCULATION OF EXCESS 
ATTENUATION. 

Sigma FRACTION OF ENERGY BACKSCATTERED FROM GIVEN RANGE. 


Speed_sound SPEED OF SOUND IN AIR AT INPUT TEMPERATURE. 


Sumpow_back 


| 

! 

| 

| 

! 

; 

i 

i 

i 

ee | INPUT TEMPERATURE IN DEGREES KELVIN. 

| Temp INPUT TEMPERATURE IN DEGREES CELSIUS. 

| Temps VARIABLE STRING USED IN FUNCTION YES. 

i 

! Tstar INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 
| IN SUBPROGRAM ATTENUATION. 

| Var USED TO SELECT THE VARIABLE THE OPERATOR WISHED TO 
CHANGE BEFORE MAKING ANOTHER RUN. 

| Wat pres ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
| OPERATOR. 

es THIRD ORDER INDEX USED IN VARIOUS LOOPS. 

| = X$ STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 
| RESPONSE TO YES OR NO QUESTION. 

| Xmax VALUE OF LARGEST VALUE OF X FOR VARIOUS PLOTS. 

| = =Xmin VALUE OF SMALLEST VALUE OF X FOR VARIOUS PLOTS. 

| Xrange VALUE OF RANGE OF X VALUES FOR VARIOUS PLOTS. 

! Ylabel$ LABEL ON Y AXIS PASSED TO SUBPROGRAM Semi_log FOR PLOTTING. 
! Ymax LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

! Ymin SMALLEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

| Yrange RANGE OF Y VARIABLES FOR VARIOUS PLOTS. 

JF INTERMEDIATE VALUE IN CALCULATION OF EXCESS ATTENUATION. 
| 


!'DECLARE VARIABLES 
INTEGER I 
DIM Pow_back( 1500) ,Pow_ret( 15@@) , Range( 1500) 
DIM Cne2Z( 1500) ,Ran( 50,6) 


PLOTTER IS 70S. nRGE. 
LINE Vee 


{INPUT ATMOSPHERIC DATA 

INPUT “ENTER TEMPERATURE IN DEGREES CELCIUS" , Temp 

INPUT "ENTER ATMOSPHERIC PRESSURE IN MILLIBARS" ,Atom_pres 

INPUT "ENTER WATER VAPOR PRESSURE IN MILLIBARS” ,Wat_pres 

| 

| 

! INPUT ECHOSOUNDER DATA 

INPUT “ENTER ECHOSOUNDER PULSE LENGTH IN MILLISECONDS=" ,Pulse 
INPUT “ENTER THE ECHOSOUNDER ANTENNA DIAMETER IN METERS” ,Ant_diam 
6=.40 IANTENNA EFFECTIVE APERATURE FACTOR 


ore 


Ba. 2S 


* 


SELECT THE PROFIEE FOR THE TEMPERATURE STRUCTURE 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
I 

PRINT 
PRINT 
Mess_ 
WHILE 


"YOU HAVE A CHOICE OF TEMPERATURE STRUCTURE PROFILES" 

 “PROE ILE 102] BEW USED.” 

“YOUR SELECTIONS ARE" 

i 1 A TEMPERARURE STRUCTURE PROFILE BASED ON DATA AS “ 
‘ PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 WHICH " 
i Givyesene Helen 10 THE —!.1b PROFILE’ 

i‘ Ve THE TEMPERATURE STRUCTURE PROFILE ABOVE BUT" 

ee WITH A FACTOR OF TWO TO APPROXIMATE LOOKING" 

: UP A THERMAL PLOOM" 

: 3 A TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS" 
. PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 FOR" 

. NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT " 

a IS EXP(-.@@1#HEIGHT ABOVE 65 METERS AND HEIGHT" 

4 10 THE -1.46 BELOW 65S METERS” 

2 4 A TEMPERATURE STRUCTURE PROFILE BASED ON A" 

zy HEYGHY 10 THE 473. 
up= | 

Mess_up=1 


INPUT ENTER THE DESTRES PROFILE (1 OR 2 OR 3 OR 4)” ,Profile 
IF Profile=4 THEN 

INPUT “ENTER THE HEIGHT OF THE INVERSION IN METERS", Inver 
END IF 
IF Profivbesi THEN 


Mess_up=® 
EESE 
IF Profile=Z THEN 
Mess_up=@ 
ELSE 
IF Profile=3 THEN 
Mess_up=@ 
ELSE 
IF Profile=4 THEN 
Mess_up=® 
ELSE 
PRINT Profile,” WAS NOT ONE OF THE OPTIONS! !!!" 
Mess_up=1 
END IF 
END IF 
ENDO 


END IF 


END WHILE 
OUTPUT KBD; "K"; 


Again=1 
WHILE Again=! 


FOR J=1 TO I 
Pow _back(J)=@ 
Pow_ret(J)=@ 
Range( J)=@ 
Cne2(J)=@ 
NEXT J 
| 
ICONVERT TEMPERATURE TO KELVIN 
T=Temp+273 
| 


!'CALCULATE SPEED OF SOUND BASED ON INPUT TEMPERATURE( CELCIUS) 
Speed_sound=290.05#(T)*.5 

ICALCULATE THE SPEED OF SOUND AT @ DEGREES CELCIUS 
C=70.05*273 2s 

| 


I={ [INDEX FOR INNER LOOP (POWER RETURNED TO NOISE) 

Ind=1 !INDEX FOR SECOND LOOP (EFFICIENCY) 

Ifreq=! ! INDEX FOR THIRD LOOP (FREQUENCY ) 

Pow_back(@)=@ !INITIALIZE VARIABLE FOR POWER BACK SCATERRED 


Interval=( Speed_sound#Pulse*!.E-3)/2 
Ant_area=PI*(Ant_diam/2)°2 
! 
fINITIALIZE PLOTTING PARAMETERS SO CAN COMPARE OURING RUN TO SCALE 
PEGS 
Ymin= 1000 
Ymax=@ 
ICALCULATE RETURNED SIGNAL POWER UNTIL IT IS LESS THAN BACKGROUND 
'FOR EACH FREQUENCY. 
Range( @)=@ 
FOR Freg=5@0 TO 1500 STEP 500 
FOR Pow _trans=50@ TO S@@ STEP 5@ 
Bn=2*Freq#( 1-1/(6/Speed_sound+1!))!BANOWIOTH FOR 3M/S VERTICAL VELOCITY 
Noise=1.38E-23#Bn#( T)+2.E-14! MINIMUM DETECTABLE SIGNAL AS 
| JOHNSON NOISE( NEGLIGABLE)+ ESTIMATED BACKGROUND 
ICALCULATE THE ATTENUATION COEFFICIENT BASED ON ATMOSPHERIC DATA 
CALL Attenuation(Atom_pres,Atten, Freq, Temp,Wat_pres) 
K=2*PI*Freq/Speed_sound | WAVENUMBER 
I={ 
Sumpow_back=@ 
REPEAT 
Range( I )=Range(I-1)+2 
SELECT Profile 
CASE 
'THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 


Loe 


'AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398. 
!'WAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
!{TULAROSA BASIN, NEW MEXICO. 
!AN ADDITONAL FACTOR OF 1.8 WAS INCLUDED AS PER DIFFERENT 
'AVERAGING TIME. 

Cte2=2.12*Range(I)*(-1.16) 
GAS EZ 
'THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
1AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398. 
!'WAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
!TULAROSA BASIN, NEW MEXICO. 
'AN ADDITONAL FACTOR OF 1.8 WAS INCLUDED AS PER DIFFERENT 
{AVERAGING TIME AND AN ADDITIONAL FACTOR OF TWO WAS 
!INCLUDED TO APPROXIMATE LOOKING UP A THERMAL PLOOM 

Cte2=2#2.12*Range(I)*(-1.16) 
GASE 34 
'THIS TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS" 
IPRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 FOR" 
'NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT “ 
'1TS EXP(-.@@1*#HEIGHT ABOVE 6S METERS AND HEIGHT" 
'TO THE -1.46 BELOW 6S METERS" 

IF Range(I)<6S THEN 

Cte2=75.5#Range(I)*(-1.46) 
ELSE 
CteZ=3.66E-2*EXP( -.@01#Range(I)) 

ENDS Ir 
CASE 4 
!CALCULATE THE TEMPERATURE STRUCTURE FACTOR BASED ON ASSUMPTION 
| THAT IT’S PROPORTIONAL TO HEIGHT TO THE -4/3 POWER 
| THAT THE SURFACE VERTICAL HEAT FLUX IS .@95 C.m./sec 
(THE DEPENDENCE WITH HEIGHT ABOVE THE INVERSION I PULLED FROM 
rMyeLEFT EAR 
!EQUATION FROM NEFF, 1975 
SELECT Range(I)/Inver 

CASE. “eg 
C3=((.024)#(7T)*(.667)) 
Cte2=C3#(Range(I))*¢(-1.33) 

GASE ao ue 
Cte2=Cte2 
C1=Cte2 
R1l=Range( I) 

Case. | 10 1.3 
Cte2=10*((LGTCC3*Range( I) )-LGT(C3#R1))4+LGT( C1) ) 
R2=Range( I) 

C2=Cte2 
CASE sEKSE 
Gite? =C 34 Ramge( 1)°(-1.33)=RZ2°>(-1.33))+C2 
END SELECT 
END SELECT 
ICALCULATE VELOCITY STRUCTURE FACTOR. FORMULA FROM GAYNOR 77 


Ensilon=( .2866/Range(I))#*( 1+.07#( Range(I))*(.6))*¢1.5) 
CveZ=2*Epsilon*( .667) 

ICALCULATE ACOUSTIC REFRACTIVE INDEX STRUCTURE FACTOR. 
[FORMULA FROM TETARSKI 

Cne2Z(I)=(Cte2/(2.98E+S) )+( CveZ2/(C#C)) 

ICALCULATE THE FRACTION OF POWER BACKSCATTERED FROM INTERVAL. 
!EQUATION FROM NEFF, 1975 

Sigma=( .@039*#(K*%( 1/3) )*Cte2)/(7T)%Z 

I 


ICALCULATE THE EXCESS ATTENUATION. EQUATION FROM CLIFFORD, 1980 
ITHE EXCESS ATTENUATION IS Exc_att 
IF I1<S@ THEN 
Remainder =@ 
Cisse 
IF I<15@ THEN 
Remainder=Range(I) MODULO 1@ 
EES 
Remainder=Range(I) MODULO 2® 
END IF 
END IF 
IF Remainder=@ THEN 
Rho=® 
L=Q 
H=@ 
Roe=Range( I ) 1'CONSTANT IN INTEGRAL 
R=@ 
FOR J=@ TO 2+*I 
F=Cne2(INT(J/2+1)) 
F=F #( 1-R/(Roe)*(1.67)4+¢(R/( Roe) )*(1.67)) 
IF J>@ THEN 
IF J<2*I THEN 
IF INTC J/2)=J/2 THEN 


L=L+F 
F=Q 
EGse 
H=H+F 
F=Q 
END IF 
END IF 
END IF 
Rho=RhotF 
R=R+1 
NEXT J 
Rho=Rhot4+#L+2+*H 


Rho=(((Rho#.33)*#K#*K#1.46)°%(-.6)) 
N=(Ant_diam/Rho)%2 
IF N<=1 THEN 
Ze=1/( 1+N) 
ECSE 
Ze=1.5/(14+N) 
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STEP OF |2Sea25EE CLIFRORD (digg 
END IF 
Exc_att=Ze+#Ze 
ei Losi 2 
Exc_att=Exc_att 
END IF 
| 
ICALCULATE THE POWER BACKSCATTERED 
Pow_back(I)=(Pow_trans*Eff-Sumpow_back )#EXP(-Atten*#Range(I ) ) 
Pow_back(I)=Pow_back(I)#Interval*Exc_att#Sigma 
Sumpow_back=Sumpow_back+Pow_back( I) 
ICALCULATE THE POWER BACKSCATTERED TO THE ANTENNA 
Pow _ret( I )=Pow_back(I)#EXP(-Atten*#Range(I1))*Ant_areatG#Ef f/Range(I)*2 
PRINT “RANGE=" ,Range( I) 
PRINT "POWER RETURNED=" ,Pow_ret(I ) 
PRINT “FREQUENCY=" ,Freq 
T=I+1 
UNTIL Pow_ret(I-1)>@ AND Pow_ret(I-1)<=Noise 
Ran(Ind,Ifreq)=Range(I-2) 
IF Ran(Ind,Ifreq)>Ymax THEN 
Ymax=Ran(Ind,I freq) 
END IF 
IF Ran(Ind,Ifreg)<Ymin THEN 
Ymin=Ran(Ind,I freq) 
END IF 
Ind=Indt1 
Again=1 
FOR J=@ TO I 
Pow_back(J)=@ 
Pow_ret( J)=@ 
Range( J)=@ 
Cne2(J)=@ 
NEXT J 
NEXT Pow_trans 
Ifreq=I freqt! 
Ind=1 
NEXT Freq 


{ 

OUTPUT KBD;"K"; 

PRINT “INPUT CONDITIONS" 
| 


PRINT " " 


PRINTSUSING “K°; “TEMPERATURE= letae CELCTIUS”" 
PRINT USING "“K";"ATMOSPHERIC PRESSURE= ",Atom_pres,' mb" 
PRINT USING "K";"WATER VAPOR PRESSURE= " Wat_pres, "mb" 
PRINT USING "K";"PULSE LENGTH= ",Pulse,"ms" 
PRINTSUSTNGs Ks; TRANSDUCER EFFICIENCYs “Eff 

PRINT USING "“K";"POWER TRANSMITTED= ",Pow_trans," WATTS" 


PRINTS eMeerRAlURE STRUCTURE PROFILE USED “,Profile 


Ore 


IF Profile=4 THEN 
PRINT "INVERSION HEIGHT ", Inver 


END IF 
SELECT Profile 
CASE 1 
PRINT " 1 A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
PRINT " AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 " 
PRINT " WHICH GIVES A HEIGHT TO THE -1.16 PROFILE" 
CASE 2 
PRINT “ Z A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
PRINT “ AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 " 
PRINT " WHICH GIVES A HEIGHT TO THE -1.16 PROFILE” 
PRIN BUT WITH A FACTOR OF TWO TO APPROXIMATE" 
PRINT “ LOOKING UP A THERMAL PLOOM" 
CASE 3 
PRINT " 3 A TEMPERATURE STRUCTURE PROFILE BASED ON DATA" 
PRINT " AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398" 
PRINT " FOR NIGHT CONDITIONS. THE DEPENDENCE WITH" 
PRINT “ HEIGHT IS EXP(-.@@01*HEIGHT ABOVE 65 METERS" 
PRINT AND HEIGHT TO THE -1.46 BELOW 65 METERS" 
CASE 4 
PRINT " 4 A TEMPERATURE STRUCTURE PROFILE BASED ON A" 
PRINT " HEIGHT TO THE -4/3 “™ 
END SELECT 


PRINT USING "K";"THE MINIMUM DETECTABLE SIGNAL WAS SET AT" ,Noise,"WATTS" 
! 


PRINT " " 

PRINT “HIT CONTINUE TO CONTINUE" 
PAUSE 

OUTPUT KEDe aK. ; 

! 

GRAPHICS ON 

VIEWPORT 15,1290,190,7@ 

| 

! 

I1PLOT RANGE VERSUS POWER TRANSMITTED FOR VARIOUS FREQUENCIES 
Xmin=@ 

Xmax =60@ 

Ymin=PROUND( Ymin-5@, Z) 
Ymax =PROUND( Ymaxt5@, Z) 
Dx=100 

Dy=10@ 

Xrange=Xmax-Xmin 
Yrange=Ymax-Ymin 

WINDOW Xmin, Xmax , Ymin, Ymax 
AXES Dx ,Dy,Xmin,Ymin,1,1 
CLIP GBF 
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| 
HEAeEL PLOFE 
CSI7E 4..6 
LDIR @ 
LORG 5 
MOVE .5*Xrange,Ymax#*!.@5 
LABEL “RANGE VERSUS POWER TRANSMITTED" 
|LABEL HORIZONTAL AXES 
LDIR @ 
LORG S 
FOR J=@ TO Xmax STEP Dx 
Gaize 4.26 
MOVE J,Ymin-.@5*Yrange 
LABEL J 
NEXT J 
MOVE .5*Xrange, Ymin-.1*Yrange 
SSi7E =~ .6 
LABEL “POWER TRANSMITTED (WATTS)" 
1LABEL VERTICAL AXES 
LORG 8 
FOR J=Ymin TO Ymax STEP Dy 
CSIZE 4.2.6 
MOVE Xmin-.@@025* Xrange, J 
LABEL USING "K";J 
NEXT J 
EOIR PI Z 
LORG 6 
MOVE Xmin-.15*Xrange, .5#Yranget+tYmin 
eS ZE sae > 
LABEL “RANGE ( METERS)” 
CLIP ON 
Ind=1 
Ifreq=! 
FOR Freq=50@ TO 150@ STEP 5@@ 
FOR Pow_trans=5@ TO 5S@@ STEP 5S@ 
IF INT(Ifreq/2)=Ifreq/2 THEN 
EINE TYPE™S 
EESE 
RENE TYPE 
END IF 
PLOT Pow_trans ,Ran(Ind,I freq) 
Ind=Ind+t 
NEXT Pow_trans 
LDIR @ 
LORG 2 
CSlZE 73... 6 
MOVE Pow_trans-5@,Ran(Ind-! ,I freq) 
LINE TYRE 1 
LABEL USING ok ' FPreg Hz” 
Ifreq=I freqt! 


Ind=1 
NEXT Freq 
! 
PRINT "HIT CONTINUE TO CONTINUE” 
PAUSE 
GCLEAR 
OUTPUT KBO;"K"; 
! 
LINPUT "WOULD YOU LIKE TO MAKE ANOTHER RUN CENTER Y OR N)?" XS 
Again=FNYes( X$) 
Mess_up=! 
WHILE Mess_up=1 
SELECT Again 
CASE 1 
LINPUT “WOULD YOU LIKE TO CHANGE A VARIABLE (ENTER Y OR N)?" | X$% 
New_va=FNYes( X$) 
Ano_change=1 
WHILE Ano_change=! 
SELECT New_va 


Crisey ! 
PRINT “ VARIABLE CURRENT VALUE" 
PRINT USING "K";"1 TEMPERATURE /, temp. Ceiba is” 


PRINT USING "K";"2 ATOMOSPHERIC PRESSURE ",Atom_pres,"mb" 
PRINT USING "“K";"3 WATER VAPOR PRESSURE ",Wat_pres,'mb" 


PRENT 
PRINT USING "“K";"4 ANTENNA DIAMETER " Ant_diam,”" m." 
PRINT USING °K”; "SSeGUb SE SEEiotG " ,Pulse,= Me 


PRINT USING “K"3;"6 “POWER TRANSMITTED ~ ,Poubtrarceeeworls 
PRINT USING “K"3"7 ATMOSPHERIC PROBTEE SS Ser ot ite 

PRINT " " 

PRINT "ENTER THE NUMBER OF THE VARIABLE YOU WISH TO" 

PRINT “CHANGE" 


INPUT Var 
SELECT Var 
CASE 
INPUT "TEMPERATURE=" , Temp 
CASE 
INPUT "ATMOSPHERIC PRESSURE IN mb=",Atom_pres 
Chicka 
INPUT “WATER VAPOR PRESSURE IN mb=",Wat_pres 
CASES4 
INPUT "ANTENNA DIAMETER IN m.=",Ant_diam 
CASeS 
INPUT “PULSE LENGTH IN ms=",Pulse 
CASE & 
INPUT "POWER TRANSMITTED IN WATTS=" Pow_trans 
GCASe 7 
PRINT ” 1 ==> Ct°2ePROP TBE TOF (2G - 1] ee 
PRINT " FROM WALTERS/KUNDEL 1981" 


kos 


PRINT “~ Z ==> SAME AS ONE BUT WITH FACTOR” 


PRINT = OF TWO TO APPROXIMATE LOOKING" 
PRINT “ UP A THERMAL PLOOM" 

PRINT = 3 2==> Ct°>2 PROGELE_OF BxKP( - TCI)” 
PRINT “ FROM WALTERS/KUNDEL 1981." 


PRINTS 49==>°C1°2Z PROPIEE "OF Z°(-473)" 
INPUT “ENTERPNUMBER OF DESIRED PROFILE” ,Profile 
IF Profile=4 THEN 
INPUT “HEIGHT OF INVERSION IN METERS=" ,Inver 
END IF 
CASE. Eece 
PRINT Var ,"IS NOT ONE OF THE OPTIONS" 
ENG® SELECE 
LINPUT "MADE ANOTHER CHANGEC ENTER Y OR N)?", XS 
New_va=FNYes( X$) 
Mess_up=2Z 
CASE 
Mess_up=2 
Ano_change=2 
CASE SESE 
PRINT “YOUR CHOICES WHERE Y OR N !IITI" 
Ano_change=!1 
LINPUT “CHANGE A VARIABLE (CENTER Y OR N)?" XS 
New_va=FNYes( X$) 
END SELECT 
END WHILE 
CASE sz 
Mess_up=2 
CASE ee ese 
PRINT “YOUR CHOICES WHERE Y OR N Jlliit" 
Mess_up=! 
LINPUT “MAKE ANOTHER RUN CENTER Y OR N)?" XS 
Again=FNYes( X$) 
ENDS SEESET 
END WHILE 
END WHILE 
PRINT “THAT’S ALL, FOLKS" 
end 
| 
| 
ICALCULTE THE ATTENUATION 
sub Attenuation(Atom_pres,Atten,Freq, Temp, Wat_pres) 
'THIS SUBPROGRAM CALCULATES THE ATTENUATION OF SOUND 
'IN AIR BASED UPON EQUATIONS IN NEFF 1975 
| 
| INPUT ATMOSPHERIC PRESSURE IN MILLIBARS 
| FREQUENCY OF SOUND WAVE IN HERTZ 
| Tenens rne UN DEGREES sGEEeLUs 
! WATER-VAPOR PRESSURE IN MILLIBARS 


fos 


OUTPUT ATTENUATION IN 1/METERS 
{ 


VARIABLES 
Atom_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 
Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 


SUBPROGRAM ATTENUATION. 

Att_max VARIABLE IN SUBPROGRAM ATTENUATION. IT IS THE 
ATTENUATION AT THE FREQUENCY OF THE MAXIMUM 
ATTENUATION FOR THE INPUT CONDITIONS. 


F VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 
ATTENUATION. 

Freq INPUT FREQUENCY OF ECHOSOUNDER. 

H VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 

Pstar VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 

Temp INPUT TEMPERATURE IN DEGREE S@GEL STUE. 

Tstar INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 


IN SUBPROGRAM ATTENUATION. 
Wat pres ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
OPERATOR. 
H=100*Wat_pres/Atom_pres 
Tstar=(1.8*Temp+492)/519 
Pstar=Atom_pres/1@14 
Fmax=( 10+6600#H+444Q00#H#H) #Pstar/Tstar*.8 
Att_max=.@078#Fmax#Tstar*(-2.5)#EXP(7.77#(1-1/Tstar) ) 
F=Freqg/Fmax 
Atten=(Att_max/ 304.8 )+(( 1 18*F )° 24S 2*E er / Clr) 7? oo 
Atten=(Attent!.74E-10#Freg*Freg)/4. 35 
SUBEND 
def FNYes( XS) 
' THIS FUNCTION INTERPRETS THE OPERATORS RESPONSES TO YES NO QUESTIONS 
| INPUT X$ 
| 
OUTE UT FNYes 


{ 
l 
l 
| 
| 
| 
{ 
{ Fmax FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
l 
| 
| 
| 
t 
l 


IVARTeBEES 
| TempS VARIABLE STRING USED IN FUNCTION yes= 
| XS STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 


| RESPONSE 710 YES OR NO QUESWIOGN: 
DIM TempS{ 1] 
Temp$(1,1)=TRIMS( XS) 
SELECT Temps 
CASE Yaa ve 
RETURN | 
CASE Nite 
RETURN 2 
CASE 


xe 


RETURN 1 
GAoE ELSE 

RETURN -2 
END SELECT 


FNEND 
| 
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APERNET 
COMPUTER PROGRAM 6 


1 PULEER ROBERT) PROG_6 NOISE 
\ 


110 SEP 85 

| 

‘ 

beatae tee PURPOSEs eaten 

l 

THIS PROGRAM WILL TAKE PARAMETERS OF AN ECHOSOUNDER AND 

1ESTIMATE THE RANGE AS A FUNCTION OF BACKGROUND NOISE 

!'THE FOLLOWING INPUTS ARE REQUIRED: 

ATMOSPHERIC DATA 

1 XATOMOSPHERIC PRESSURE in millibars 

Z)SELECT GNE OF FOUR AIMGSPHERICG PROFIGCES FORTTHE (EnPeEnnlURe 
STRUCTURE PROFILE. 
a)FOR ONE PROFILE THE HEIGHT OF THE INVERSION LAYER IS ALSO INPUT 

3) TEMPERATURE SIN BEGREES Carolus 

4)WATER VAPOR PRESSURE IN millibars 


ECHOSOUNDER DATA 

5 )FREQUENCY 

B)ANTENNA DIAMETER IN METERS 

7) PULSE LENGTH OF THE TRANSMITTED AC OUscITG ENERGY 
8)POWER TRANSMITTED 


THE PROGRAM OUTPUTS THE FOLLOWING GRAPH TO AN EXTERNAL PLOTTER: 
1)RANGE(m.) VERSUS BACKGROUND NOISE IN WATTS 


#2242 VARIABLES #2424 


Again TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 
WISHES TO MAKE ANOUTHER RUN. 

Ano_change TEST VALUE IN WHILE LOOP TO DETERMINE IF OPERATOR 
WISHES TO CHANGE THE VALUE OF INPUT DATA BEFORE MAKING 
ANOUTHER RUN. 

Ant_area ECHOSOUNDER ANTENNA AREA CALCULATED FROM INPUT OF 
ANTENNA DIAMETER. 

Ant_diam INPUT OF ECHOSOUNDER ANTENNA DIAMETER IN METERS. 

Atom_pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 


Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 
SUBPROGRAM ATTENUATION. 

© VELOCITY OF SOUND  CAECUEATION AleZEROTCEESIUS. 

Cl VARIABLE USED TO SCALE INVERSION HEIGHT IN 
TEMPERATURE STRUCTURE PROFILE THREE. 

Gs VARIABLE USED IN CALCULATING THE FOURTHSC leZsEROE)] EE 


CneZ(*) ARRAY OF VALUES OF THE ACOUSTIC SREPRIGiI VE sine 
PARAMETER. CALCULATED BASED ON SELEGIT ON OF 
PROFILE FOR TEMPERATURE STRUCTURE SPoaRAnelER 


! 
| 
! 
! 
| 
| 
! 
| 
! 
! 
! 
| 
| 
! 
! 
{ 
! 
! 
! 
! 
! 
| AND ASSUMED PROFILE FOR THE YEEOEI Ty siRUe TURE 
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_— comm 0h —_ — — — —_ —— cm ca —> = «tn — cm — ——<— coum — —_— ox comm — —_ c= — —_— .o—_ —_ — em —_ comm ices om comme oo fou comm vam ‘*aup coup °Gm jcows s@up oicem cGmp «m= tau 


CteZ 


CveZ 


Dx 

Dy 
Epsilon 
et 

Eff 


EX C_"att 
F 


F max 
Freq 
G 
H 
I 
IND 


Interval 
Inver 


New_va 


Noise 


PARAMETER. 
VALUE OF THE TEMPERATURE STRUCTURE PARAMETER. 
jaewORERAVOR SELECTS A PROFILE FROM SEVERAL PROVIDED: 
VALUE OF THE VELOCITY STRUCTURE PARAMETER. 
VALUES BASED ON CALCULATION USING ASSUMED DISSAPATION 
RATE. 
STEP SIZesrGR X AxIS EOR VARIOUS PLOTS. 
STEP “SIZEeOR Y AXIS FOR VARIOUS PLOTS. 
DISSAPATION RATE USED IN CALCULATION OF Cvez. 
TRANSMISSION EFFICIENCY OF ECHOSOUNDER. 
EFFICIENCY OF ECHOSOUNDER WHEN RECEIVING ACOUSTIC 
BREKSCATIER. 
EXCESS "ATTENUATION" AT GIVEN RANGE. 
VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 
FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
ATTENUATION. 
INPUT FREQUENCY OF ECHOSOUNDER. 
ANTENNA EFFECTIVE APERATURE FACTOR. 
VARIABLE USED IN INTEGRATION FOR EXCESS ATTENUATION. 
MAIN INDEX FOR LOOP TO CHECK RETURNED POWER UNTIL 
POWER WAS LESS THAN MINIMUM DETECTABLE SIGNAL. 
INDEX USED FOR LOOP FOR DIFFERENT NOISE CUT OFF LEVELS 
MAXIMUM EFFECTIVE-SCATTERING VOLUME THICKNESS. 
HEIGHT OF INVERSION LAYER. 
FIRST ORDER INDEX FOR ASSORTED LOOPS. 
WAVENUMBER 
SECOND ORDER INDEX FOR ASSORTED LOOPS. 
TEST VALUE TO CHECK RESPONSES WHEN OPERATOR IS INPUTTING 
RESPONSES. 
TEST VALUE IN CALCULATION OF COHERENCE LENTH IN CALCULATION 
OF EXCESS ATTENUATION. 
TEST VALUE IN CHECKING IF OPERATOR WISHES TO CHANGE 
A VARIABLE BEFORE A NEW RUN. 
ASSUMED MINIMUM DETECTABLE SIGNAL. 


Pow_back(*) POWER BACKSCATTERED FROM GIVEN RANGE. 
Pow_ret(#) POWER BACKSCATTERED TO ECHOSOUNDER FROM GIVEN RANGE. 


Pow_trans 


Profile 
Pstar 


Pulse 


R 
R1 
Range( # ) 


Remainder 


Rage 


INPUT OF POWER SUPPLIED TO TRANSDUCER OF ECHOSOUNDER. 

OPERATOR INPUT OF CteZ PROFILE FROM AVAILABLE PROFILES. 

VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 

TRANSMITTED PULSE LENGTH IN MILLISECONDS. 
OPERATOR. 

RANGE FOR INPUT TO INTEGRAL OF EXCESS ATTENUATION. 

VARIABLE USED IN SCALING INVERSION HEIGHT FOR PROFILE THREE. 

ARRAY OF RANGE VALUES. 

REMAINDER OF MODULO FUNCTION USED TO DECREES THE 
NUMBER OF PASSES THROUGH THE INTEGRATION FOR EXCESS 
ATTENUATION. 

VALUE OF RANGE IN INTEGRAL FOR EXCESS ATTENUATION. 


INPUT BY 
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Rho CORRELATION LENGTH USED IN CALCULATION OF EXCESS 


ATTENUATION. 

Sigma FRACTION OF ENERGY BACKSCATTERED FROM GIVEN RANGE. 

Speed_sound SPEED OF SOUND IN AIR AT INPUT TEMPERATURE. 

T INPUT TEMPERATURE IN DEGREES KELVIN. 

Temp INPUT TEMPERATURE IN DEGREES CELSIUS. 

Sumpow_back 

Temps VARIABLE STRING USED IN FUNCTION YES. 

Tstar INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 
IN SUBPROGRAM ATTENUATION. 

Var USED TO SELECT THE VARIABLE THE OPERATOR WISHED TO 


| 
| 
| 
| 
| 
| 
} 
| 
| 
| 
| 
| 
! CHANGE BEFORE MAKING ANOTHER RUN. 
| Wat pres ATMOSPHERIC WATER PRESSURE IN MILLIBARS. INPUT BY 
| 
| 
| 
| 
! 
! 
l 
! 
j 
l 
| 


COPERATOR. 

x THIRD ORDER INDEX USED IN VARIOUS LOOPS. 

XS STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 
RESPONSE TO YES OR NO QUESTION. 

Xmax VALUE OF LARGEST VALUE OF X FOR VARIOUS PLOTS. 

Xmin VALUE OF SMALLEST VALUE OF X FOR VARIOUS PLOTS. 

Xrange VALUE OF RANGE OF X VALUES FOR VARIOUS PLOTS. 

Ylabel$ LABEL ON Y AXIS PASSED TO SUBPROGRAM Semi_log FOR PLOTTING. 

Ymax LARGEST VALUE OF Y VARIABLE FOR VARIOUS PLOTS. 

Ymin SMALLEST VALUE OF Y YARTABEE FOR eVARIOUS Erors. 

Yrange RANGE OF Y VARIABLES FOR VARIOUS PLOTS. 

Le INTERMEDIATE VALUE IN CALCULATION OF EXCESS ATTENUATION. 


'DECLARE VARIABLES 

INTEGER I 

DIM Pow_back (5000) ,Pow_ret( S000) , Range‘ 5000) 

DIM CneZ{( 5000) ,Ran( SO) 

| 

PLOTTER IS 70S." HPGL™ 

ERNE tT vee 

!INPUT ATMOSPHERIC DATA 

INPUT "ENTER TEMPERATURE IN DEGREES CELCIUS" , Temp 

INPUT "ENTER ATMOSPHERIC PRESSURE IN MILLIBARS" ,Atom_pres 
INPUT "ENTER WATER VAPOR PRESSURE IN MILLIBARS" ,Wat_pres 

|! INPUT ECHOSOUNDER DATA 

INPUT "ENTER ECHOSOUNDER FREQUENCY IN HERTZ" ,Freg 

INPUT "ENTER ECHOSOUNDER PULSE LENGTH IN MILLISECONDS=" ,Pulse 
INPUT “ENTER THE ECHOSOUNDER ANTENNA DIAMETER IN METERS" ,Ant_diam 
INPUT "ENTER THE POWER TRANSMITTED IN WATTS" ,Pow_trans 

G=.40 'ANTENNA EFFECTIVE APERATURE FACTOR 


EfFf=.75 
| 
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!'SELECT THE PROFILE FOR THE TEMPERATURE STRUCTURE 


PRINT "YOU HAVE A CHOICE OF TEMPERATURE STRUCTURE PROFILES" 

PRINT "PROFILE TO BE USED." 

PRINT "YOUR SELECTIONS ARE" 

PRINT " 1 A TEMPERARURE STRUCTURE PROFILE BASED ON DATA AS “ 
PRINT “ PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 WHICH “ 
EReNT GIVES A HEIGHT TO THE -1. lies PROF IE” 

PRINT “ * 

PRINT " Z THE TEMPERATURE STRUCTURE PROFILE ABOVE BUT" 
PRINT " WITH A FACTOR OF TWO TO APPROXIMATE LOOKING" 
eelNT “ UP A THERMAL PLOOM" 

BRINT * ~ 

PRINT " 2 A TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS" 
PRENT “ PRESENTED IN WALTERS/KUNDEL 198! PAGE 398 FOR" 
PRINT “ NIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT " 
PRINT ~ TS EXP(-.@@1*HEIGHT ABOVE 6S METERS AND HEIGHT" 
PRINT "“ TOSTHE -laa6 BELOW 6S METERS” 

ERINT * ” 

PRINT " 4 A TEMPERATURE STRUCTURE PROFILE BASED ON A” 
PRINT ~ HETGHI 10, THE -473 ~ 

| 

PRINT " " 

PRINT " " 

Mess_up=! 


WHILE Mess_up=! 
INPUT "ENTER THE DESIRED PROFILE (1 
IF Profile=4 THEN 
INPUT “ENTER THE HEIGHT OF THE INVERSION IN METERS" Inver 

END IF 
IF Profile=! 

Mess_up=@ 
ELSE 

IF Profile=Z THEN 

Mess_up=@ 


GR Z90R 3 OR 4)" Profile 


THEN 


ELSE 
IF Profile=3 THEN 
Mess_up=@ 
ELSE 
IF Profile=4 THEN 
Mess_up=9 
ELSE 
PRINT Profile," WAS NOT ONE OF THE OPTIONS! !!!" 
Mess_up=! 
END IF 
END IF 
END IF 
END IF 
ENO WHILE 
OUTPUT KBD; "K"; 
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| 

{ 

Again=1 

WHILE Again=! 

FOR J=1 T0 Il 
Pow_back(J)=@ 
Pow_ret(J)=@ 
Range( J)=@ 
CneZ(J)=@ 

NEXT J 

ICONVERT TEMPERATURE TO KELVIN 

T=Temp+273 

| 


ICALCULATE SPEED OF SOUND BASED ON INPUT TEMPERATURE(CELCIUS) 
Speed_sound=2@0.@5#(T)*.5 

ICALCULATE THE SPEED OF SOUND AT @ DEGREES CELCIUS 
C=20.05*273° .5 


i=] !TINDEX FOR INNER LOOP (POWER RETURNED TO NOISE) 


Ind=! !INDEX FOR SECOND LOOP (NOISE) 
Pow_back(@)=@ lINITIALIZE VARIABLE FOR POWER BACKSCATERRED 


Interval =( Speed_sound#Pulse#l.E-3)/2 
Ant_area=PI*(Ant_diam/2Z)*2 
| 
ICALCULATE RETURNED SIGNAL POWER UNTIL IT IS LESS THAN BACKGROUND 
!FOR EACH FREQUENCY. 
Range(@)=@ 
Noise=1.E-18 
FOR Ind=1 TO 5 
CALL Attenuation(Atom_pres,Atten,Freq,Temp,Wat_pres) 
K=Z2*PI#Freq/Speed_sound | WAVENUMBER 
[=1 
REPEG 
Range( I )=Range(I-1)+2Z 
Sumpow_back =@ 
SELECT Prof me 
CASea 
{THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
'AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398. 
IWAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
I! TULAROSA BASIN, NEW MEXICO. 
1AN ADDITONAL FACTOR OF 1.8 WAS INCLUDED AS PER DIFFERENT 
!\AVERAGING TIME. 
CteZ=2.12#Range(I)*(-1.16) 
CASE 2 
'THIS TEMPERATURE STRUCTURE PROFILE IS BASED UPON DATA 
'AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398. 
IWAS FOR MIDDAY AND CLEAR WEATHER ABOVE THE 
ITULAROSA BASIN, NEW MEXICO. 
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l1AN ADDITONAL FACTOR OF 1.8 WAS INCLUDED AS PER DIFFERENT 
!AVERAGING TIME AND AN ADDITIONAL FACTOR OF TWO WAS 
!'INCLUDED TO APPROXIMATE LOOKING UP A THERMAL PLOOM 
CteZ2=2*Z2.12*Range(I)*(-1.16) 
CASES 
!THIS TEMPERATURE STRUCTURE PROFILE BASED ON DATA AS" 
IPRESENTED IN WALTERS/KUNDEL 198t PAGE 398 FOR” 
INIGHT CONDITIONS. THE DEPENDENCE WITH HEIGHT " 
11S EXP(-.@@1*#HEIGHT ABOVE 65 METERS AND HEIGHT" 
'TO THE -1.46 BELOW 6S METERS” 
IF Range(I)<65 THEN 
CteZ=75.5*Range(I)*(-1.46) 
BSE 
Cte2=3.66E-2*EXP(-.001*Range(I)) 
END IF 
CASE 4 
1CALCULATE THE TEMPERATURE STRUCTURE FACTOR BASED ON ASSUMPTION 
| THAT IT*S PROPORTIONAL TO HEIGHT TO THE -4/3 POWER 
| tHal tee sUREACE VERTECAMSREnY FEUX 15 .095 C.m./s€c 
'THE DEPENDENCE WITH HEIGHT ABOVE THE INVERSION I PULLED FROM 
MY -LERI@eAR 
('EQUATION FROM NEFF, 1975 
SELECT Range(I)/Inver 
CASEMG=2 
C3=((.024)*(T)*(.667)) 
Cte2=C3*(Range(I))*(-1.33) 
CASE .9 70 | 
Cte2=Ctez 
C1=Cte2 
Rt=Range( I) 
GASE | TOP les 
Cte2=10*((LGT(C3#Range( I) )-LGT(C3#Rt) )+LGT(Ct)) 
R2=Range( I ) 
C2=Cte2 
CASE wel oe 
CteZ=C3*(Range(I)*(-1.33)-R2*(-1.33))4+C€2 
Ene SELECT 
END SELECT 
ICALCULATE VELOCITY STRUCTURE FACTOR. FORMULA FROM GAYNOR 77 
Epsilon=(.2866/Range( I) )#( 1+.07#(Range(I))*(.6))*(1.5) 
CveZ=Z2*Epsilon*( .667) 
ICALCULATE ACOUSTIC REFRACTIVE INDEX STRUCTURE FACTOR. 
!1FORMULA FROM TETARSKI 
Cne2(1I)=(Cte2/(2.98E+5) )+( Cve2/(C#C) ) 
1CALCULATE THE FRACTION OF POWER BACKSCATTERED FROM INTERVAL. 
EQUATION FROM NEFF, 1975 
Sigma=( .0039*(K*°(1/3))#Cte2)/(T)*2 
| 


ICALCULATE THE EXCESS ATTENUATION. EQUATION FROM CLIFFORD, 198@ 
1'THE EXCESS ATTENUATION IS Exc_att 


omy 


IF I1<S@ THEN 
Remainder =@ 
ELSE 
IF I1<15@ THEN 
Remainder=Range(I) MODULO 102 
BESE 
Remainder=Range(I) MODULO 20 
ENO IF 
END IF 
IF Remainder=0 THEN 
Rho=® 
L=® 
H=@ 
Rge=Range( I ) !'CONSTANT IN INTEGRAL 
R=@ 
FOR J=@ TO Z2#I 
F=Cne2Z( INT(J/2+1)) 
F=F#( 1-R/(Rge)*(1.67)4+(R/(Rge))*(1.67)) 
IF J>@ THEN 
IF J<2*I THEN 
IF INT(J/2)=J/2 THEN 
eet 
F=@ 
ELSE 
H=H+F 
F=Q 
END IF 
END IF 
ENO IF 
Rho=Rhot+F 
R=R+1 
NEXT J 
Rho=Rhot4#L+Z2+#H 
Rho=(((Rho#.33)*K*#*K*1.46)°(-.6)) 
N=(Ant_diam/Rho)*2 
IF N<=1 THEN 
Ze=1/(14+N) 
ESE 
Ze=1.5/(14N) 
(STERSOF 1. S== > SEE ChIFFORO wil? 
ENO IF 
Exc_att=ZertZe 
ERSE 
Exc_att=Exc_att 
END IF 
| 
(CALCULATE THE POWER BACKSCATTERED 
Pow _back(I)=(Pow_trans*Eff-Sumpow_back )*EXP(-Atten*Range(I)) 
Pow _back(I)=Pow_back(I)#Interval*Exc_att#Sigma 
Sumpow_back =Sumpow_back+Pow_back( I) 
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ICALCULATE THE POWER BACKSCATTERED TO THE ANTENNA 
Pow_ret(1)=Pow_back(I)*#EXP( -Atten*Range(I))#Ant_area#G*Eff/Range(1)*Z 
PRINT "RANGE=",Range( I) 
PRINT "POWER RETURNED=" ,Pow_ret( I) 
PRINT “NOISE=",Noise 
I=I +] 
UNTIL Pow_ret(I-t1)>@ AND Pow_ret(I-1)<=Noise 
Ran( Ind) =Range( I-Z ) 
Again=1 
FOR J=@ TO I 
Pow_back(J)=@ 
Pow_ret(J)=@ 
Range( J)=@ 
CneZ( J)=0 
NEXT J 
Noise=Noiset#1@ 
NEXT Ind 
| 
| 
OUTPeUie hao. Ks 
PRINT “INPUT CONDITIONS" 
| 


PRN ieee” 


PRINT USING "K";" TEMPERATURE= peteme.. CELCIUS” 
PRINT USING "K";"ATMOSPHERIC PRESSURE= ",Atom_pres,' mb” 
PRINT USING "K";"WATER VAPOR PRESSURE= ",Wat_pres, mb" 
PRINT USING "K";"FREQUENCY= meq, waliz” 

PRINT USING "K"; "PULSE LENGTH= ",Pulse,"ms” 

PRINT USING "K";""TRANSDUCER EFFICIENCY= " Eff 

PRINT USING “K";"POWER TRANSMI TTED= "Pow _trans," WATTS" 


PRINT "TEMPERATURE STRUCTURE PROFILE USED "Profile 
IF Profile=4 THEN 
PRINT "INVERSION HEIGHT " Inver 


END IF 
SELECT Provile 
CASE | 
PRINT " 1 A TEMPERARURE STRUCTURE PROFILE BASED ON DATA” 
PRINT “ AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 " 
PRINT " UBICheGLVYES A HElGHi= +0 THE ai. 16 PROFIL ike 
CASEme 
PRINT ™ 2 A TEMPERARURE STRUCTURE PROFILE BASED ON DATA" 
PRINT " AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398 " 
PRINT " WHICH GIVES A HEIGHT TO THE -1.16 PROFILE” 
PRINT " BUT WITH A FACTOR OF TWO TO APPROXIMATE" 
PRINT " LOOKING UP A THERMAL PLOOM" 
Chsees 
PRINT " 3 A TEMPERATURE STRUCTURE PROFILE BASED ON DATA" 
PRINT " AS PRESENTED IN WALTERS/KUNDEL 1981 PAGE 398" 
PRINT " FOR NIGHT CONDITIONS. THE DEPENDENCE WITH" 
PRINT " METS OO lene lon) ABOVE 6S METERS. 
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PRINT “ AND HEIGHT TO THE -1.46 BELOW 65 METERS” 


CASE 4 
PRINT “ 4 A TEMPERATURE STRUCTURE 
PRINT “ HEIGHT TO THE -4/3 “ 
END SELECT 
PRINT " “* 


t 
t 
PRINT “HIT CONTINUE TO CONTINUE" 
PAUSE 
OUTPUT KBO;"K"; 
| 
GRAPHICS ON 
VIEWPORT 15,120,10,7@ 
I'PLOT RANGE VERSUS NOISE 
Xmin=-18 
Xmax=-13 
Ymax=4 
Ymin=2 
Dx=1 
Dy=i 
Yrange=Z 
Xrange=5 
WINDOW Xmin, Xmax , Ymin, Ymax 
AXES Dx ,Dy,Xmin,Ymin,t, | 
CLIP OFF 
| 
| CABELS Pion 
CSI ZE 4326 
LDIR @ 
LORG 4 
MOVE Xmin+.5*xXrange, Ymint+1i*Yrange 
LABEL “RANGE vs NOISE” 
!LABEL HORIZONTAL AXES 
LDIR @ 
ECRGzS 
FOR J=Xmin TO Xmax STEP Dx 
CSIZE 4236 
MOVE J-.@13*#Xrange, Ymin- .@67*Yrange 
LABEL ig 
MOVE J+.013*xXrange, Ymin- .0@33* Yrange 
CSI ZEs2Z 
LABEL J 
NEXT J 
MOVE Xmin+.5*Xrange, Ymin-.12*Yrange 
LORG 5S 
CSIZE 45 


oe 


PROFILE BASED ON A" 


LABEL “NOISE (WATTS)” 
!LABEL VERTICAL AXES 
LORG 8 
FOR J=Ymin TO Ymax STEP Dy 
CSIZE 4,.6 
MOVE Xmin-.@14*xXrange,J 
LABEL 190 
CSE 7 
MOVE Xmin- .9025*xXrange,J+.93*Yrange 
LABEL J 
NEXT J 
POrR Pl7Z 
LORG 6 
MOVE Xmin-.1*Xrange, Ymint+.5*Yrange 
CSIZE 4,.6 
LABEL "RANGE (METERS)” 
CLIP ON 
Ind=1 
Noise=1.E-18 
FOR Ind=1 TO 5 
PLOT LGT( Noise) ,LGT(Ran( Ind) ) 
Noise=Noise#1Q 
NEXT Ind 
PRINT “HIT CONTINUE TO CONTINUE” 
PAUSE 
GCLEAR 
OUTPUT KBD: °K"; 
! 
; 
LINPUT "WOULD YOU LIKE TO MAKE ANOTHER RUN (ENTER Y OR N)?" XS 
Again=FNYes( X$) 
Mess_up=]1 
WHILE Mess_up=] 
SELECT Again 
CASE 1] 

LINPUT "WOULD YOU LIKE TO CHANGE A VARIABLE (ENTER Y OR N)?”" XS 
New_va=FNYes( X$) 
Ano_change=! 

WHILE Ano_change=! 
SELECT New_va 


CASES) 
PRINT VARIABLE CURRENT VALUE" 
PRINT USING "K"3;"1 TEMPERATURE ", temp, CELSIUS: 


PRINT USING "K"3;"2 ATOMOSPHERIC PRESSURE ° ,Atom_pres,''mb" 
PRINT USING “K"3;"3 WATER VAPOR PRESSURE ",Wat_pres, mb" 


earn ~~ 
PRINT USING “K";"4 ANTENNA DIAMETER  AnkeGtam. Mm. 
Grin USING =K°3"S PULSE LENGTH " Pulee Ne 


PRINT USING "K";"6 POWER TRANSMITTED " ,Pow_trans,”" WATTS" 


So | 


PRINT 
PRINT 
PRINT 


USING "K";"% ATMOSPHERIC PROFILE 


“ENTER THE NUMBER OF THE VARIABLE YOU WISH TO" 


“ ProTrire 


PRINT “CHANGE" 
INPUT Var 
SELECT Var 
CASE! 
INPUT “TEMPERATURE=" , Temp 
CASeae- 
INPUT “ATMOSPHERIC PRESSURE IN mb="Atom_pres 
CASE eo 
INPUT “WATER VAPOR PRESSURE IN mb=",Wat_pres 
CASE 4 
INPUT “ANTENNA DIAMETER IN m.=",Ant_diam 
CASE 5S 
INPUT “PULSE LENGTH IN ms=",Pulse 
CASE, b 
INPUT “POWER TRANSMITTED IN WATTS=" ,Pow_trans 
CASE 7 
PRINT " | @&> Gt OZeEROFITEE JOR. Jope 
PRENT = FROM WALTERS/KUNDEL 1981" 
PRINT " 2 ==> SAME AS ONE BUT WITH FACTOR" 
PRINT " OF TWO TO APPROXIMATE LOOKING” 
PRINT “ UP A THERMAL PLOOM" 
PRINT * 3 ==> Ct°Z PROFILE OF EXP(-.O0l22). 
PRINT " FROM WALTERS/KUNDEL 1981." 
PRINT “ 4 ==> Gi Z FROEITES OPez. 4 4 


INPUT “ENTER NUMBER OF DESIRED PROFILE” ,Profile 
IF Profile=4 THEN 
INPUT "HEIGHT OF INVERSION IN METERS=" Inver 
ENO IF 
CASEREL SE 
PRINT Var,"IS NOT ONE OF THE OPTIONS" 
ENO SELECT 
LINPUT “MADE ANOTHER CHANGECENTER Y OR N)?" | XS 
New_va=FNYes( X$) 
Mess_up=2 
CASE Z 
Mess_up=2 
Ano_change=2 
EIESE 
PRINT "YOUR CHOICES WHERE Y OR N 
Ano_change=]1 
LINPUT "CHANGE A VARIABLE (ENTER Y OR N)?", XS 
New_va=FNYes( XS) 
END (SEES I 
END WHILE 
CASE =a 
Mess_up=2 
CASE ese 


CASE 
ee 


-_ *£ s @# 8 


pz 


PRINT “YOUR CHOICES WHERE Y OR N Jlliii” 
Mess_up=] 
LINPUT "MAKE ANOTHER RUN CENTER Y OR N)?", XS 
Again=FNYes( XS) 
END SELECT 
END WHILE 
END WHILE 
PRINT “THAT’S ALL, FOLKS" 
end 
| 
| 
ICALCULTE THE ATTENUATION 
sub Attenuation(Atom_pres,Atten,Freq,Temp,Wat_pres) 
1'THIS SUBPROGRAM CALCULATES THE ATTENUATION OF SOUND 
1IN AIR BASED UPON EQUATIONS IN NEFF 1975 
{ 
INPUT ATMOSPHERIC PRESSURE IN MILLIBARS 
FREQUENCY OF SOUND WAVE IN HERTZ 
TEMPERATURE IN DEGREES CELCIUS 
WATER-VAPOR PRESSURE IN MILLIBARS 


OUTPUT ATTENUATION IN 1!/METERS 


!'VARIABLES 
Atom pres INPUT OF ATMOSPHERIC PRESSURE IN mb. 
Atten ATTENUATION OF ACOUSTIC WAVE. CALCULATED IN 


| 

I 

| 

| 

| 

| 

l 

| 

i 

| SUBPROGRAM ATTENUATION. . 

| Att _max VARIABLE IN SUBPROGRAM ATTENUATION. IT IS THE 
| ATTENUATION AT THE FREQUENCY OF THE MAXIMUM 
| ATTENUATION FOR THE INPUT CONDITIONS. 

l 
| 
| 
| 
| 
I 
| 
| 
| 
I 
I 
| 


E VARIABLE USED IN SUBPROGRAM ATTENUATION. IS THE 
RATIO OF FREQUENCY TO FREQUENCY AT MAXIMUM ATTENUATION. 

F max FREQUENCY OF MAXIMUM ATTENUATION. USED IN SUBPROGRAM 
ATTENUATION. 

Freq INPUT FREQUENCY OF ECHOSOUNDER. 

H VARIABLE USED IN INTECRATION FOR EXCESS ATTENUATION. 

Pstar VARIABLE USED IN SUBPROGRAM ATTENUATION FOR INTERMEDIATE 
CALCULATION. 

Temp POEMS lemPeRATURE IN DEGREES CELSIUS. 

Tstar INTERMEDIATE VALUE USED IN CALCULATION OF ATTENUATION 


IN SUBPROGRAM ATTENUATION. 
Wat_pres Alin@seeerre WATER PRESSURE IN MILLIBARS. INPUT BY 
OPERATOR. 
H=10@0*#Wat_pres/Atom_pres 
Tstar=(1.88 Tempi492)/519 
Pstar=Atom_pres/1Q014 
Fmax=( 1046600*H+44400+#H#H)*#Pstar/Tstar*.8 
Att_max=.Q0Q0784#Fmax#Tstar*(-2.S)#EXP(7.77#(1-3/Tstar) ) 
F=Freq/Fmax 
Atten=(Att_max/304.8)#((.188F)*2+(28FeF/C1+FeF 720° .5 


Pore 


Atten=(Attent+!.74E-10*#Freq#Freq)/4.35 
SUBEND 

! 

def FNYes( XS) 


!'THIS FUNCTION INTERPRETS THE OPERATORS RESPONSES TO YES NO QUESTIONS 


INPUT X$ 

! 

! OUTPUT FNYes 

!'VARIABLES 

! Temps VARIABLE STRING USED IN FUNCTION YES. 

| XS STRING PASSED TO FUNCTION YES TO DETERMINE OPERATOR 


RESPONSE TO YES OR NO QUESTION. 
DIM Temp$[ 1] 
Temp${1,1]=TRIMS( XS) 
SELECT Temp$ 
CASE] ioe 

RETURN 1 
CASE “Ne Sane 

RETURN 2 
Cacemaue 

RETURN 1 
CASE ELSE 

RETURN -2 
END SELECT 

FNEND 
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